1 #include "isajet/pilot.h"
2 SUBROUTINE SIGINT(F,Z,A1S,B1S,A2S,B2S)
4 C F(N+1) = INT(-Z,Z)(DX X**N/((A1+B1*X)*A2+B2*X)))
5 C F(8) = F(9) = 0 (DUMMY VALUES)
8 #if defined(CERNLIB_DOUBLE)
9 DOUBLE PRECISION A1,B1,A2,B2,A,B,C,Z,F,A1S,B1S,A2S,B2S
19 C SPECIAL CASE: X**N/(A1*A2)
20 IF(B1.EQ.0..AND.B2.EQ.0.) THEN
23 F(3)=2.*Z**3/(3.*A1*A2)
25 F(5)=2.*Z**5/(5.*A1*A2)
27 F(7)=2.*Z**7/(7.*A1*A2)
31 C SPECIAL CASE: X**N/(A+BX)
32 IF(B1.EQ.0..OR.B2.EQ.0.) THEN
42 F(2)=-A*LOG((A+Z)/(A-Z))+2.*Z
44 F(3)=A**2*LOG((A+Z)/(A-Z))-2.*A*Z
46 F(4)=-A**3*LOG((A+Z)/(A-Z))+2.*A**2*Z+2.*Z**3/3.
48 F(5)=A**4*LOG((A+Z)/(A-Z))-2.*A**3*Z-2.*A*Z**3/3.
50 F(6)=-A**5*LOG((A+Z)/(A-Z))+2.*A**4*Z+2.*A**2*Z**3/3.+2.*Z**5/5.
52 F(7)=A**6*LOG((A+Z)/(A-Z))-2.*A**5*Z-2.*A**3*Z**3/3.
63 C SPECIAL CASE: X**N/(A+B*X)**2
68 F(2)=-2.*A*Z/(A**2-Z**2)+LOG((A+Z)/(A-Z))
70 F(3)=(4.*A**2*Z-2.*Z**3)/(A**2-Z**2)-2.*A*LOG((A+Z)/(A-Z))
72 F(4)=(4.*A*Z**3-6.*A**3*Z)/(A**2-Z**2)+3.*A**2*LOG((A+Z)/(A-Z))
74 F(5)=(-16.*A**2*Z**3/3.+8.*A**4*Z-2.*Z**5/3.)/(A**2-Z**2)
75 $ -4.*A**3*LOG((A+Z)/(A-Z))
77 F(6)=(4.*A*Z**5/3.+20.*A**3*Z**3/3.-10.*A**5*Z)/(A**2-Z**2)
78 $ +5*A**4*LOG((A+Z)/(A-Z))
80 F(7)=(-8.*A**2*Z**5/5.-8.*A**4*Z**3+12.*A**6*Z-2.*Z**7/5.)
81 $ /(A**2-Z**2)-6.*A**5*LOG((A+Z)/(A-Z))
87 F(1)=(-LOG((A1+Z)/(A1-Z))+LOG((A2+Z)/(A2-Z)))/(A1-A2)
89 F(2)=(A1*LOG((A1+Z)/(A1-Z))-A2*LOG((A2+Z)/(A2-Z)))/(A1-A2)
91 F(3)=(-A1**2*LOG((A1+Z)/(A1-Z))+A2**2*LOG((A2+Z)/(A2-Z)))/(A1-A2)
94 F(4)=(A1**3*LOG((A1+Z)/(A1-Z))-A2**3*LOG((A2+Z)/(A2-Z)))/(A1-A2)
97 F(5)=(-A1**4*LOG((A1+Z)/(A1-Z))+A2**4*LOG((A2+Z)/(A2-Z)))/(A1-A2)
98 $+2.*Z*(A1*A2+A1**2+A2**2)+2.*Z**3/3.
100 F(6)=(A1**5*LOG((A1+Z)/(A1-Z))-A2**5*LOG((A2+Z)/(A2-Z)))/(A1-A2)
101 $+2.*Z*(-A1*A2**2-A1**2*A2-A1**3-A2**3)+2.*Z**3/3.*(-A1-A2)
103 F(7)=(-A1**6*LOG((A1+Z)/(A1-Z))+A2**6*LOG((A2+Z)/(A2-Z)))/(A1-A2)
104 $+2.*Z*(A1*A2**3+A1**2*A2**2+A1**3*A2+A1**4+A2**4)
105 $+2.*Z**3/3.*(A1*A2+A1**2+A2**2)+2.*Z**5/5.