5 IMPLICIT REAL*8(a-
h,o-
z)
14 s=dlog(dlog(qq/dlam**2)/dlog(q0**2/dlam**2))
17 at1=0.419+0.004*
s-0.007*
s**2
18 at2=3.460+0.724*
s-0.066*
s**2
19 gmud=4.40-4.86*
s+1.33*
s**2
20 at3=0.763-0.237*
s+0.026*
s**2
21 at4=4.00+0.627*
s-0.019*
s**2
22 gmd=-0.421*
s+0.033*
s**2
24 cas=1.265-1.132*
s+0.293*
s**2
25 as=-0.372*
s-0.029*
s**2
26 bs=8.05+1.59*
s-0.153*
s**2
27 aphs=6.31*
s-0.273*
s**2
28 btas=-10.5*
s-3.17*
s**2
38 cag=1.56-1.71*
s+0.638*
s**2
39 ag=-0.949*
s+0.325*
s**2
40 bg=6.0+1.44*
s-1.05*
s**2
41 aphg=9.0-7.19*
s+0.255*
s**2
42 btag=-16.5*
s+10.9*
s**2
47 at2=3.33+0.753*
s-0.076*
s**2
48 gmud=6.03-6.22*
s+1.56*
s**2
49 at3=0.761-0.232*
s+0.023*
s**2
50 at4=3.83+0.627*
s-0.019*
s**2
51 gmd=-0.418*
s+0.036*
s**2
53 cas=1.67-1.92*
s+0.582*
s**2
54 as=-0.273*
s-0.164*
s**2
55 bs=9.15+0.530*
s-0.763*
s**2
57 btas=-101.0*
s+44.7*
s**2
58 gms=223.0*
s-117.0*
s**2
67 cag=0.879-0.971*
s+0.434*
s**2
69 bg=4.0+1.23*
s-0.254*
s**2
70 aphg=9.0-5.64*
s-0.817*
s**2
71 btag=-7.54*
s+5.50*
s**2
72 gmg=-0.596*
s+1.26*
s**2
74 300 b12=dexp(
gmre(at1)+
gmre(at2+1.d0)-
gmre(at1+at2+1.d0))
76 cnud=3.d0/b12/(1.d0+gmud*at1/(at1+at2+1.d0))
77 cnd=1.d0/b34/(1.d0+gmd*at3/(at3+at4+1.d0))
82 fud1=cnud*x1**at1*(1.d0-x1)**at2*(1.d0+gmud*x1)
83 fs1=cas*x1**as*(1.d0-x1)**bs*(1.d0+aphs*x1
84 & +btas*x1**2+gms*x1**3)
85 f(1,3)=cnd*x1**at3*(1.d0-x1)**at4*(1.d0+gmd*x1)+fs1/6.d0
86 f(1,1)=fud1-
f(1,3)+fs1/3.d0
91 f(1,7)=cag*x1**ag*(1.d0-x1)**bg*(1.d0+aphg*x1
92 & +btag*x1**2+gmg*x1**3)
94 fud2=cnud*x2**at1*(1.d0-x2)**at2*(1.d0+gmud*x2)
95 fs2=cas*x2**as*(1.d0-x2)**bs*(1.d0+aphs*x2
96 & +btas*x2**2+gms*x2**3)
97 f(2,3)=cnd*x2**at3*(1.d0-x2)**at4*(1.d0+gmd*x2)+fs2/6.d0
98 f(2,1)=fud2-
f(2,3)+fs2/3.d0
103 f(2,7)=cag*x2**ag*(1.d0-x2)**bg*(1.d0+aphg*x2
104 & +btag*x2**2+gmg*x2**3)
108 aax=1.193*alog(float(
ihnt2(1)))**0.16666666
109 rrx=aax*(x1**3-1.2*x1**2+0.21*x1)+1.0
110 & +1.079*(float(
ihnt2(1))**0.33333333-1.0)
111 & /dlog(
ihnt2(1)+1.0d0)*dsqrt(x1)*dexp(-x1**2/0.01)
112 IF(ip_crs.EQ.1 .OR.ip_crs.EQ.3) rrx=dexp(-x1**2/0.01)
118 aax=1.193*alog(float(
ihnt2(3)))**0.16666666
119 rrx=aax*(x2**3-1.2*x2**2+0.21*x2)+1.0
120 & +1.079*(float(
ihnt2(3))**0.33333-1.0)
121 & /dlog(
ihnt2(3)+1.0d0)*dsqrt(x2)*dexp(-x2**2/0.01)
122 IF(ip_crs.EQ.2 .OR. ip_crs.EQ.3) rrx=dexp(-x2**2/0.01)