8 common/lujets/
n,
k(9000,5),
p(9000,5),
v(9000,5)
17 p1=sqrt(e1**2-
p(jl,5)**2)
18 p3=sqrt(e3**2-
p(jl+1,5)**2)
20 IF(p1.GT.0..AND.p3.GT.0.) cbet=(
p(jl,5)**2
21 & +
p(jl+1,5)**2+2.*e1*e3-
s*(1.-x2))/(2.*p1*p3)
22 IF(abs(cbet).GT.1.0) cbet=
max(-1.,min(1.,cbet))
27 psi=.5*
ulangl(p1**2+p3**2*cos(2.*bet),-p3**2*sin(2.*bet))
32 ELSE IF(p3.GT.p1)
THEN
33 psi=.5*
ulangl(p3**2+p1**2*cos(2.*bet),-p1**2*sin(2.*bet))
46 p(jl+2,1)=pt3*sin(del)
47 p(jl+2,2)=-pt3*cos(del)
50 p(jl+1,1)=-
p(jl,1)-
p(jl+2,1)
51 p(jl+1,2)=-
p(jl,2)-
p(jl+2,2)
52 p(jl+1,3)=-
p(jl,3)-
p(jl+2,3)