9 common/ludat1/mstu(200),paru(200),mstj(200),parj(200)
15 f3re=
pyhispen((ga-1.)/(ga+be-1.),0.,1)-
18 & (
log(be)**2-
log(be-1.)**2)/2.+
log(ga)*
log((ga+be-1.)/be)+
19 &
log(ga-1.)*
log((be-1.)/(ga+be-1.))
21 ELSEIF(
eps.LT.1.)
THEN
22 f3re=
pyhispen((ga-1.)/(ga+be-1.),0.,1)-
25 &
log(ga/(1.-ga))*
log((ga+be-1.)/(be-ga))
26 f3im=-paru(1)*
log((ga+be-1.)/(be-ga))
28 rsq=
eps/(
eps-1.+(2.*be-1.)**2)
29 rcthe=rsq*(1.-2.*be/
eps)
30 rsthe=
sqrt(rsq-rcthe**2)
31 rcphi=rsq*(1.+2.*(be-1.)/
eps)
32 rsphi=
sqrt(rsq-rcphi**2)
38 & (
phi-the)*(
phi+the-paru(1))
43 IF(ireim.EQ.1)
pyhii3au=2./(2.*be-1.)*f3re
44 IF(ireim.EQ.2)
pyhii3au=2./(2.*be-1.)*f3im