4 SUBROUTINE luradk(ECM,MK,PAK,THEK,PHIK,ALPK)
7 common/ludat1/mstu(200),paru(200),mstj(200),parj(200)
12 &parj(163)*
log((
xx-szm)**2+szw**2)+parj(164)*atan((
xx-szm)/szw)
17 IF(parj(160).LT.
rlu(0))
RETURN
22 xku=
min(parj(136),1.-(2.*parj(127)/ecm)**2)
23 IF(mstj(102).LE.1)
THEN
24 100 xk=1./(1.+(1./xkl-1.)*((1./xku-1.)/(1./xkl-1.))**
rlu(0))
25 IF(1.+(1.-xk)**2.LT.2.*
rlu(0)) goto 100
29 szm=1.-(parj(123)/ecm)**2
30 szw=parj(123)*parj(124)/ecm**2
34 fxkr=fxkl+
rlu(0)*(fxku-fxkl)
46 IF(nxk.LT.15.AND.fxku-fxkl.GT.fxkd) goto 110
47 xk=xkl+(xku-xkl)*(fxkr-fxkl)/(fxku-fxkl)
53 120 cthm=pme*(2./pme)**
rlu(0)
54 IF(1.-(xk**2*cthm*(1.-0.5*cthm)+2.*(1.-xk)*pme/
max(pme,
55 &cthm*(1.-0.5*cthm)))/(1.+(1.-xk)**2).LT.
rlu(0)) goto 120
57 IF(
rlu(0).GT.0.5) cthe=-cthe
58 sthe=
sqrt(
max(0.,(cthm-pme)*(2.-cthm)))
64 IF(0.5*(2.-xk*(1.-cthe))**2/((2.-xk)**2+(xk*cthe)**2).GT.
66 alpk=asin(sgn*sthe*(xk-sgn*(2.*
sqrt(1.-xk)-2.+xk)*cthe)/
67 &(2.-xk*(1.-sgn*cthe)))