8 common/ludat1/mstu(200),paru(200),mstj(200),parj(200)
13 & 1.000000
e+00, -5.000000
e-01, 1.666667
e-01,
14 & 0.000000
e+00, -3.333333
e-02, 0.000000
e+00,
15 & 2.380952
e-02, 0.000000
e+00, -3.333333
e-02,
16 & 0.000000
e+00, 7.575757
e-02, 0.000000
e+00,
17 &-2.531135
e-01, 0.000000
e+00, 1.166667
e+00/
21 IF(abs(1.-xre).LT.1.
e-6.AND.abs(xim).LT.1.
e-6)
THEN
22 IF(ireim.EQ.1)
pyhispen=paru(1)**2/6.
27 xmod=sqrt(xre**2+xim**2)
28 IF(xmod.LT.1.
e-6)
THEN
34 xarg=
sign(acos(xre/xmod),xim)
40 algxim=xarg-
sign(paru(1),xarg)
41 sp0re=-paru(1)**2/6.-(algxre**2-algxim**2)/2.
54 xmod=sqrt(xre**2+xim**2)
55 xarg=
sign(acos(xre/xmod),xim)
58 sp0re=sp0re+sgn*(paru(1)**2/6.-(algxre*algyre-algxim*algyim))
59 sp0im=sp0im-sgn*(algxre*algyim+algxim*algyre)
65 xmod=sqrt(xre**2+xim**2)
66 xarg=
sign(acos(xre/xmod),xim)
75 termre=(savere*zre-saveim*zim)/float(
i+1)
76 termim=(savere*zim+saveim*zre)/float(
i+1)
80 100 spim=spim+b(
i)*termim
82 IF(ireim.EQ.1)
pyhispen=sp0re+sgn*spre
83 IF(ireim.EQ.2)
pyhispen=sp0im+sgn*spim