VBF: coefmom.f

File coefmom.f, 72.0 KB (added by trac, 7 years ago)
Line 
1
2      PROGRAM MAIN
3C THIS PROGRAM COMPUTES THE MOMENTS OF THE COEFFICIENT FUNCTIONS
4C THE COEFFICIENT FUNCTIONS ARE PRESENTED IN THE DIS AND THE MSBAR SCHEME .
5
6      IMPLICIT REAL*8(A-H,O-Z)
7      DOUBLE PRECISION LAMBDA,NF
8      COMMON/PAR/NF,Q2,FS,RS,DN
9      COMMON/PAR1/LAMBDA,Q02
10      COMMON/DIS/NDIS
11      COMMON/CHOICE/MODE
12      COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
13      EXTERNAL F30R,F31R,F31S,F32R,F32S
14      EXTERNAL F20R,F21R,F21S,F22R,F22S
15      EXTERNAL F1G1,F1G2,F1Q2
16csm     
17      external FG21R,FG22R,FPS22R,FL1R,FL2R,FPSL2R,FLG1R,FLG2R
18
19C CHOOSE THE MOMENT DN=N  M^(N) = \int_0^1 X^N C_{i,k}(X,Q^2/\MU^2)
20
21      write(*,*)
22     # 'enter number of flavours nf'
23      read (*,*) NF
24      write(*,*)
25     # 'enter Mellin moment N'
26      read (*,*) DNNN
27
28      DN = DNNN - 1.d0
29csm      DN=5.0D0
30
31C CHOOSE FA=1.0D0 IF ONE TAKES THE ZEROTH MOMENT DN=0.0D0
32C CHOOSE FA=2.0D0 IF ONE TAKES DN > 0.0D0
33
34      FA=2.0D0
35      KEY=2
36      EPS=1.0D-10
37      MAX=10000
38C HERE FOLLOW THE COLOUR FACTORS CF, CF^2, CA.CF, CF.T_f
39
40      CF=4.0D0/3.0D0
41C     CF=0.0D0
42      CA=3.0D0
43C     CA=0.0D0
44      CF2=16.0D0/9.0D0
45C     CF2=0.0D0
46      CACF=4.0D0
47C     CACF=0.0D0
48      CFTF=2.0D0/3.0D0
49C     CFTF=0.0D0
50
51C CHOOSE DIS OR DISP SCHEME
52C DISP SCHEME IS A SCHEME IN WHICH THE MOMENTUM CONSERVATION SUM RULE
53C IS SATISFIED
54
55C DISP SCHEME  : NDIS=1
56C DIS SCHEME   : NDIS=2
57
58      NDIS=0
59
60      PI=3.14159265359D0
61
62C-------------------------------------------------------------------
63C HERE FOLLOWS THE VIRTUAL PHOTON MOMENTUM IN GEV
64      Q=DSQRT(10.0D0)
65csm      WRITE(*,21) Q
66   21 FORMAT(' ','THE VIRTUAL PHOTON MOMENTUM=',D20.10)
67      Q2=Q*Q
68C HERE FOLLOWS THE FACTORIZATION SCALE IN GEV \MU=FS
69      FS=Q
70csm      WRITE(*,23) FS
71   23 FORMAT(' ','THE FACTORIZATION SCALE=',D20.10)
72C HERE FOLLOWS THE RENORMALIZATION SCALE IN GEV
73      RS=Q
74csm      WRITE(*,24) RS
75   24 FORMAT(' ','THE RENORMALIZATION SCALE=',D20.10)
76
77C HERE FOLLOWS THE NUMBER OF FLAVOURS
78csm      NF=4.0D0
79      WRITE(*,25) NF
80   25 FORMAT(' ','nf  = ',g20.2)
81
82      DNN = DN + 1.0d0
83      WRITE(*,325) DNN
84  325 FORMAT(' ','N   = ',g20.2)
85
86      R20R=DAIND(0.0D0,1.0D0,F20R,EPS,KEY,MAX,KOUNT,EST)
87      WRITE(*,29) R20R
88   29 FORMAT(' ','c2qq0   = ',D20.10)
89      R21R=DAIND(0.0D0,1.0D0,F21R,EPS,KEY,MAX,KOUNT,EST)
90      IF(FA.GT.1.1D0) THEN
91      R21S=DAIND(0.0D0,1.0D0,F21S,EPS,KEY,MAX,KOUNT,EST)
92      ELSE
93      R21S=0.0D0
94      ENDIF
95      F2QR1=R21R
96      F2QL1=R21S
97      F2QVS1=COR1(0.0D0)
98      F2Q1=F2QR1+F2QL1+F2QVS1
99      WRITE(*,33) F2Q1
100   33 FORMAT(' ','c2qq1   = ',D20.10)
101csm beg added
102      RG21R=DAIND(0.0D0,1.0D0,FG21R,EPS,KEY,MAX,KOUNT,EST)
103      RG21S=0.0D0
104      F2GR1=RG21R
105      F2GL1=RG21S
106      F2G1=F2GR1+F2GL1
107      WRITE(*,333) F2G1
108  333 FORMAT(' ','c2gg1   = ',D20.10)
109csm end added
110   35 FORMAT(' ','c3qq1   = ',D20.10)
111      R22R=DAIND(0.0D0,1.0D0,F22R,EPS,KEY,MAX,KOUNT,EST)
112      IF(FA.GT.1.1D0) THEN
113      R22S=DAIND(0.0D0,1.0D0,F22S,EPS,KEY,MAX,KOUNT,EST)
114      ELSE
115      R22S=0.0D0
116      ENDIF
117      F2QR2=R22R
118      F2QL2=R22S
119      F2QVS2=COR2(0.0D0)
120      F2Q2=F2QR2+F2QL2+F2QVS2
121csm      WRITE(*,*) F2QR2,F2QL2, F2QVS2
122      WRITE(*,37) F2Q2
123   37 FORMAT(' ','c2qq2   = ',D20.10)
124csm beg added
125      RPS22R=DAIND(0.0D0,1.0D0,FPS22R,EPS,KEY,MAX,KOUNT,EST)
126      RPS22S=0.0D0
127      F2PSR2=RPS22R
128      F2PSL2=RPS22S
129      F2PS2=F2PSR2+F2PSL2
130      WRITE(*,337) F2PS2
131 337  FORMAT(' ','c2qqps2 = ',D20.10)
132      RG22R=DAIND(0.0D0,1.0D0,FG22R,EPS,KEY,MAX,KOUNT,EST)
133      RG22S=0.0D0
134      F2GR2=RG22R
135      F2GL2=RG22S
136      F2G2=F2GR2+F2GL2
137      WRITE(*,338) F2G2
138 338  FORMAT(' ','c2gg2   = ',D20.10)
139
140      RL1R=DAIND(0.0D0,1.0D0,FL1R,EPS,KEY,MAX,KOUNT,EST)
141      FLQR1=RL1R
142      FLQL1=0.d0
143      FLQ1=FLQR1+FLQL1
144      WRITE(*,331) FLQ1
145 331  FORMAT(' ','cLqq1   = ',D20.10)
146
147      RL2R=DAIND(0.0D0,1.0D0,FL2R,EPS,KEY,MAX,KOUNT,EST)
148      FLQR2=RL2R
149      FLQL2=0.d0
150      FLQ2=FLQR2+FLQL2
151      WRITE(*,332) FLQ2
152 332  FORMAT(' ','cLqq2   = ',D20.10)
153
154      RPSL2R=DAIND(0.0D0,1.0D0,FPSL2R,EPS,KEY,MAX,KOUNT,EST)
155      RPSL2S=0.0D0
156      FLPSR2=RPSL2R
157      FLPSL2=RPSL2S
158      FLPS2=FLPSR2+FLPSL2
159      WRITE(*,334) FLPS2
160 334  FORMAT(' ','cLqqps2 = ',D20.10)
161
162      RGL1R=DAIND(0.0D0,1.0D0,FLG1R,EPS,KEY,MAX,KOUNT,EST)
163      RGL1S=0.0D0
164      FLGR1=RGL1R
165      FLGL1=RGL1S
166      FLG1=FLGR1+FLGL1
167      WRITE(*,335) FLG1
168 335  FORMAT(' ','cLgg1   = ',D20.10)
169
170      RGL2R=DAIND(0.0D0,1.0D0,FLG2R,EPS,KEY,MAX,KOUNT,EST)
171      RGL2S=0.0D0
172      FLGR2=RGL2R
173      FLGL2=RGL2S
174      FLG2=FLGR2+FLGL2
175      WRITE(*,336) FLG2
176 336  FORMAT(' ','cLgg2   = ',D20.10)
177
178csm end added
179
180      WRITE(*,325) DN
181      R30R=DAIND(0.0D0,1.0D0,F30R,EPS,KEY,MAX,KOUNT,EST)
182      WRITE(*,31) R30R
183   31 FORMAT(' ','c3qq0   = ',D20.10)
184      R31R=DAIND(0.0D0,1.0D0,F31R,EPS,KEY,MAX,KOUNT,EST)
185      IF(FA.GT.1.1D0) THEN
186      R31S=DAIND(0.0D0,1.0D0,F31S,EPS,KEY,MAX,KOUNT,EST)
187      ELSE
188      R31S=0.0D0
189      ENDIF
190      F3QR1=R31R
191      F3QL1=R31S
192      F3QVS1=COR1(0.0D0)
193      F3Q1=F3QR1+F3QL1+F3QVS1
194      WRITE(*,35) F3Q1
195      R32R=DAIND(0.0D0,1.0D0,F32R,EPS,KEY,MAX,KOUNT,EST)
196      IF(FA.GT.1.1D0) THEN
197      R32S=DAIND(0.0D0,1.0D0,F32S,EPS,KEY,MAX,KOUNT,EST)
198      ELSE
199      R32S=0.0D0
200      ENDIF
201      F3QR2=R32R
202      F3QL2=R32S
203      F3QVS2=COR2(0.0D0)
204      F3Q2=F3QR2+F3QL2+F3QVS2
205      WRITE(*,39) F3Q2
206   39 FORMAT(' ','c3qq2   = ',D20.10)
207*
208* csm polarized coefficient functions
209*
210      R1G1=DAIND(0.0D0,1.0D0,F1G1,EPS,KEY,MAX,KOUNT,EST)
211      G1G1=R1G1
212csm      WRITE(*,41) G1G1
213   41 FORMAT(' ','THE ORDER ALPHAS CONTRIBUTION TO G1 =',D20.10)
214      R1G2=DAIND(0.0D0,1.0D0,F1G2,EPS,KEY,MAX,KOUNT,EST)
215      G1G2=R1G2
216csm      WRITE(*,43) G1G2
217   43 FORMAT(' ','ORDER ALPHAS**2 GLUON CONTRIBUTION TO G1 =',D20.10)
218      R1Q2=DAIND(0.0D0,1.0D0,F1Q2,EPS,KEY,MAX,KOUNT,EST)
219      G1Q2=R1Q2
220csm      WRITE(*,44) G1Q2
221   44 FORMAT(' ','ORDER ALPHAS**2 QUARK CONTRIBUTION TO G1 =',D20.10)
222*
223*
224      END
225
226
227C THIS IS THE WILSON COEFFICIENT CORRESPONDING TO THE PROCESS :
228C V + Q -> Q + G  (SINGLET AND NONSINGLET) , F2(X,Q^2) , MSBAR SCHEME .
229
230      DOUBLE PRECISION FUNCTION C2NS1(Y)
231      IMPLICIT REAL*8(A-H,O-Z)
232      DOUBLE PRECISION NF
233      COMMON/PAR/NF,Q2,FS,RS,DN
234      COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
235      FS2=FS*FS
236      DL1=DLOG(Q2/FS2)
237      DL=DLOG(Y)
238      DLM=DLOG(1.0D0-Y)
239      DM=1.0D0/(1.0D0-Y)
240      A=2.0D0*(1.0D0+Y)*(DL-DLM-DL1)-4.0D0*DM*DL+6.0D0+4.0D0*Y
241      C2NS1=CF*A
242      RETURN
243      END
244
245
246C THIS IS THE WILSON COEFFICIENT CORRESPONDING TO THE PROCESS :
247C V + Q -> Q + G  (SINGLET AND NONSINGLET) , F3(X,Q^2) , MSBAR SCHEME .
248
249
250      DOUBLE PRECISION FUNCTION C3NS1(Y)
251      IMPLICIT REAL*8(A-H,O-Z)
252      DOUBLE PRECISION NF
253      COMMON/PAR/NF,Q2,FS,RS,DN
254      COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
255      FS2=FS*FS
256      DL1=DLOG(Q2/FS2)
257      DL=DLOG(Y)
258      DLM=DLOG(1.0D0-Y)
259      DM=1.0D0/(1.0D0-Y)
260      A=2.0D0*(1.0D0+Y)*(DL-DLM-DL1)-4.0D0*DM*DL+4.0D0+2.0D0*Y
261      C3NS1=CF*A
262      RETURN
263      END
264
265C THIS IS THE WILSON COEFFICIENT CORRESPONDING TO THE PROCESS :
266C V + Q -> Q + G  (SINGLET AND NONSINGLET) , F2(X,Q^2) AND F3(X,Q^2) ,
267C MSBAR SCHEME .THIS IS THE SOFT GLUON CONTRIBUTION.
268
269
270
271      DOUBLE PRECISION FUNCTION SOFT1(Y)
272      IMPLICIT REAL*8(A-H,O-Z)
273      DOUBLE PRECISION NF
274      COMMON/PAR/NF,Q2,FS,RS,DN
275      COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
276      FS2=FS*FS
277      DL1=DLOG(Q2/FS2)
278      DLM=DLOG(1.0D0-Y)
279      DM=1.0D0/(1.0D0-Y)
280      A=4.0D0*(DLM+DL1)-3.0D0
281      SOFT1=CF*A*DM
282      RETURN
283      END
284
285C THIS IS THE WILSON COEFFICIENT CORRESPONDING TO THE PROCESS :
286C V + Q -> Q + G  (SINGLET AND NONSINGLET) , F2(X,Q^2)
287C MSBAR SCHEME . THIS IS THE SOFT PLUS VIRTUAL GLUON CONTRIBUTION.
288
289
290
291      DOUBLE PRECISION FUNCTION BACOR1(Z)
292      IMPLICIT REAL*8(A-H,O-Z)
293      DOUBLE PRECISION NF
294      COMMON/PAR/NF,Q2,FS,RS,DN
295      COMMON/CHOICE/MODE
296      COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
297      ZETA2=1.6449340668D0
298      FS2=FS*FS
299      DL1=DLOG(Q2/FS2)
300      CALL STRUCT(Z,FS,MODE,UV,DV,SS,ST,CH,BOT,GL)
301      SIGMA=UV+DV+4.0D0*SS+2.0D0*(ST+CH)
302      FA=SIGMA
303      DLM=DLOG(1.0D0-Z)
304      C1=2.0D0*DLM**2-3.0D0*DLM-9.0D0-4.0D0*ZETA2
305      C2=4.0D0*DLM+3.0D0
306      BACOR1=4.0D0*(C1+C2*DL1)*FA/3.0D0
307      RETURN
308      END
309
310      DOUBLE PRECISION FUNCTION COR1(Z)
311      IMPLICIT REAL*8(A-H,O-Z)
312      DOUBLE PRECISION NF
313      COMMON/PAR/NF,Q2,FS,RS,DN
314      COMMON/CHOICE/MODE
315      COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
316      ZETA2=1.6449340668D0
317      FS2=FS*FS
318      DL1=DLOG(Q2/FS2)
319      DLM=DLOG(1.0D0-Z)
320      C1=2.0D0*DLM**2-3.0D0*DLM-9.0D0-4.0D0*ZETA2
321      C2=4.0D0*DLM+3.0D0
322      COR1=CF*(C1+C2*DL1)
323      RETURN
324      END
325
326C THIS IS THE WILSON COEFFICIENT CORRESPONDING TO THE PROCESS :
327C V + Q -> Q + G  (NONSINGLET) , F3(X,Q^2) , MSBAR SCHEME .
328C THIS IS THE SOFT PLUS VIRTUAL GLUON CONTRIBUTION.
329
330
331
332      DOUBLE PRECISION FUNCTION COR3(Z)
333      IMPLICIT REAL*8(A-H,O-Z)
334      DOUBLE PRECISION NF
335      COMMON/PAR/NF,Q2,FS,RS,DN
336      COMMON/CHOICE/MODE
337      COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
338      ZETA2=1.6449340668D0
339      FS2=FS*FS
340      DL1=DLOG(Q2/FS2)
341      CALL STRUCT(Z,FS,MODE,UV,DV,SS,ST,CH,BOT,GL)
342      FA=UV+DV
343      DLM=DLOG(1.0D0-Z)
344      C1=2.0D0*DLM**2-3.0D0*DLM-9.0D0-4.0D0*ZETA2
345      C2=4.0D0*DLM+3.0D0
346      COR3=CF*(C1+C2*DL1)*FA
347      RETURN
348      END
349
350C THIS IS THE WILSON COEFFICIENT CORRESPONDING TO THE PROCESS :
351C V + Q -> Q + G  (NONSINGLET) , F3(X,Q^2) , DIS SCHEME .
352
353
354      DOUBLE PRECISION FUNCTION D3NS1(Y)
355      IMPLICIT REAL*8(A-H,O-Z)
356      DOUBLE PRECISION NF
357      COMMON/PAR/NF,Q2,FS,RS,DN
358      COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
359      FS2=FS*FS
360      DL1=DLOG(Q2/FS2)
361      A=-2.0D0*(1.0D0+Y)*DL1-2.0D0-2.0D0*Y
362      D3NS1=CF*A
363      RETURN
364      END
365
366
367C THIS IS THE WILSON COEFFICIENT CORRESPONDING TO THE PROCESS :
368C V + Q -> Q + G  (SINGLET AND NONSINGLET) , F2(X,Q^2) AND F3(X,Q^2) ,
369C DIS SCHEME .THIS IS THE SOFT GLUON CONTRIBUTION.
370
371
372
373      DOUBLE PRECISION FUNCTION DSOFT1(Y)
374      IMPLICIT REAL*8(A-H,O-Z)
375      DOUBLE PRECISION NF
376      COMMON/PAR/NF,Q2,FS,RS,DN
377      COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
378      FS2=FS*FS
379      DL1=DLOG(Q2/FS2)
380      DM=1.0D0/(1.0D0-Y)
381      A=4.0D0*DL1
382      DSOFT1=CF*A*DM
383      RETURN
384      END
385
386
387C THIS IS THE WILSON COEFFICIENT CORRESPONDING TO THE PROCESS :
388C V + Q -> Q + G  (SINGLET AND NONSINGLET) , F2(X,Q^2)
389C DIS SCHEME . THIS IS THE SOFT PLUS VIRTUAL GLUON CONTRIBUTION.
390
391
392
393      DOUBLE PRECISION FUNCTION DCOR1(Z)
394      IMPLICIT REAL*8(A-H,O-Z)
395      DOUBLE PRECISION NF
396      COMMON/PAR/NF,Q2,FS,RS,DN
397      COMMON/CHOICE/MODE
398      COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
399      FS2=FS*FS
400      DL1=DLOG(Q2/FS2)
401      CALL STRUCT(Z,FS,MODE,UV,DV,SS,ST,CH,BOT,GL)
402      SIGMA=UV+DV+4.0D0*SS+2.0D0*(ST+CH)
403      FA=SIGMA
404      DLM=DLOG(1.0D0-Z)
405      A=(4.0D0*DLM+3.0D0)*DL1
406      DCOR1=CF*A*FA
407      RETURN
408      END
409
410
411C THIS IS THE WILSON COEFFICIENT CORRESPONDING TO THE PROCESS :
412C V + Q -> Q + G  (NONSINGLET) , F3(X,Q^2) , DIS SCHEME .
413C THIS IS THE SOFT PLUS VIRTUAL GLUON CONTRIBUTION.
414
415
416      DOUBLE PRECISION FUNCTION DCOR3(Z)
417      IMPLICIT REAL*8(A-H,O-Z)
418      DOUBLE PRECISION NF
419      COMMON/PAR/NF,Q2,FS,RS,DN
420      COMMON/CHOICE/MODE
421      COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
422      FS2=FS*FS
423      DL1=DLOG(Q2/FS2)
424      CALL STRUCT(Z,FS,MODE,UV,DV,SS,ST,CH,BOT,GL)
425      FA=UV+DV
426      DLM=DLOG(1.0D0-Z)
427      A=(4.0D0*DLM+3.0D0)*DL1
428      DCOR3=CF*A*FA
429      RETURN
430      END
431
432
433C THIS IS THE WILSON COEFFICIENT CORRESPONDING TO THE PROCESS :
434C V + G -> Q + QB  , F2(X,Q^2) , MSBAR SCHEME .
435
436
437      DOUBLE PRECISION FUNCTION C2G1(Y)
438      IMPLICIT REAL*8(A-H,O-Z)
439      DOUBLE PRECISION NF
440      COMMON/PAR/NF,Q2,FS,RS,DN
441      COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
442      FS2=FS*FS
443      DL1=DLOG(Q2/FS2)
444      DL=DLOG(Y)
445      DLM=DLOG(1.0D0-Y)
446      A1=(2.0D0-4.0D0*Y+4.0D0*Y*Y)*(DLM+DL1-DL)+16.0D0*Y*(1.0D0-Y)
447     1-2.0D0
448      C2G1=NF*A1
449      RETURN
450      END
451
452
453C THIS IS THE WILSON COEFFICIENT CORRESPONDING TO THE PROCESSES :
454C V + Q -> Q + G + G ; V + Q -> Q + Q + QB WITH NONIDENTICAL QUARKS
455C IN THE FINAL STATE (NONSINGLET) , F2(X,Q^2) , MSBAR SCHEME .
456
457      DOUBLE PRECISION FUNCTION C2NSP(Y)
458      IMPLICIT REAL*8(A-H,O-Z)
459      COMPLEX*16 WGPLG
460      DOUBLE PRECISION NF
461      COMMON/PAR/NF,Q2,FS,RS,DN
462      COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
463      FS2=FS*FS
464      RS2=RS*RS
465      DL1=DLOG(Q2/FS2)
466      DL2=DL1*DL1
467      EL1=DLOG(FS2/RS2)
468      S121MX=WGPLG(1,2,1.0D0-Y)
469      S211MX=WGPLG(2,1,1.0D0-Y)
470      S21MX=WGPLG(2,1,-Y)
471      S111MX=WGPLG(1,1,1.0D0-Y)
472      S11MX=WGPLG(1,1,-Y)
473      ZETA2=1.6449340668D0
474      ZETA3=1.2020569031D0
475      DM=1.0D0/(1.0D0-Y)
476      DL=DLOG(Y)
477      DLP=DLOG(1.0D0+Y)
478      DLM=DLOG(1.0D0-Y)
479      A1=DL**3*(3.0D0+3.0D0*Y-8.0D0*DM/3.0D0)+
480     1DL**2*DLM*(-16.0D0-16.0D0*Y+24.0D0*DM)+
481     2DL**2*(16.0D0+36.0D0*Y*Y*Y/5.0D0+24*Y*Y+
482     314.0D0*Y-3.0D0*DM)+
483     4DL*DLM**2*(16.0D0+16.0D0*Y-28.0D0*DM)
484      A2=DL*DLM*(-40.0D0-56.0D0*Y+24.0D0*DM)+
485     2DL*DLP*(40.0D0-72.0D0*Y*Y*Y/5.0D0-48.0D0*Y*Y+8.0D0*Y)+
486     38.0D0*DL*DLP/Y/Y/5.0D0-8.0D0*DL/Y/5.0D0+
487     4ZETA2*DL*(-28.0D0-28.0D0*Y+48.0D0*DM)+
488     5DL*S111MX*(-4.0D0-4.0D0*Y)+
489     6DL*S11MX*(16.0D0+64.0D0*Y-48.0D0*DM)
490      A3=DL*(-146.0D0/5.0D0+72.0D0*Y*Y/5.0D0-356.0D0*Y/5.0D0+
491     161.0D0*DM)+
492     2DLM**3*(-4.0D0-4.0D0*Y)+
493     3DLM**2*(14.0D0+18.0D0*Y)+
494     4ZETA2*DLM*(24.0D0-24.0D0*Y)+
495     5DLM*S111MX*(8.0D0+8.0D0*Y-8.0D0*DM)+
496     6DLM*(-32.0D0+84.0D0*Y)+
497     7ZETA2*(-28.0D0-72.0D0*Y*Y*Y/5.0D0-48.0D0*Y*Y-12.0D0*Y)
498      A4=ZETA3*(-32.0D0-32.0D0*Y)+
499     1S111MX*(-8.0D0-24.0D0*Y-12.0D0*DM)+
500     2S211MX*(-16.0D0-16.0D0*Y+24.0D0*DM)+
501     3S121MX*(4.0D0+52.0D0*Y-24.0D0*DM)+
502     4S11MX*(40.0D0-72.0D0*Y*Y*Y/5.0D0-48.0D0*Y*Y+8.0D0*Y)+
503     58.0D0*S11MX/Y/Y/5.0D0+8.0D0/Y/5.0D0+
504     6S21MX*(-32.0D0-128.0D0*Y+96.0D0*DM)+72.0D0*ZETA3*DM+
505     772.0D0*Y*Y/5.0D0+38.0D0/5.0D0-543.0D0*Y/5.0D0
506      A5=DL*(6.0D0+6.0D0*Y-8.0D0*DM)+
507     1DLM*(-8.0D0-8.0D0*Y)-
508     210.0D0-2.0D0*Y
509      A6=DL**2*(-8.0D0-8.0D0*Y+8.0D0*DM)+
510     1DL*DLM*(28.0D0+28.0D0*Y-48.0D0*DM)+
511     2DL*(-8.0D0-24.0D0*Y-12.0D0*DM)+
512     3DLM**2*(-12.0D0-12.0D0*Y)+
513     4DLM*(16.0D0+24.0D0*Y)+
514     5ZETA2*(16.0D0+16.0D0*Y)+
515     6S111MX*(4.0D0+4.0D0*Y)+
516     722.0D0+44.0D0*Y
517      FCF=CF2*(A1+A2+A3+A4+A5*DL2+A6*DL1)
518      B1=DL**3*(1.0D0+Y-2.0D0*DM)+
519     1DL**2*DLM*(-2.0D0-2.0D0*Y+4.0D0*DM)+
520     2DL**2*(43.0D0/6.0D0-18.0D0*Y*Y*Y/5.0D0-12.0D0*Y*Y+
521     367.0D0*Y/6.0D0-55.0D0*DM/3.0D0)+
522     4DL*DLM*(-32.0D0/3.0D0-32.0D0*Y/3.0D0+88.0D0*DM/3.0D0)+
523     5DL*DLP*(-20.0D0+36.0D0*Y*Y*Y/5.0D0+24.0D0*Y*Y-4.0D0*Y)
524     6-4.0D0*DL*DLP/Y/Y/5.0D0+4.0D0*DL/Y/5.0D0
525      B2=DL*S111MX*(-4.0D0-4.0D0*Y+8.0D0*DM)+
526     1DL*S11MX*(-8.0D0-32.0D0*Y+24.0D0*DM)+
527     2DL*(604.0D0/15.0D0-36.0D0*Y*Y/5.0D0+1474.0D0*Y/15.0D0-
528     3239.0D0*DM/3.0D0)+
529     4DLM**2*(11.0D0/3.0D0+11.0D0*Y/3.0D0)+
530     5ZETA2*DLM*(24.0D0*Y)+
531     6DLM*S111MX*(-4.0D0-4.0D0*Y+8.0D0*DM)
532      B3=DLM*(16.0D0/9.0D0-740.0D0*Y/9.0D0)+
533     1ZETA2*(-28.0D0/3.0D0+36.0D0*Y*Y*Y/5.0D0+24.0D0*Y*Y
534     2-52.0D0*Y/3.0D0)+ZETA3*(-2.0D0-2.0D0*Y)+
535     3S111MX*(-10.0D0/3.0D0-10.0D0*Y/3.0D0+44.0D0*DM/3.0D0)+
536     4S211MX*(12.0D0+12.0D0*Y-24.0D0*DM)+
537     5S121MX*(-8.0D0-32.0D0*Y+24.0D0*DM)
538      B4=S11MX*(-20.0D0+36.0D0*Y*Y*Y/5.0D0+24.0D0*Y*Y-4.0D0*Y)
539     1-4.0D0*S11MX/Y/Y/5.0D0-4.0D0/Y/5.0D0+
540     2S21MX*(16.0D0+64.0D0*Y-48.0D0*DM)-36.0D0*ZETA3*DM+
541     31522.0D0/135.0D0-36.0D0*Y*Y/5.0D0+18733.0D0*Y/135.0D0
542      B5=11.0D0/3.0D0+11.0D0*Y/3.0D0
543      B6=DL**2*(-2.0D0-2.0D0*Y+4.0D0*DM)+
544     1DL*(-32.0D0/3.0D0-32.0D0*Y/3.0D0+88.0D0*DM/3.0D0)+
545     2DLM*(22.0D0/3.0D0+22.0D0*Y/3.0D0)+
546     3ZETA2*(4.0D0+4.0D0*Y)-
547     492.0D0/9.0D0-506.0D0*Y/9.0D0
548      B7=(22.0D0/3.0D0+22.0D0*Y/3.0D0)*(DLM-DL)+
549     144.0D0*DL*DM/3.0D0-22.0D0-44.0D0*Y/3.0D0
550      FCA=CACF*(B1+B2+B3+B4+B5*DL2+2.0D0*B5*DL1*EL1+B6*DL1+
551     1B7*EL1)
552      C1=DL**2*(-5.0D0/3.0D0-5.0D0*Y/3.0D0+10.0D0*DM/3.0D0)+
553     1DL*DLM*(8.0D0/3.0D0+8.0D0*Y/3.0D0-16.0D0*DM/3.0D0)+
554     2DL*(-26.0D0/3.0D0-38.0D0*Y/3.0D0+38.0D0*DM/3.0D0)+
555     3DLM**2*(-2.0D0/3.0D0-2.0D0*Y/3.0D0)+
556     4DLM*(32.0D0/9.0D0+68.0D0*Y/9.0D0)+
557     5ZETA2*(4.0D0/3.0D0+4.0D0*Y/3.0D0)+
558     6S111MX*(4.0D0/3.0D0+4.0D0*Y/3.0D0-8.0D0*DM/3.0D0)
559     7-158.0D0/27.0D0-488.0D0*Y/27.0D0
560      C2=-2.0D0/3.0D0-2.0D0*Y/3.0D0
561      C3=DL*(8.0D0/3.0D0+8.0D0*Y/3.0D0-16.0D0*DM/3.0D0)+
562     1DLM*(-4.0D0/3.0D0-4.0D0*Y/3.0D0)+
563     232.0D0/9.0D0+68.0D0*Y/9.0D0
564      C4=(-4.0D0/3.0D0-4.0D0*Y/3.0D0)*(DLM-DL)-
565     18.0D0*DL*DM/3.0D0+4.0D0+8.0D0*Y/3.0D0
566      FNF=2.0D0*NF*CFTF*(C1+C2*DL2+2.0D0*C2*DL1*EL1+C3*DL1+C4*EL1)
567      C2NSP=FCF+FCA+FNF
568      RETURN
569      END
570
571C THIS IS THE WILSON COEFFICIENT CORRESPONDING TO THE PROCESSES :
572C V + Q -> Q + G + G ; V + Q -> Q + Q + QB WITH NONIDENTICAL QUARKS
573C IN THE FINAL STATE (NONSINGLET) , F3(X,Q^2) , MS BAR SCHEME .
574
575      DOUBLE PRECISION FUNCTION C3NSP(Y)
576      IMPLICIT REAL*8(A-H,O-Z)
577      COMPLEX*16 WGPLG
578      DOUBLE PRECISION NF
579      COMMON/PAR/NF,Q2,FS,RS,DN
580      COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
581      FS2=FS*FS
582      RS2=RS*RS
583      DL1=DLOG(Q2/FS2)
584      DL2=DL1*DL1
585      EL1=DLOG(FS2/RS2)
586      S121MX=WGPLG(1,2,1.0D0-Y)
587      S211MX=WGPLG(2,1,1.0D0-Y)
588      S21MX=WGPLG(2,1,-Y)
589      S111MX=WGPLG(1,1,1.0D0-Y)
590      S11MX=WGPLG(1,1,-Y)
591      ZETA2=1.6449340668D0
592      ZETA3=1.2020569031D0
593      DM=1.0D0/(1.0D0-Y)
594      DL=DLOG(Y)
595      DLP=DLOG(1.0D0+Y)
596      DLM=DLOG(1.0D0-Y)
597      A1=DL**3*(3.0D0+3.0D0*Y-8.0D0*DM/3.0D0)+
598     1DL**2*DLM*(-16.0D0-16.0D0*Y+24.0D0*DM)+
599     2DL**2*(16.0D0+4*Y*Y+22.0D0*Y-3.0D0*DM)+
600     3DL*DLM**2*(16.0D0+16.0D0*Y-28.0D0*DM)
601      A2=DL*DLM*(-28.0D0-44.0D0*Y+24.0D0*DM)+
602     1DL*DLP*(-8.0D0-8.0D0*Y*Y-8.0D0*Y-8.0D0/Y)+
603     2ZETA2*DL*(-28.0D0-28.0D0*Y+48.0D0*DM)+
604     3DL*S111MX*(-4.0D0-4.0D0*Y)+
605     4DL*S11MX*(32.0D0+16.0D0*Y-48.0D0*DM)
606      A3=DL*(-26.0D0-32.0D0*Y+61.0D0*DM)+
607     1DLM**3*(-4.0D0-4.0D0*Y)+
608     2DLM**2*(10.0D0+14.0D0*Y)+
609     3ZETA2*DLM*(8.0D0+24.0D0*Y)+
610     4DLM*S111MX*(8.0D0+8.0D0*Y-8.0D0*DM)+
611     5DLM*(10.0D0+6.0D0*Y)+
612     6ZETA2*(-28.0D0-8.0D0*Y*Y-28.0D0*Y)
613      A4=ZETA3*(-32.0D0-32.0D0*Y)+
614     1S111MX*(-4.0D0-20.0D0*Y-12.0D0*DM)+
615     2S211MX*(-16.0D0-16.0D0*Y+24.0D0*DM)+
616     3S121MX*(20.0D0+4.0D0*Y-24.0D0*DM)+
617     4S11MX*(-8.0D0-8.0D0*Y*Y-8.0D0*Y-8.0D0/Y)+
618     6S21MX*(-64.0D0-32.0D0*Y+96.0D0*DM)+72.0D0*ZETA3*DM-
619     741.0D0-10.0D0*Y
620      A5=DL*(6.0D0+6.0D0*Y-8.0D0*DM)+
621     1DLM*(-8.0D0-8.0D0*Y)-
622     210.0D0-2.0D0*Y
623      A6=DL**2*(-8.0D0-8.0D0*Y+8.0D0*DM)+
624     1DL*DLM*(28.0D0+28.0D0*Y-48.0D0*DM)+
625     2DL*(-4.0D0-20.0D0*Y-12.0D0*DM)+
626     3DLM**2*(-12.0D0-12.0D0*Y)+
627     4DLM*(8.0D0+16.0D0*Y)+
628     5ZETA2*(16.0D0+16.0D0*Y)+
629     6S111MX*(4.0D0+4.0D0*Y)+
630     716.0D0+38.0D0*Y
631      FCF=CF2*(A1+A2+A3+A4+A5*DL2+A6*DL1)
632      B1=DL**3*(1.0D0+Y-2.0D0*DM)+
633     1DL**2*DLM*(-2.0D0-2.0D0*Y+4.0D0*DM)+
634     2DL**2*(31.0D0/6.0D0-2.0D0*Y*Y+31.0D0*Y/6.0D0-
635     355.0D0*DM/3.0D0)+
636     4DL*DLM*(-32.0D0/3.0D0-32.0D0*Y/3.0D0+88.0D0*DM/3.0D0)+
637     5DL*DLP*(4.0D0+4.0D0*Y*Y+4.0D0*Y+4.0D0/Y)
638      B2=DL*S111MX*(-4.0D0-4.0D0*Y+8.0D0*DM)+
639     1DL*S11MX*(-16.0D0-8.0D0*Y+24.0D0*DM)+
640     2DL*(32.0D0+66.0D0*Y-239.0D0*DM/3.0D0)+
641     3DLM**2*(11.0D0/3.0D0+11.0D0*Y/3.0D0)+
642     4ZETA2*DLM*(8.0D0)+
643     5DLM*S111MX*(-4.0D0-4.0D0*Y+8.0D0*DM)
644      B3=DLM*(-134.0D0/9.0D0-314.0D0*Y/9.0D0)+
645     1ZETA2*(-16.0D0/3.0D0+4.0D0*Y*Y-16.0D0*Y/3.0D0)+
646     2ZETA3*(-2.0D0-2.0D0*Y)+
647     3S111MX*(-10.0D0/3.0D0-10.0D0*Y/3.0D0+44.0D0*DM/3.0D0)+
648     4S211MX*(12.0D0+12.0D0*Y-24.0D0*DM)+
649     5S121MX*(-16.0D0-8.0D0*Y+24.0D0*DM)
650      B4=S11MX*(4.0D0+4.0D0*Y*Y+4.0D0*Y+4.0D0/Y)+
651     1S21MX*(32.0D0+16.0D0*Y-48.0D0*DM)-36.0D0*ZETA3*DM+
652     2725.0D0/27.0D0+1820.0D0*Y/27.0D0
653      B5=11.0D0/3.0D0+11.0D0*Y/3.0D0
654      B6=DL**2*(-2.0D0-2.0D0*Y+4.0D0*DM)+
655     1DL*(-32.0D0/3.0D0-32.0D0*Y/3.0D0+88.0D0*DM/3.0D0)+
656     2DLM*(22.0D0/3.0D0+22.0D0*Y/3.0D0)+
657     3ZETA2*(4.0D0+4.0D0*Y)-
658     426.0D0/9.0D0-440.0D0*Y/9.0D0
659      B7=(22.0D0/3.0D0+22.0D0*Y/3.0D0)*(DLM-DL)+
660     144.0D0*DL*DM/3.0D0-44.0D0/3.0D0-22.0D0*Y/3.0D0
661      FCA=CACF*(B1+B2+B3+B4+B5*DL2+2.0D0*B5*DL1*EL1+B6*DL1+
662     1B7*EL1)
663      C1=DL**2*(-5.0D0/3.0D0-5.0D0*Y/3.0D0+10.0D0*DM/3.0D0)+
664     1DL*DLM*(8.0D0/3.0D0+8.0D0*Y/3.0D0-16.0D0*DM/3.0D0)+
665     2DL*(-6.0D0-10.0D0*Y+38.0D0*DM/3.0D0)+
666     3DLM**2*(-2.0D0/3.0D0-2.0D0*Y/3.0D0)+
667     4DLM*(20.0D0/9.0D0+56.0D0*Y/9.0D0)+
668     5ZETA2*(4.0D0/3.0D0+4.0D0*Y/3.0D0)+
669     6S111MX*(4.0D0/3.0D0+4.0D0*Y/3.0D0-8.0D0*DM/3.0D0)
670     7-116.0D0/27.0D0-302.0D0*Y/27.0D0
671      C2=-2.0D0/3.0D0-2.0D0*Y/3.0D0
672      C3=DL*(8.0D0/3.0D0+8.0D0*Y/3.0D0-16.0D0*DM/3.0D0)+
673     1DLM*(-4.0D0/3.0D0-4.0D0*Y/3.0D0)+
674     220.0D0/9.0D0+56.0D0*Y/9.0D0
675      C4=(-4.0D0/3.0D0-4.0D0*Y/3.0D0)*(DLM-DL)-
676     18.0D0*DL*DM/3.0D0+8.0D0/3.0D0+4.0D0*Y/3.0D0
677      FNF=2.0D0*NF*CFTF*(C1+C2*DL2+2.0D0*C2*DL1*EL1+C3*DL1+C4*EL1)
678      C3NSP=FCF+FCA+FNF
679C      C3NSP=FNF
680      RETURN
681      END
682
683C THIS IS THE WILSON COEFFICIENT CORRESPONDING TO THE PROCESSES :
684C V + Q -> Q + G + G ; V + Q -> Q + Q + QB WITH NONIDENTICAL QUARKS
685C IN THE FINAL STATE (NONSINGLET) ,F2(X,Q^2) AND F3(X,Q^2) ,
686C MSBAR SCHEME .THIS IS THE SOFT GLUON CONTRIBUTION.
687
688      DOUBLE PRECISION FUNCTION SOFT2(Y)
689      IMPLICIT REAL*8(A-H,O-Z)
690      DOUBLE PRECISION NF
691      COMMON/PAR/NF,Q2,FS,RS,DN
692      COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
693      FS2=FS*FS
694      RS2=RS*RS
695      DL1=DLOG(Q2/FS2)
696      DL2=DL1*DL1
697      EL1=DLOG(FS2/RS2)
698      ZETA2=1.6449340668D0
699      ZETA3=1.2020569031D0
700      DLM=DLOG(1.0D0-Y)
701      DM=1.0D0/(1.0D0-Y)
702      FCF=8.0D0*DLM**3-18.0D0*DLM**2-DLM*(32.0D0*ZETA2
703     1+27.0D0)-8.0D0*ZETA3+36.0D0*ZETA2+25.5D0+
704     2(16.0D0*DLM+12.0D0)*DL2+(24.0D0*DLM**2-12.0D0*DLM-32.0D0*
705     3ZETA2-45.0D0)*DL1
706      FCA=-22.0D0*DLM**2/3.0D0+DLM*(-8.0D0*ZETA2+367.0D0/9.0D0)+
707     140.0D0*ZETA3+44.0D0*ZETA2/3.0D0-3155.0D0/54.0D0-
708     222.0D0*(DL2+2.0D0*DL1*EL1)/3.0D0+(-44.0D0*DLM/3.0D0-8.0D0*
709     3ZETA2+367.0D0/9.0D0)*DL1+(-44.0D0*DLM/3.0D0+11.0D0)*EL1
710      FNF=4.0D0*DLM**2/3.0D0-58.0D0*DLM/9.0D0-8.0D0*ZETA2/3.0D0
711     1+247.0D0/27.0D0+4.0D0*(DL2+2.0D0*DL1*EL1)/3.0D0+
712     2(8.0D0*DLM/3.0D0-58.0D0/9.0D0)*DL1+(8.0D0*DLM/3.0D0-2.0D0)
713     3*EL1
714      SOFT2=DM*(CF2*FCF+CACF*FCA+2.0D0*NF*CFTF*FNF)
715      RETURN
716      END
717
718C THIS IS THE WILSON COEFFICIENT CORRESPONDING TO THE PROCESSES :
719C V + Q -> Q + G + G ; V + Q -> Q + Q + QB WITH NONIDENTICAL QUARKS
720C IN THE FINAL STATE (NONSINGLET) ,F2(X,Q^2) , MS BAR SCHEME .
721C THIS IS THE SOFT PLUS VIRTUAL GLUON CONTRIBUTION.
722
723
724      DOUBLE PRECISION FUNCTION BACOR2(Z)
725      IMPLICIT REAL*8(A-H,O-Z)
726      DOUBLE PRECISION NF
727      COMMON/PAR/NF,Q2,FS,RS,DN
728      COMMON/CHOICE/MODE
729      COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
730      ZETA2=1.6449340668D0
731      ZETA3=1.2020569031D0
732      FS2=FS*FS
733      RS2=RS*RS
734      DL1=DLOG(Q2/FS2)
735      DL2=DL1*DL1
736      EL1=DLOG(FS2/RS2)
737      CALL STRUCT(Z,FS,MODE,UV,DV,SS,ST,CH,BOT,GL)
738      SIGMA=UV+DV+4.0D0*SS+2.0D0*(ST+CH)
739      FA=SIGMA
740      DLM=DLOG(1.0D0-Z)
741      A1=2.0D0*DLM**4-6.0D0*DLM**3-(16.0D0*ZETA2+13.5D0)*DLM**2+
742     1(-8.0D0*ZETA3+36.0D0*ZETA2+25.5D0)*DLM+6.0D0*ZETA2**2-
743     278.0D0*ZETA3+69.0D0*ZETA2+331/8.0D0
744      A2=8.0D0*DLM**2+12.0D0*DLM+4.5D0-8.0D0*ZETA2
745      A3=8.0D0*DLM**3-6.0D0*DLM**2+(45.0D0-32.0D0*ZETA2)*DLM-25.5D0
746     1-12.0D0*ZETA2+40.0D0*ZETA3
747      FCF=CF2*(A1+A2*DL2+A3*DL1)
748      B1=-22.0D0*DLM**3/9.0D0+(-4.0D0*ZETA2+367.0D0/18.0D0)*
749     1DLM**2+(40.0D0*ZETA3+44.0D0*ZETA2/3.0D0-3155.0D0/54.0D0)*DLM+
750     271.0D0*ZETA2**2/5.0D0+140.0D0*ZETA3/3.0D0-251.0D0*ZETA2/3.0D0-
751     35464.0D0/72.0D0
752      B2=-22.0D0*DLM/3.0D0-5.5D0
753      B3=-22.0D0*DLM**2/3.0D0+(367.0D0/9.0D0-8.0D0*ZETA2)*DLM+
754     1215.0D0/6.0D0+88.0D0*ZETA2/3.0D0-12.0D0*ZETA3
755      B4=-22.0D0*DLM**2/3.0D0+11.0D0*DLM+33.0D0+44.0D0*ZETA2/3.0D0
756      FCA=CACF*(B1+B2*DL2+2.0D0*B2*DL1*EL1+B3*DL1+B4*EL1)
757      C1=4.0D0*DLM**3/9.0D0-29.0D0*DLM**2/9.0D0+(-8.0D0*ZETA2/3.0D0+
758     1247.0D0/27.0D0)*DLM+4.0D0*ZETA3/3.0D0+38.0D0*ZETA2/3.0D0+
759     2457.0D0/36.0D0
760      C2=4.0D0*DLM/3.0D0+1.0D0
761      C3=4.0D0*DLM**2/3.0D0-58.0D0*DLM/9.0D0-19.0D0/3.0D0-
762     116.0D0*ZETA2/3.0D0
763      C4=4.0D0*DLM**2/3.0D0-2.0D0*DLM-6.0D0-8.0D0*ZETA2/3.0D0
764      FNF=2.0D0*NF*CFTF*(C1+C2*DL2+2.0D0*C2*DL1*EL1+C3*DL1+C4*EL1)
765      BACOR2=(FCF+FCA+FNF)*FA
766      RETURN
767      END
768
769      DOUBLE PRECISION FUNCTION COR2(Z)
770      IMPLICIT REAL*8(A-H,O-Z)
771      DOUBLE PRECISION NF
772      COMMON/PAR/NF,Q2,FS,RS,DN
773      COMMON/CHOICE/MODE
774      COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
775      ZETA2=1.6449340668D0
776      ZETA3=1.2020569031D0
777      FS2=FS*FS
778      RS2=RS*RS
779      DL1=DLOG(Q2/FS2)
780      DL2=DL1*DL1
781      EL1=DLOG(FS2/RS2)
782      DLM=DLOG(1.0D0-Z)
783      A1=2.0D0*DLM**4-6.0D0*DLM**3-(16.0D0*ZETA2+13.5D0)*DLM**2+
784     1(-8.0D0*ZETA3+36.0D0*ZETA2+25.5D0)*DLM+6.0D0*ZETA2**2-
785     278.0D0*ZETA3+69.0D0*ZETA2+331/8.0D0
786      A2=8.0D0*DLM**2+12.0D0*DLM+4.5D0-8.0D0*ZETA2
787      A3=8.0D0*DLM**3-6.0D0*DLM**2+(45.0D0-32.0D0*ZETA2)*DLM-25.5D0
788     1-12.0D0*ZETA2+40.0D0*ZETA3
789      FCF=CF2*(A1+A2*DL2+A3*DL1)
790      B1=-22.0D0*DLM**3/9.0D0+(-4.0D0*ZETA2+367.0D0/18.0D0)*
791     1DLM**2+(40.0D0*ZETA3+44.0D0*ZETA2/3.0D0-3155.0D0/54.0D0)*DLM+
792     271.0D0*ZETA2**2/5.0D0+140.0D0*ZETA3/3.0D0-251.0D0*ZETA2/3.0D0-
793     35465.0D0/72.0D0
794      B2=-22.0D0*DLM/3.0D0-5.5D0
795      B3=-22.0D0*DLM**2/3.0D0+(367.0D0/9.0D0-8.0D0*ZETA2)*DLM+
796     1215.0D0/6.0D0+88.0D0*ZETA2/3.0D0-12.0D0*ZETA3
797      B4=-22.0D0*DLM**2/3.0D0+11.0D0*DLM+33.0D0+44.0D0*ZETA2/3.0D0
798      FCA=CACF*(B1+B2*DL2+2.0D0*B2*DL1*EL1+B3*DL1+B4*EL1)
799      C1=4.0D0*DLM**3/9.0D0-29.0D0*DLM**2/9.0D0+(-8.0D0*ZETA2/3.0D0+
800     1247.0D0/27.0D0)*DLM+4.0D0*ZETA3/3.0D0+38.0D0*ZETA2/3.0D0+
801     2457.0D0/36.0D0
802      C2=4.0D0*DLM/3.0D0+1.0D0
803      C3=4.0D0*DLM**2/3.0D0-58.0D0*DLM/9.0D0-19.0D0/3.0D0-
804     116.0D0*ZETA2/3.0D0
805      C4=4.0D0*DLM**2/3.0D0-2.0D0*DLM-6.0D0-8.0D0*ZETA2/3.0D0
806      FNF=2.0D0*NF*CFTF*(C1+C2*DL2+2.0D0*C2*DL1*EL1+C3*DL1+C4*EL1)
807      COR2=FCF+FCA+FNF
808C      COR2=FNF
809      RETURN
810      END
811
812
813C THIS IS THE WILSON COEFFICIENT CORRESPONDING TO THE PROCESSES :
814C V + Q -> Q + G + G ; V + Q -> Q + Q + QB WITH NONIDENTICAL QUARKS
815C IN THE FINAL STATE (NONSINGLET) ,F3(X,Q^2) , MS BAR SCHEME .
816C THIS IS THE SOFT PLUS VIRTUAL GLUON CONTRIBUTION.
817
818
819      DOUBLE PRECISION FUNCTION COR4(Z)
820      IMPLICIT REAL*8(A-H,O-Z)
821      DOUBLE PRECISION NF
822      COMMON/PAR/NF,Q2,FS,RS,DN
823      COMMON/CHOICE/MODE
824      COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
825      ZETA2=1.6449340668D0
826      ZETA3=1.2020569031D0
827      FS2=FS*FS
828      RS2=RS*RS
829      DL1=DLOG(Q2/FS2)
830      DL2=DL1*DL1
831      EL1=DLOG(FS2/RS2)
832      DLM=DLOG(1.0D0-Z)
833      A1=2.0D0*DLM**4-6.0D0*DLM**3-(16.0D0*ZETA2+13.5D0)*DLM**2+
834     1(-8.0D0*ZETA3+36.0D0*ZETA2+25.5D0)*DLM+6.0D0*ZETA2**2-
835     278.0D0*ZETA3+69.0D0*ZETA2+331/8.0D0
836      A2=8.0D0*DLM**2+12.0D0*DLM+4.5D0-8.0D0*ZETA2
837      A3=8.0D0*DLM**3-6.0D0*DLM**2-(45.0D0+32.0D0*ZETA2)*DLM-25.5D0
838     1-12.0D0*ZETA2+40.0D0*ZETA3
839      FCF=CF2*(A1+A2*DL2+A3*DL1)
840      B1=-22.0D0*DLM**3/9.0D0+(-4.0D0*ZETA2+367.0D0/18.0D0)*
841     1DLM**2+(40.0D0*ZETA3+44.0D0*ZETA2/3.0D0-3155.0D0/54.0D0)*DLM+
842     271.0D0*ZETA2**2/5.0D0+140.0D0*ZETA3/3.0D0-251.0D0*ZETA2/3.0D0-
843     35465.0D0/72.0D0
844      B2=-22.0D0*DLM/3.0D0-5.5D0
845      B3=-22.0D0*DLM**2/3.0D0+(367.0D0/9.0D0-8.0D0*ZETA2)*DLM+
846     1215.0D0/6.0D0+88.0D0*ZETA2/3.0D0-12.0D0*ZETA3
847      B4=-22.0D0*DLM**2/3.0D0+11.0D0*DLM+33.0D0+44.0D0*ZETA2/3.0D0
848      FCA=CACF*(B1+B2*DL2+2.0D0*B2*DL1*EL1+B3*DL1+B4*EL1)
849      C1=4.0D0*DLM**3/9.0D0-29.0D0*DLM**2/9.0D0+(-8.0D0*ZETA2/3.0D0+
850     1247.0D0/27.0D0)*DLM+4.0D0*ZETA3/3.0D0+38.0D0*ZETA2/3.0D0+
851     2457.0D0/36.0D0
852      C2=4.0D0*DLM/3.0D0+1.0D0
853      C3=4.0D0*DLM**2/3.0D0-58.0D0*DLM/9.0D0-19.0D0/3.0D0-
854     116.0D0*ZETA2/3.0D0
855      C4=4.0D0*DLM**2/3.0D0-2.0D0*DLM-6.0D0-8.0D0*ZETA2/3.0D0
856      FNF=2.0D0*NF*CFTF*(C1+C2*DL2+2.0D0*C2*DL1*EL1+C3*DL1+C4*EL1)
857      COR4=(FCF+FCA+FNF)
858      RETURN
859      END
860
861C THIS IS THE WILSON COEFFICIENT CORRESPONDING TO THE PROCESSES :
862C V + Q -> Q + G + G ; V + Q -> Q + Q + QB WITH NONIDENTICAL QUARKS
863C IN THE FINAL STATE (NONSINGLET) , F3(X,Q^2) , DIS SCHEME .
864
865      DOUBLE PRECISION FUNCTION D3NSP(Y)
866      IMPLICIT REAL*8(A-H,O-Z)
867      COMPLEX*16 WGPLG
868      DOUBLE PRECISION NF
869      COMMON/PAR/NF,Q2,FS,RS,DN
870      COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
871      FS2=FS*FS
872      RS2=RS*RS
873      DL1=DLOG(Q2/FS2)
874      DL2=DL1*DL1
875      EL1=DLOG(FS2/RS2)
876      S121MX=WGPLG(1,2,1.0D0-Y)
877      S21MX=WGPLG(2,1,-Y)
878      S111MX=WGPLG(1,1,1.0D0-Y)
879      S11MX=WGPLG(1,1,-Y)
880      ZETA2=1.6449340668D0
881      DM=1.0D0/(1.0D0-Y)
882      DL=DLOG(Y)
883      DLP=DLOG(1.0D0+Y)
884      DLM=DLOG(1.0D0-Y)
885      A1=DL**2*(2-36.0D0*Y*Y*Y/5.0D0-20*Y*Y+10.0D0*Y)+
886     1DL*DLM*(8.0D0+8.0D0*Y)+
887     2DL*DLP*(-48.0D0+72.0D0*Y*Y*Y/5.0D0+40.0D0*Y*Y-16.0D0*Y-
888     38.0D0/Y)-8.0D0*DL*DLP/Y/Y/5.0D0+8.0D0*DL/Y/5.0D0+
889     4DL*S11MX*(16.0D0-48.0D0*Y)+
890     5DL*(-34.0D0/5.0D0-72.0D0*Y*Y/5.0D0+186.0D0*Y/5.0D0)+
891     6ZETA2*DLM*(-16.0D0+48.0D0*Y)
892      A2=DLM*(36.0D0-84.0D0*Y)+
893     1ZETA2*(-8.0D0+72.0D0*Y*Y*Y/5.0D0+40.0D0*Y*Y-24.0D0*Y)+
894     2S111MX*(8.0D0+8.0D0*Y)+
895     3S121MX*(16.0D0-48.0D0*Y)+
896     4S11MX*(-48.0D0+72.0D0*Y*Y*Y/5.0D0+40.0D0*Y*Y-16.0D0*Y-
897     58.0D0/Y)-8.0D0*S11MX/Y/Y/5.0D0-8.0D0/Y/5.0D0+
898     6S21MX*(-32.0D0+96.0D0*Y)-
899     7283.0D0/5.0D0-72.0D0*Y*Y/5.0D0+353.0D0*Y/5.0D0
900      A3=DL*(6.0D0+6.0D0*Y-8.0D0*DM)+
901     1DLM*(-8.0D0-8.0D0*Y)-
902     210.0D0-2.0D0*Y
903      A4=DL**2*(-2.0D0-2.0D0*Y)+
904     1DL*DLM*(8.0D0+8.0D0*Y-16.0D0*DM)+
905     2DL*(4.0D0-4.0D0*Y-12.0D0*DM)+
906     3DLM*(-8.0D0-8.0D0*Y)-
907     426.0D0+14.0D0*Y
908      FCF=CF2*(A1+A2+A3*DL2+A4*DL1)
909      B1=DL**2*(-2.0D0+18.0D0*Y*Y*Y/5.0D0+10.0D0*Y*Y-6.0D0*Y)+
910     1DL*DLP*(24.0D0-36.0D0*Y*Y*Y/5.0D0-20.0D0*Y*Y+8.0D0*Y+
911     24.0D0/Y)+4.0D0*DL*DLP/Y/Y/5.0D0-4.0D0*DL/Y/5.0D0+
912     3DL*S11MX*(-8.0D0+24.0D0*Y)+
913     4DL*(-124.0D0/15.0D0+36.0D0*Y*Y/5.0D0-484.0D0*Y/15.0D0)+
914     5ZETA2*DLM*(8.0D0-24.0D0*Y)
915      B2=DLM*(-50.0D0/3.0D0+142.0D0*Y/3.0D0)+
916     1ZETA2*(4.0D0-36.0D0*Y*Y*Y/5.0D0-20.0D0*Y*Y+12.0D0*Y)+
917     2S121MX*(-8.0D0+24.0D0*Y)+
918     3S11MX*(24.0D0-36.0D0*Y*Y*Y/5.0D0-20.0D0*Y*Y+8.0D0*Y+
919     44.0D0/Y)+4.0D0*S11MX/Y/Y/5.0D0+4.0D0/Y/5.0D0+
920     5S21MX*(16.0D0-48.0D0*Y)+
921     6701.0D0/45.0D0+36.0D0*Y*Y/5.0D0-3211.0D0*Y/45.0D0
922      B3=11.0D0/3.0D0+11.0D0*Y/3.0D0
923      B4=DL**2*(-2.0D0-2.0D0*Y+4.0D0*DM)+
924     1DL*(-32.0D0/3.0D0-32.0D0*Y/3.0D0+88.0D0*DM/3.0D0)+
925     2DLM*(22.0D0/3.0D0+22.0D0*Y/3.0D0)+
926     3ZETA2*(4.0D0+4.0D0*Y)-
927     426.0D0/9.0D0-440.0D0*Y/9.0D0
928      B5=22.0D0/3.0D0+22.0D0*Y/3.0D0
929      FCA=CACF*(B1+B2+B3*DL2+2.0D0*B3*DL1*EL1+B4*DL1+
930     1B5*EL1)
931      C1=DL*(8.0D0/3.0D0+8.0D0*Y/3.0D0)+
932     1DLM*(-4.0D0/3.0D0-4.0D0*Y/3.0D0)+
933     214.0D0/9.0D0+62.0D0*Y/9.0D0
934      C2=-2.0D0/3.0D0-2.0D0*Y/3.0D0
935      C3=DL*(8.0D0/3.0D0+8.0D0*Y/3.0D0-16.0D0*DM/3.0D0)+
936     1DLM*(-4.0D0/3.0D0-4.0D0*Y/3.0D0)+
937     220.0D0/9.0D0+56.0D0*Y/9.0D0
938      C4=-4.0D0/3.0D0-4.0D0*Y/3.0
939      FNF=2.0D0*NF*CFTF*(C1+C2*DL2+2.0D0*C2*DL1*EL1+C3*DL1+C4*EL1)
940      D3NSP=FCF+FCA+FNF
941      RETURN
942      END
943
944C THIS IS THE WILSON COEFFICIENT CORRESPONDING TO THE PROCESSES :
945C V + Q -> Q + G + G ; V + Q -> Q + Q + QB WITH NONIDENTICAL QUARKS
946C IN THE FINAL STATE (NONSINGLET) ,F2(X,Q^2) AND F3(X,Q^2) ,
947C DIS SCHEME .THIS IS THE SOFT GLUON CONTRIBUTION.
948
949      DOUBLE PRECISION FUNCTION DSOFT2(Y)
950      IMPLICIT REAL*8(A-H,O-Z)
951      DOUBLE PRECISION NF
952      COMMON/PAR/NF,Q2,FS,RS,DN
953      COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
954      FS2=FS*FS
955      RS2=RS*RS
956      DL1=DLOG(Q2/FS2)
957      DL2=DL1*DL1
958      EL1=DLOG(FS2/RS2)
959      DLM=DLOG(1.0D0-Y)
960      DM=1.0D0/(1.0D0-Y)
961      FCF=(16.0D0*DLM+12.0D0)*DL2
962      FCA=22.0D0*(DL2+2.0D0*DL1*EL1)/3.0D0
963      FNF=4.0D0*(DL2+2.0D0*DL1*EL1)/3.0D0
964      DSOFT2=DM*(CF2*FCF+CACF*FCA+2.0D0*NF*CFTF*FNF)
965      RETURN
966      END
967
968C THIS IS THE WILSON COEFFICIENT CORRESPONDING TO THE PROCESSES :
969C V + Q -> Q + G + G ; V + Q -> Q + Q + QB WITH NONIDENTICAL QUARKS
970C IN THE FINAL STATE (NONSINGLET) , F2(X,Q^2) , DIS SCHEME .
971C THIS IS THE SOFT PLUS VIRTUAL GLUON CONTRIBUTION.
972
973
974      DOUBLE PRECISION FUNCTION DCOR2(Z)
975      IMPLICIT REAL*8(A-H,O-Z)
976      DOUBLE PRECISION NF
977      COMMON/PAR/NF,Q2,FS,RS,DN
978      COMMON/CHOICE/MODE
979      COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
980      ZETA2=1.6449340668D0
981      FS2=FS*FS
982      RS2=RS*RS
983      DL1=DLOG(Q2/FS2)
984      DL2=DL1*DL1
985      EL1=DLOG(FS2/RS2)
986      CALL STRUCT(Z,FS,MODE,UV,DV,SS,ST,CH,BOT,GL)
987      SIGMA=UV+DV+4.0D0*SS+2.0D0*(ST+CH)
988      FA=SIGMA
989      DLM=DLOG(1.0D0-Z)
990      A1=8.0D0*DLM**2+12.0D0*DLM+4.5D0-8.0D0*ZETA2
991      FCF=CF2*A1*DL2
992      B1=-22.0D0*DLM/3.0D0-5.5D0
993      FCA=CACF*(B1*DL2+2.0D0*B1*DL1*EL1)
994      C1=4.0D0*DLM/3.0D0+1.0D0
995      FNF=2.0D0*NF*CFTF*(C1*DL2+2.0D0*C1*DL1*EL1)
996      DCOR2=(FCF+FCA+FNF)*FA
997      RETURN
998      END
999
1000C THIS IS THE WILSON COEFFICIENT CORRESPONDING TO THE PROCESSES :
1001C V + Q -> Q + G + G ; V + Q -> Q + Q + QB WITH NONIDENTICAL QUARKS
1002C IN THE FINAL STATE (NONSINGLET) , F3(X,Q^2) , DIS SCHEME .
1003C THIS IS THE SOFT PLUS VIRTUAL GLUON CONTRIBUTION.
1004
1005
1006      DOUBLE PRECISION FUNCTION DCOR4(Z)
1007      IMPLICIT REAL*8(A-H,O-Z)
1008      DOUBLE PRECISION NF
1009      COMMON/PAR/NF,Q2,FS,RS,DN
1010      COMMON/CHOICE/MODE
1011      COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
1012      ZETA2=1.6449340668D0
1013      FS2=FS*FS
1014      RS2=RS*RS
1015      DL1=DLOG(Q2/FS2)
1016      DL2=DL1*DL1
1017      EL1=DLOG(FS2/RS2)
1018      CALL STRUCT(Z,FS,MODE,UV,DV,SS,ST,CH,BOT,GL)
1019      FA=UV+DV
1020      DLM=DLOG(1.0D0-Z)
1021      A1=8.0D0*DLM**2+12.0D0*DLM+4.5D0-8.0D0*ZETA2
1022      FCF=CF2*A1*DL2
1023      B1=-22.0D0*DLM/3.0D0-5.5D0
1024      FCA=CACF*(B1*DL2+2.0D0*B1*DL1*EL1)
1025      C1=4.0D0*DLM/3.0D0+1.0D0
1026      FNF=2.0D0*NF*CFTF*(C1*DL2+2.0D0*C1*DL1*EL1)
1027      DCOR4=(FCF+FCA+FNF)*FA
1028      RETURN
1029      END
1030
1031
1032
1033
1034C THIS IS THE WILSON COEFFICIENT CORRESPONDING TO THE PROCESS:
1035C V + Q -> Q + Q + QB WITH IDENTICAL QUARKS IN THE FINAL STATE
1036C (NONSINGLET) , F2(X,Q^2) , MSBAR SCHEME .
1037
1038
1039      DOUBLE PRECISION FUNCTION C2NSM(Y)
1040      IMPLICIT REAL*8(A-H,O-Z)
1041      COMPLEX*16 WGPLG
1042      DOUBLE PRECISION NF
1043      COMMON/PAR/NF,Q2,FS,RS,DN
1044      COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
1045      FS2=FS*FS
1046      DL1=DLOG(Q2/FS2)
1047      S121MX=WGPLG(1,2,1.0D0-Y)
1048      S12MX=WGPLG(1,2,-Y)
1049      S211MX=WGPLG(2,1,1.0D0-Y)
1050      S21MX=WGPLG(2,1,-Y)
1051      Z=(1.0D0-Y)/(1.0D0+Y)
1052      S21Z=WGPLG(2,1,Z)
1053      S21MZ=WGPLG(2,1,-Z)
1054      S111MX=WGPLG(1,1,1.0D0-Y)
1055      S11MX=WGPLG(1,1,-Y)
1056      ZETA2=1.6449340668D0
1057      ZETA3=1.2020569031D0
1058      DP=1.0D0/(1.0D0+Y)
1059      DL=DLOG(Y)
1060      DLP=DLOG(1.0D0+Y)
1061      DLM=DLOG(1.0D0-Y)
1062      B1=DL**3*(2.0D0-2.0D0*Y-4.0D0*DP)+
1063     1DL**2*DLM*(-4.0D0+4.0D0*Y+8.0D0*DP)+
1064     2DL**2*DLP*(-16.0D0+40.0D0*Y+40.0D0*DP)+
1065     3DL**2*(-4.0D0+36.0D0*Y*Y*Y/5.0D0-24.0D0*Y*Y-16.0D0*Y)+
1066     4DL*DLM*DLP*(16.0D0-16.0D0*Y-32.0D0*DP)
1067      B2=DL*DLM*(8.0D0+8.0D0*Y)+
1068     1DL*DLP**2*(-48.0D0*Y-16.0D0*DP)+
1069     2DL*DLP*(32.0D0-72.0D0*Y*Y*Y/5.0D0+48.0D0*Y*Y+32.0D0*Y)+
1070     38.0D0*DL*DLP/Y/Y/5.0D0-8.0D0*DL/Y/5.0D0+
1071     4DL*S111MX*(-8.0D0+8.0D0*Y+16.0D0*DP)+
1072     5DL*S11MX*(-16.0D0+16.0D0*Y+32.0D0*DP)+
1073     6DL*(14.0D0/5.0D0+72.0D0*Y*Y/5.0D0-146.0D0*Y/5.0D0
1074     7-16.0D0*DP)+
1075     8ZETA2*DLM*(8.0D0-8.0D0*Y-16.0D0*DP)
1076      B3=DLM*S11MX*(16.0D0-16.0D0*Y-32.0D0*DP)+
1077     1DLM*(16.0D0-16.0D0*Y)+
1078     2ZETA2*DLP*(-48.0D0*Y-16.0D0*DP)+
1079     3DLP*S11MX*(-96.0D0*Y-32.0D0*DP)+
1080     4S21MZ*(16.0D0-16.0D0*Y-32.0D0*DP)+
1081     5S21Z*(-16.0D0+16.0D0*Y+32.0D0*DP)+
1082     6ZETA2*(-4.0D0-72.0D0*Y*Y*Y/5.0D0+48.0D0*Y*Y+20.0D0*Y)+
1083     7ZETA3*(48.0D0*Y+16.0D0*DP)
1084      B4=S111MX*(8.0D0+8.0D0*Y)+
1085     1S211MX*(16.0D0-16.0D0*Y-32.0D0*DP)+
1086     2S121MX*(-8.0D0+8.0D0*Y+16.0D0*DP)+
1087     3S11MX*(32.0D0-72.0D0*Y*Y*Y/5.0D0+48.0D0*Y*Y+32.0D0*Y)+
1088     48.0D0*S11MX/Y/Y/5.0D0+8.0D0/Y/5.0D0+
1089     5S21MX*(48.0D0*Y+16.0D0*DP)+
1090     6S12MX*(-96.0D0*Y-32.0D0*DP)+
1091     772.0D0*Y*Y/5.0D0-162.0D0/5.0D0+82.0D0*Y/5.0D0
1092      B5=DL**2*(-4.0D0+4.0D0*Y+8.0D0*DP)+
1093     1DL*DLP*(16.0D0-16.0D0*Y-32.0D0*DP)+
1094     2DL*(8.0D0+8.0D0*Y)+
1095     3ZETA2*(8.0D0-8.0D0*Y-16.0D0*DP)+
1096     4S11MX*(16.0D0-16.0D0*Y-32.0D0*DP)+
1097     516.0D0-16.0D0*Y
1098      COLOR=CF2-0.5D0*CACF
1099      C2NSM=COLOR*(B1+B2+B3+B4+B5*DL1)
1100      RETURN
1101      END
1102
1103C THIS IS THE WILSON COEFFICIENT CORRESPONDING TO THE PROCESS:
1104C V + Q -> Q + Q + QB WITH IDENTICAL QUARKS IN THE FINAL STATE
1105C (NONSINGLET) , F3(X,Q^2) , MSBAR SCHEME .
1106
1107
1108      DOUBLE PRECISION FUNCTION C3NSM(Y)
1109      IMPLICIT REAL*8(A-H,O-Z)
1110      COMPLEX*16 WGPLG
1111      DOUBLE PRECISION NF
1112      COMMON/PAR/NF,Q2,FS,RS,DN
1113      COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
1114      FS2=FS*FS
1115      DL1=DLOG(Q2/FS2)
1116      S121MX=WGPLG(1,2,1.0D0-Y)
1117      S12MX=WGPLG(1,2,-Y)
1118      S211MX=WGPLG(2,1,1.0D0-Y)
1119      S21MX=WGPLG(2,1,-Y)
1120      Z=(1.0D0-Y)/(1.0D0+Y)
1121      S21Z=WGPLG(2,1,Z)
1122      S21MZ=WGPLG(2,1,-Z)
1123      S111MX=WGPLG(1,1,1.0D0-Y)
1124      S11MX=WGPLG(1,1,-Y)
1125      ZETA2=1.6449340668D0
1126      ZETA3=1.2020569031D0
1127      DP=1.0D0/(1.0D0+Y)
1128      DL=DLOG(Y)
1129      DLP=DLOG(1.0D0+Y)
1130      DLM=DLOG(1.0D0-Y)
1131      B1=DL**3*(2.0D0-2.0D0*Y-4.0D0*DP)+
1132     1DL**2*DLM*(-4.0D0+4.0D0*Y+8.0D0*DP)+
1133     2DL**2*DLP*(-24.0D0+16.0D0*Y+40.0D0*DP)+
1134     3DL**2*(-8.0D0-4.0D0*Y*Y-4.0D0*Y)+
1135     4DL*DLM*DLP*(16.0D0-16.0D0*Y-32.0D0*DP)
1136      B2=DL*DLM*(8.0D0+8.0D0*Y)+
1137     1DL*DLP**2*(16.0D0-16.0D0*DP)+
1138     2DL*DLP*(8.0D0*Y*Y+8.0D0/Y)+
1139     3DL*S111MX*(-8.0D0+8.0D0*Y+16.0D0*DP)+
1140     4DL*S11MX*(-16.0D0+16.0D0*Y+32.0D0*DP)+
1141     5DL*(14.0D0+30.0D0*Y-16.0D0*DP)+
1142     6ZETA2*DLM*(8.0D0-8.0D0*Y-16.0D0*DP)
1143      B3=DLM*S11MX*(16.0D0-16.0D0*Y-32.0D0*DP)+
1144     1DLM*(16.0D0-16.0D0*Y)+
1145     2ZETA2*DLP*(16.0D0-16.0D0*DP)+
1146     3DLP*S11MX*(32.0D0-32.0D0*DP)+
1147     4S21MZ*(16.0D0-16.0D0*Y-32.0D0*DP)+
1148     5S21Z*(-16.0D0+16.0D0*Y+32.0D0*DP)+
1149     6ZETA2*(4.0D0+8.0D0*Y*Y-4.0D0*Y)+
1150     7ZETA3*(-16.0D0+16.0D0*DP)
1151      B4=S111MX*(8.0D0+8.0D0*Y)+
1152     1S211MX*(16.0D0-16.0D0*Y-32.0D0*DP)+
1153     2S121MX*(-8.0D0+8.0D0*Y+16.0D0*DP)+
1154     3S11MX*(8.0D0*Y*Y+8.0D0/Y)+
1155     4S21MX*(-16.0D0+16.0D0*DP)+
1156     5S12MX*(32.0D0-32.0D0*DP)+
1157     630.0D0-30.0D0*Y
1158      B5=DL**2*(-4.0D0+4.0D0*Y+8.0D0*DP)+
1159     1DL*DLP*(16.0D0-16.0D0*Y-32.0D0*DP)+
1160     2DL*(8.0D0+8.0D0*Y)+
1161     3ZETA2*(8.0D0-8.0D0*Y-16.0D0*DP)+
1162     4S11MX*(16.0D0-16.0D0*Y-32.0D0*DP)+
1163     516.0D0-16.0D0*Y
1164      C3NSM=(CF2-0.5D0*CACF)*(B1+B2+B3+B4+B5*DL1)
1165      RETURN
1166      END
1167
1168C THIS IS THE WILSON COEFFICIENT CORRESPONDING TO THE PROCESS:
1169C V + Q -> Q + Q + QB WITH IDENTICAL QUARKS IN THE FINAL STATE
1170C (NONSINGLET) , F3(X,Q^2) , DIS SCHEME .
1171
1172
1173      DOUBLE PRECISION FUNCTION D3NSM(Y)
1174      IMPLICIT REAL*8(A-H,O-Z)
1175      COMPLEX*16 WGPLG
1176      DOUBLE PRECISION NF
1177      COMMON/PAR/NF,Q2,FS,RS,DN
1178      COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
1179      FS2=FS*FS
1180      DL1=DLOG(Q2/FS2)
1181      S12MX=WGPLG(1,2,-Y)
1182      S21MX=WGPLG(2,1,-Y)
1183      S11MX=WGPLG(1,1,-Y)
1184      ZETA2=1.6449340668D0
1185      ZETA3=1.2020569031D0
1186      DP=1.0D0/(1.0D0+Y)
1187      DL=DLOG(Y)
1188      DLP=DLOG(1.0D0+Y)
1189      B1=DL**2*DLP*(-8.0D0-24.0D0*Y)+
1190     1DL**2*(-4.0D0-36.0D0*Y*Y*Y/5.0D0+20.0D0*Y*Y+12.0D0*Y)+
1191     2DL*DLP**2*(16.0D0+48.0D0*Y)+
1192     3DL*DLP*(-32.0D0+72.0D0*Y*Y*Y/5.0D0-40.0D0*Y*Y-32.0D0*Y+
1193     48.0D0/Y)-8.0D0*DL*DLP/Y/Y/5.0D0+8.0D0*DL/Y/5.0D0+
1194     5DL*(56.0D0/5.0D0-72.0D0*Y*Y/5.0D0+296.0D0*Y/5.0D0)+
1195     6ZETA2*DLP*(16.0D0+48.0D0*Y)
1196      B2=DLP*S11MX*(32.0D0+96.0D0*Y)+
1197     1ZETA2*(8.0D0+72.0D0*Y*Y*Y/5.0D0-40.0D0*Y*Y-24.0D0*Y)+
1198     2ZETA3*(-16.0D0-48.0D0*Y)+
1199     3S11MX*(-32.0D0+72.0D0*Y*Y*Y/5.0D0-40.0D0*Y*Y-32.0D0*Y+
1200     48.0D0/Y)-8.0D0*S11MX/Y/Y/5.0D0-8.0D0/Y/5.0D0+
1201     5S21MX*(-16.0D0-48.0D0*Y)+
1202     6S12MX*(32.0D0+96.0D0*Y)+
1203     7312.0D0/5.0D0-72.0D0*Y*Y/5.0D0-232.0D0*Y/5.0D0
1204      B3=DL**2*(-4.0D0+4.0D0*Y+8.0D0*DP)+
1205     1DL*DLP*(16.0D0-16.0D0*Y-32.0D0*DP)+
1206     2DL*(8.0D0+8.0D0*Y)+
1207     3ZETA2*(8.0D0-8.0D0*Y-16.0D0*DP)+
1208     4S11MX*(16.0D0-16.0D0*Y-32.0D0*DP)+
1209     516.0D0-16.0D0*Y
1210      D3NSM=(CF2-0.5D0*CACF)*(B1+B2+B3*DL1)
1211      RETURN
1212      END
1213
1214
1215C THIS IS THE WILSON COEFFICIENT CORRESPONDING TO THE PROCESS :
1216C V + Q -> Q + QB + Q (SINGLET) , F2(X,Q^2) , MSBAR SCHEME .
1217
1218
1219      DOUBLE PRECISION FUNCTION C2S(Y)
1220      IMPLICIT REAL*8(A-H,O-Z)
1221      COMPLEX*16 WGPLG
1222      DOUBLE PRECISION NF
1223      COMMON/PAR/NF,Q2,FS,RS,DN
1224      COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
1225      FS2=FS*FS
1226      DL1=DLOG(Q2/FS2)
1227      DL2=DL1*DL1
1228      S211MX=WGPLG(2,1,1.0D0-Y)
1229      S111MX=WGPLG(1,1,1.0D0-Y)
1230      S11MX=WGPLG(1,1,-Y)
1231      ZETA2=1.6449340668D0
1232      DL=DLOG(Y)
1233      DLP=DLOG(1.0D0+Y)
1234      DLM=DLOG(1.0D0-Y)
1235      A1=DL**3*(10.0D0/3.0D0+10.0D0*Y/3.0D0)+
1236     1DL**2*DLM*(-8.0D0-8.0D0*Y)+
1237     2DL**2*(-1.0D0-32.0D0*Y*Y/3.0D0+15.0D0*Y)+
1238     3DL*DLM**2*(4.0D0+4.0D0*Y)+
1239     4DL*DLM*(16.0D0*Y*Y)+
1240     5DL*DLP*(-16.0D0-16.0D0*Y*Y/3.0D0-16.0D0*Y-16.0D0/Y/3.0D0)
1241      A2=DL*ZETA2*(-8.0D0-8.0D0*Y)+
1242     1DL*S111MX*(-8.0D0-8.0D0*Y)+
1243     2DL*(56.0D0-128.0D0*Y*Y/9.0D0-88.0D0*Y/3.0D0)+
1244     3DLM**2*(2.0D0-8.0D0*Y*Y/3.0D0-2.0D0*Y+8.0D0/Y/3.0D0)+
1245     5DLM*S111MX*(8.0D0+8.0D0*Y)+
1246     6DLM*(-104.0D0/3.0D0-32.0D0*Y*Y/9.0D0+80.0D0*Y/3.0D0
1247     7+104.0D0/Y/9.0D0)
1248      A3=ZETA2*(-4.0D0+16.0D0*Y*Y/3.0D0-12.0D0*Y-32.0D0/Y/3.0D0)+
1249     1S111MX*(4.0D0+32.0D0*Y*Y/3.0D0-4.0D0*Y+16.0D0/Y/3.0D0)+
1250     2S211MX*(-8.0D0-8.0D0*Y)+
1251     3S11MX*(-16.0D0-16.0D0*Y*Y/3.0D0-16.0D0*Y-16.0D0/Y/3.0D0)+
1252     4158.0D0/9.0D0+448.0D0*Y*Y/27.0D0-422.0D0*Y/9.0D0+
1253     5344.0D0/Y/27.0D0
1254C      A3=A3+32.0D0*ZETA2/Y/3.0D0-344.0D0/Y/27.0D0
1255      A4=DL*(4.0D0+4.0D0*Y)+
1256     12.0D0-2.0D0*Y-8.0D0*Y*Y/3.0D0+8.0D0/Y/3.0D0
1257      A5=DL**2*(-8.0D0-8.0D0*Y)+
1258     1DL*DLM*(8.0D0+8.0D0*Y)+
1259     2DL*(16.0D0*Y*Y)+
1260     3DLM*(4.0D0-4.0D0*Y-16.0D0*Y*Y/3.0D0+16.0D0/Y/3.0D0)+
1261     4S111MX*(8.0D0+8.0D0*Y)-104.0D0/3.0D0+80.0D0*Y/3.0D0-
1262     532.0D0*Y*Y/9.0D0+104.0D0/Y/9.0D0
1263      C2S=2.0D0*NF*CFTF*(A1+A2+A3+A4*DL2+A5*DL1)
1264      RETURN
1265      END
1266
1267
1268C THIS IS THE WILSON COEFFICIENT CORRESPONDING TO THE PROCESS :
1269C V + G -> Q + QB + G , F2(X,Q^2) , MSBAR SCHEME .
1270
1271
1272      DOUBLE PRECISION FUNCTION C2G2(Y)
1273      IMPLICIT REAL*8(A-H,O-Z)
1274      COMPLEX*16 WGPLG
1275      DOUBLE PRECISION NF
1276      COMMON/PAR/NF,Q2,FS,RS,DN
1277      COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
1278      FS2=FS*FS
1279      RS2=RS*RS
1280      DL1=DLOG(Q2/FS2)
1281      DL2=DL1*DL1
1282      EL1=DLOG(FS2/RS2)
1283      S121MX=WGPLG(1,2,1.0D0-Y)
1284      S12MX=WGPLG(1,2,-Y)
1285      S211MX=WGPLG(2,1,1.0D0-Y)
1286      S21MX=WGPLG(2,1,-Y)
1287      Z=(1.0D0-Y)/(1.0D0+Y)
1288      S21Z=WGPLG(2,1,Z)
1289      S21MZ=WGPLG(2,1,-Z)
1290      S111MX=WGPLG(1,1,1.0D0-Y)
1291      S11MX=WGPLG(1,1,-Y)
1292      ZETA2=1.6449340668D0
1293      ZETA3=1.2020569031D0
1294      DL=DLOG(Y)
1295      DLP=DLOG(1.0D0+Y)
1296      DLM=DLOG(1.0D0-Y)
1297      A1=DL**3*(-5.0D0/3.0D0-20.0D0*Y*Y/3.0D0+10.0D0*Y/3.0D0)+
1298     1DL**2*DLM*(8.0D0+24.0D0*Y*Y-16.0D0*Y)+
1299     2DL**2*DLP*(8.0D0+8.0D0*Y*Y+16.0D0*Y)+
1300     3DL**2*(-1.5D0-48.0D0*Y*Y*Y/5.0D0-36.0D0*Y*Y+22.0D0*Y/3.0D0)+
1301     4DL*DLM**2*(-8.0D0-20.0D0*Y*Y+16.0D0*Y)+
1302     5DL*DLM*(16.0D0+72.0D0*Y*Y-56.0D0*Y)+
1303     6DL*DLP**2*(-16.0D0-16.0D0*Y*Y-32.0D0*Y)+
1304     7DL*DLP*(48.0D0+96.0D0*Y*Y*Y/5.0D0+64.0D0*Y/3.0D0)
1305      A2=8.0D0*DL*DLP/Y/Y/15.0D0-8.0D0*DL/Y/15.0D0+
1306     1DL*ZETA2*(12.0D0+32.0D0*Y*Y-24.0D0*Y)+
1307     2DL*S111MX*(4.0D0+16.0D0*Y*Y-8.0D0*Y)+
1308     3DL*S11MX*(-16.0D0-16.0D0*Y*Y+32.0D0*Y)+
1309     4DL*(-236.0D0/15.0D0-216.0D0*Y*Y/5.0D0+113.0D0*Y/5.0D0)+
1310     5DLM**3*(10.0D0/3.0D0+20.0D0*Y*Y/3.0D0-20.0D0*Y/3.0D0)+
1311     6DLM**2*(-13.0D0-36.0D0*Y*Y+40.0D0*Y)+
1312     7ZETA2*DLM*(-16.0D0*Y*Y)
1313      A3=DLM*S111MX*(8.0D0+8.0D0*Y*Y-16.0D0*Y)+
1314     1DLM*(14.0D0+24.0D0*Y*Y-40.0D0*Y)+
1315     2ZETA2*DLP*(-16.0D0-16.0D0*Y*Y-32.0D0*Y)+
1316     3DLP*S11MX*(-32.0D0-32.0D0*Y*Y-64.0D0*Y)+
1317     4ZETA2*(26.0D0+96.0D0*Y*Y*Y/5.0D0+72.0D0*Y*Y-176.0D0*Y/3.0D0)+
1318     5ZETA3*(52.0D0+88.0D0*Y*Y-40.0D0*Y)
1319      A4=S111MX*(-10.0D0+24.0D0*Y)+
1320     1S211MX*(-16.0D0-24.0D0*Y*Y+32.0D0*Y)+
1321     2S121MX*(-4.0D0+8.0D0*Y*Y+8.0D0*Y)+
1322     3S11MX*(48.0D0+96.0D0*Y*Y*Y/5.0D0+64.0D0*Y/3.0D0)+
1323     48.0D0*S11MX/Y/Y/15.0D0+8.0D0/Y/15.0D0+
1324     5S21MX*(48.0D0+48.0D0*Y*Y-32.0D0*Y)+
1325     6S12MX*(-32.0D0-32.0D0*Y*Y-64.0D0*Y)-
1326     7647.0D0/15.0D0-36.0D0*Y*Y/5.0D0+239.0D0*Y/5.0D0
1327      A5=DL*(-2.0D0+4.0D0*Y-8.0D0*Y*Y)+
1328     1DLM*(4.0D0-8.0D0*Y+8.0D0*Y*Y)-
1329     21.0D0+4.0D0*Y
1330      A6=DL**2*(4.0D0-8.0D0*Y+16.0D0*Y*Y)+
1331     1DL*DLM*(-12.0D0+24.0D0*Y-32.0D0*Y*Y)+
1332     2DL*(4.0D0-24.0D0*Y+40.0D0*Y*Y)+
1333     3DLM**2*(8.0D0-16.0D0*Y+16.0D0*Y*Y)+
1334     4DLM*(-14.0D0+48.0D0*Y-40.0D0*Y*Y)+
1335     5ZETA2*(-16.0D0+32.0D0*Y-32.0D0*Y*Y)+
1336     6S111MX*(4.0D0-8.0D0*Y)+
1337     718.0D0-34.0D0*Y+8.0D0*Y*Y
1338      FCF=CF*(A1+A2+A3+A4+A5*DL2+A6*DL1)
1339      B1=DL**3*(10.0D0/3.0D0+28.0D0*Y/3.0D0)+
1340     1DL**2*DLM*(-4.0D0+8.0D0*Y*Y-32.0D0*Y)+
1341     2DL**2*DLP*(4.0D0+12.0D0*Y*Y+8.0D0*Y)+
1342     3DL**2*(-1.0D0-194.0D0*Y*Y/3.0D0+88.0D0*Y)+
1343     4DL*DLM**2*(-8.0D0*Y*Y+24.0D0*Y)+
1344     5DL*DLM*DLP*(-8.0D0-16.0D0*Y*Y-16.0D0*Y)+
1345     6DL*DLM*(8.0D0+148.0D0*Y*Y-144.0D0*Y)+
1346     8DL*DLP**2*(8.0D0+8.0D0*Y*Y+16.0D0*Y)
1347      B2=DL*DLP*(-24.0D0+80.0D0*Y*Y/3.0D0-16.0D0/Y/3.0D0)+
1348     1DL*ZETA2*(16.0D0*Y*Y-48.0D0*Y)+
1349     2DL*S111MX*(-8.0D0-16.0D0*Y)+
1350     3DL*S11MX*(8.0D0+8.0D0*Y*Y+16.0D0*Y)+
1351     4DL*(58.0D0-2090.0D0*Y*Y/9.0D0+584.0D0*Y/3.0D0)+
1352     5DLM**3*(2.0D0/3.0D0+4.0D0*Y*Y/3.0D0-4.0D0*Y/3.0D0)+
1353     6DLM**2*(-2.0D0-122.0D0*Y*Y/3.0D0+36.0D0*Y+8.0D0/Y/3.0D0)+
1354     7ZETA2*DLM*(-20.0D0-32.0D0*Y*Y+24.0D0*Y)
1355      B3=DLM*S111MX*(4.0D0-8.0D0*Y*Y+40.0D0*Y)+
1356     1DLM*S11MX*(-8.0D0-16.0D0*Y*Y-16.0D0*Y)+
1357     2DLM*(-62.0D0/3.0D0+1570.0D0*Y*Y/9.0D0-454.0D0*Y/3.0D0
1358     3+104.0D0/Y/9.0D0)+
1359     4ZETA2*DLP*(8.0D0+8.0D0*Y*Y+16.0D0*Y)+
1360     5DLP*S11MX*(16.0D0+16.0D0*Y*Y+32.0D0*Y)+
1361     6S21MZ*(-8.0D0-16.0D0*Y*Y-16.0D0*Y)+
1362     7S21Z*(8.0D0+16.0D0*Y*Y+16.0D0*Y)
1363      B4=ZETA2*(4.0D0+268.0D0*Y*Y/3.0D0-80.0D0*Y-32.0D0/Y/3.0D0)+
1364     1ZETA3*(-10.0D0-12.0D0*Y*Y-12.0D0*Y)+
1365     2S111MX*(4.0D0+176.0D0*Y*Y/3.0D0-64.0D0*Y+16.0D0/Y/3.0D0)+
1366     3S211MX*(-4.0D0+8.0D0*Y*Y-72.0D0*Y)+
1367     4S121MX*(4.0D0+8.0D0*Y)+
1368     5S11MX*(-24.0D0+80.0D0*Y*Y/3.0D0-16.0D0/Y/3.0D0)+
1369     6S21MX*(-8.0D0+8.0D0*Y*Y-16.0D0*Y)+
1370     7S12MX*(16.0D0+16.0D0*Y*Y+32.0D0*Y)+344.0D0/Y/27.0D0+
1371     8239.0D0/9.0D0-4493.0D0*Y*Y/27.0D0+1072.0D0*Y/9.0D0
1372C      B4=B4+16.0D0*ZETA2/Y/3.0D0-344.0D0/Y/27.0D0
1373      B5=DL*(4.0D0+16.0D0*Y)+
1374     1DLM*(4.0D0-8.0D0*Y+8.0D0*Y*Y)+2.0D0+16.0D0*Y-
1375     262.0D0*Y*Y/3.0D0+8.0D0/Y/3.0D0
1376      B6=-22.0D0/3.0D0+44.0D0*Y/3.0D0-44.0D0*Y*Y/3.0D0
1377      B7=DL**2*(-8.0D0-24.0D0*Y)+
1378     1DL*DLM*(48.0D0*Y-16.0D0*Y*Y)+
1379     2DL*DLP*(-8.0D0-16.0D0*Y-16.0D0*Y*Y)+
1380     3DL*(-96.0D0*Y+100.0D0*Y*Y)+
1381     4DLM**2*(4.0D0-8.0D0*Y+8.0D0*Y*Y)+
1382     5DLM*(-4.0D0+80.0D0*Y-268.0D0*Y*Y/3.0D0+16.0D0/Y/3.0D0)+
1383     6ZETA2*(-8.0D0-16.0D0*Y*Y)+
1384     7S111MX*(8.0D0+32.0D0*Y)+
1385     8S11MX*(-8.0D0-16.0D0*Y-16.0D0*Y*Y)-110.0D0/3.0D0-
1386     9184.0D0*Y/3.0D0+814.0D0*Y*Y/9.0D0+104.0D0/Y/9.0D0
1387      B8=22.0D0*(1.0D0-2.0D0*Y+2.0D0*Y*Y)*(DL-DLM)/3.0D0+
1388     122.0D0/3.0D0-176.0D0*Y/3.0D0+176.0D0*Y*Y/3.0D0
1389      FCA=CA*(B1+B2+B3+B4+B5*DL2+B6*DL1*EL1+B7*DL1+B8*EL1)
1390      C1=4.0D0/3.0D0-8.0D0*Y/3.0D0+8.0D0*Y*Y/3.0D0
1391      C2=4.0D0*(1.0D0-2.0D0*Y+2.0D0*Y*Y)*(DLM-DL)/3.0D0-
1392     14.0D0/3.0D0+32.0D0*Y/3.0D0-32.0D0*Y*Y/3.0D0
1393      FNF=NF*(C1*DL1*EL1+C2*EL1)
1394      C2G2=NF*(FCF+FCA+FNF)
1395      RETURN
1396      END
1397
1398C THIS IS THE LONGITUDINAL WILSON COEFFICIENT CORRESPONDING
1399C TO THE PROCESS : V + Q -> Q + G  (SINGLET AND NONSINGLET)
1400C FL(X,Q^2 , MSBAR SCHEME .
1401
1402      DOUBLE PRECISION FUNCTION CLNS1(Y)
1403      IMPLICIT REAL*8(A-H,O-Z)
1404      COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
1405      A1=4.0D0*Y
1406      CLNS1=CF*A1
1407      RETURN
1408      END
1409
1410
1411C THIS IS THE LONGITUDINAL WILSON COEFFICIENT CORRESPONDING
1412C TO THE PROCESS : V + G -> Q + QB
1413C FL(X,Q^2) , MSBAR SCHEME .
1414
1415      DOUBLE PRECISION FUNCTION CLG1(Y)
1416      IMPLICIT REAL*8(A-H,O-Z)
1417      DOUBLE PRECISION NF
1418      COMMON/PAR/NF,Q2,FS,RS,DN
1419      A1=8.0D0*Y-8.0D0*Y*Y
1420      CLG1=NF*A1
1421      RETURN
1422      END
1423
1424
1425C THIS IS THE LONGITUDINAL WILSON COEFFICIENT CORRESPONDING TO THE PROCESSES :
1426C V + Q -> Q + G + G ; V + Q -> Q + Q + QB WITH NONIDENTICAL QUARKS
1427C IN THE FINAL STATE (NONSINGLET) , FL(X,Q^2) , MSBAR SCHEME .
1428
1429
1430
1431      DOUBLE PRECISION FUNCTION CLNSP(Y)
1432      IMPLICIT REAL*8(A-H,O-Z)
1433      COMPLEX*16 WGPLG
1434      DOUBLE PRECISION NF
1435      COMMON/PAR/NF,Q2,FS,RS,DN
1436      COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
1437      FS2=FS*FS
1438      RS2=RS*RS
1439      DL1=DLOG(Q2/FS2)
1440      CL1=DLOG(Q2/RS2)
1441      S121MX=WGPLG(1,2,1.0D0-Y)
1442      S21MX=WGPLG(2,1,-Y)
1443      S111MX=WGPLG(1,1,1.0D0-Y)
1444      S11MX=WGPLG(1,1,-Y)
1445      ZETA2=1.6449340668D0
1446      DL=DLOG(Y)
1447      DLP=DLOG(1.0D0+Y)
1448      DLM=DLOG(1.0D0-Y)
1449      A1=DL**2*(24.0D0*Y*Y*Y/5.0D0+16.0D0*Y*Y)+
1450     1DL*DLM*(-24.0D0*Y)+
1451     2DL*DLP*(32.0D0-48.0D0*Y*Y*Y/5.0D0-32.0D0*Y*Y+16.0D0/Y)+
1452     332.0D0*DL*DLP/Y/Y/5.0D0-32.0D0*DL/Y/5.0D0+
1453     4DL*S11MX*(32.0D0*Y)+
1454     5DL*(-24.0D0/5.0D0+48.0D0*Y*Y/5.0D0+16.0D0*Y/5.0D0)+
1455     6DLM**2*(8.0D0*Y)+
1456     7ZETA2*DLM*(-32.0D0*Y)
1457      A2=DLM*(8.0D0+28.0D0*Y)+
1458     1ZETA2*(-48.0D0*Y*Y*Y/5.0D0-32.0D0*Y*Y)+
1459     2S111MX*(-8.0D0*Y)+
1460     3S121MX*(32.0D0*Y)+
1461     4S11MX*(32.0D0-48.0D0*Y*Y*Y/5.0D0-32.0D0*Y*Y+16.0D0/Y)+
1462     532.0D0*S11MX/Y/Y/5.0D0+32.0D0/Y/5.0D0+
1463     6S21MX*(-64.0D0*Y)+252.0D0/5.0D0+
1464     748.0D0*Y*Y/5.0D0-502.0D0*Y/5.0D0
1465      A3=DL*(-8.0D0*Y)+
1466     1DLM*(16.0D0*Y)+
1467     28.0D0+4.0D0*Y
1468      FCF=CF2*(A1+A2+A3*DL1)
1469      B1=DL**2*(-12.0D0*Y*Y*Y/5.0D0-8.0D0*Y*Y+4.0D0*Y)+
1470     1DL*DLP*(-16.0D0+24.0D0*Y*Y*Y/5.0D0+16.0D0*Y*Y-8.0D0/Y)
1471     2-16.0D0*DL*DLP/Y/Y/5.0D0+16.0D0*DL/Y/5.0D0+
1472     3DL*S11MX*(-16.0D0*Y)+
1473     4DL*(-8.0D0/5.0D0-24.0D0*Y*Y/5.0D0+416.0D0*Y/15.0D0)+
1474     5ZETA2*DLM*(16.0D0*Y)+
1475     6DLM*(-92.0D0*Y/3.0D0)
1476      B2=ZETA2*(24.0D0*Y*Y*Y/5.0D0+16.0D0*Y*Y-8.0D0*Y)+
1477     1S121MX*(-16.0D0*Y)+
1478     2S11MX*(-16.0D0+24.0D0*Y*Y*Y/5.0D0+16.0D0*Y*Y-8.0D0/Y)
1479     3-16.0D0*S11MX/Y/Y/5.0D0-16.0D0/Y/5.0D0+
1480     4S21MX*(32.0D0*Y)-508.0D0/15.0D0
1481     5-24.0D0*Y*Y/5.0D0+4034.0D0*Y/45.0D0
1482      B3=-44.0D0*Y/3.0D0
1483      FCA=CACF*(B1+B2+B3*CL1)
1484      C1=DL*(-16.0D0*Y/3.0)+
1485     1DLM*(8.0D0*Y/3.0D0)+
1486     28.0D0/3.0D0-100.0D0*Y/9.0D0
1487      C2=8.0D0*Y/3.0D0
1488      FNF=2.0D0*NF*CFTF*(C1+C2*CL1)
1489      CLNSP=FCF+FCA+FNF
1490      RETURN
1491      END
1492
1493C THIS IS THE LONGITUDINAL WILSON COEFFICIENT CORRESPONDING TO THE PROCESS:
1494C V + Q -> Q + Q + QB WITH IDENTICAL QUARKS IN THE FINAL STATE
1495C (NONSINGLET) , FL(X,Q^2) ,  MSBAR AND DIS SCHEME .
1496
1497
1498      DOUBLE PRECISION FUNCTION CLNSM(Y)
1499      IMPLICIT REAL*8(A-H,O-Z)
1500      COMPLEX*16 WGPLG
1501      COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
1502      S12MX=WGPLG(1,2,-Y)
1503      S21MX=WGPLG(2,1,-Y)
1504      S11MX=WGPLG(1,1,-Y)
1505      ZETA2=1.6449340668D0
1506      ZETA3=1.2020569031D0
1507      DL=DLOG(Y)
1508      DLP=DLOG(1.0D0+Y)
1509      A1=DL**2*DLP*(16.0D0*Y)+
1510     1DL**2*(24.0D0*Y*Y*Y/5.0D0-16.0D0*Y*Y-8.0D0*Y)+
1511     2DL*DLP**2*(-32.0D0*Y)+
1512     3DL*DLP*(32.0D0-48.0D0*Y*Y*Y/5.0D0+32.0D0*Y*Y+32.0D0*Y-16.0D0/Y)+
1513     432.0D0*DL*DLP/Y/Y/5.0D0-32.0D0*DL/Y/5.0D0+
1514     5DL*(16.0D0/5.0D0+48.0D0*Y*Y/5.0D0-224.0D0*Y/5.0D0)+
1515     6ZETA2*DLP*(-32.0D0*Y)+
1516     7DLP*S11MX*(-64.0D0*Y)
1517      A2=ZETA2*(-48.0D0*Y*Y*Y/5.0D0+32.0D0*Y*Y+16.0D0*Y)+
1518     1ZETA3*(32.0D0*Y)+
1519     2S11MX*(32.0D0-48.0D0*Y*Y*Y/5.0D0+32.0D0*Y*Y+32.0D0*Y-16.0D0/Y)+
1520     332.0D0*S11MX/Y/Y/5.0D0+32.0D0/Y/5.0D0+
1521     4S21MX*(32.0D0*Y)+
1522     5S12MX*(-64.0D0*Y)-208.0D0/5.0D0+
1523     648.0D0*Y*Y/5.0D0+128.0D0*Y/5.0D0
1524      CLNSM=(CF2-0.5D0*CACF)*(A1+A2)
1525      RETURN
1526      END
1527
1528C THIS IS THE LONGITUDINAL WILSON COEFFICIENT CORRESPONDING TO THE PROCESS :
1529C V + Q -> Q + QB + Q (SINGLET) , FL(X,Q^2) ,  MSBAR SCHEME .
1530
1531
1532      DOUBLE PRECISION FUNCTION CLS(Y)
1533      IMPLICIT REAL*8(A-H,O-Z)
1534      COMPLEX*16 WGPLG
1535      DOUBLE PRECISION NF
1536      COMMON/PAR/NF,Q2,FS,RS,DN
1537      COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
1538      S111MX=WGPLG(1,1,1.0D0-Y)
1539      DL=DLOG(Y)
1540      DLM=DLOG(1.0D0-Y)
1541      FS2=FS*FS
1542      DL1=DLOG(Q2/FS2)
1543      A1=DL**2*(16.0D0*Y)+
1544     1DL*DLM*(-16.0D0*Y)+
1545     2DL*(16.0D0-32.0D0*Y*Y-16.0D0*Y)+
1546     3DLM*(-16.0D0+32.0D0*Y*Y/3.0D0+16.0D0/Y/3.0D0)+
1547     4S111MX*(-16.0D0*Y)+16.0D0/3.0D0+
1548     5160.0D0*Y*Y/9.0D0-64.0D0*Y/3.0D0-16.0D0/Y/9.0D0
1549      A2=DL*(-16.0D0*Y)-
1550     116.0D0+32.0D0*Y*Y/3.0D0+16.0D0/Y/3.0D0
1551      CLS=2.0D0*NF*CFTF*(A1+A2*DL1)
1552      RETURN
1553      END
1554
1555C THIS IS THE LONGITUDINAL WILSON COEFFICIENT CORRESPONDING TO THE PROCESS :
1556C V + G -> Q + QB + G , FL(X,Q^2) , MSBAR SCHEME .
1557
1558
1559      DOUBLE PRECISION FUNCTION CLG2(Y)
1560      IMPLICIT REAL*8(A-H,O-Z)
1561      COMPLEX*16 WGPLG
1562      DOUBLE PRECISION NF
1563      COMMON/PAR/NF,Q2,FS,RS,DN
1564      COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
1565      FS2=FS*FS
1566      RS2=RS*RS
1567      DL1=DLOG(Q2/FS2)
1568      EL1=DLOG(FS2/RS2)
1569      S111MX=WGPLG(1,1,1.0D0-Y)
1570      S11MX=WGPLG(1,1,-Y)
1571      ZETA2=1.6449340668D0
1572      DL=DLOG(Y)
1573      DLP=DLOG(1.0D0+Y)
1574      DLM=DLOG(1.0D0-Y)
1575      A1=DL**2*(-32.0D0*Y*Y*Y/5.0D0-32.0D0*Y/3.0D0)+
1576     1DL*DLM*(16.0D0*Y)+
1577     2DL*DLP*(64.0D0*Y*Y*Y/5.0D0-32.0D0*Y/3.0D0)+
1578     332.0D0*DL*DLP/Y/Y/15.0D0-32.0D0*DL/Y/15.0D0+
1579     4DL*(-104.0D0/15.0D0+96.0D0*Y*Y/5.0D0-208.0D0*Y/5.0D0)+
1580     5DLM*(8.0D0-32.0D0*Y*Y+24.0D0*Y)
1581      A2=ZETA2*(64.0D0*Y*Y*Y/5.0D0-32.0D0*Y/3.0D0)+
1582     1S111MX*(16.0D0*Y)+
1583     2S11MX*(64.0D0*Y*Y*Y/5.0D0-32.0D0*Y/3.0D0)+
1584     332.0D0*S11MX/Y/Y/15.0D0+32.0D0/Y/15.0D0
1585     4-128.0D0/15.0D0+336.0D0*Y*Y/5.0D0-304.0D0*Y/5.0D0
1586      A3=DL*(16.0D0*Y)+
1587     18.0D0+8.0D0*Y-16.0D0*Y*Y
1588      FCF=CF*(A1+A2+A3*DL1)
1589      B1=DL**2*(48.0D0*Y)+
1590     1DL*DLM*(32.0D0*Y*Y-96.0D0*Y)+
1591     2DL*DLP*(32.0D0*Y*Y+32.0D0*Y)+
1592     3DL*(16.0D0-208.0D0*Y*Y+128.0D0*Y)+
1593     4DLM**2*(-16.0D0*Y*Y+16.0D0*Y)+
1594     5DLM*(-16.0D0+464.0D0*Y*Y/3.0D0-144.0D0*Y+16.0D0/Y/3.0D0)
1595      B2=ZETA2*(32.0D0*Y*Y)+
1596     1S111MX*(-64.0D0*Y)+
1597     2S11MX*(32.0D0*Y*Y+32.0D0*Y)+16.0D0/3.0D0
1598     3-848.0D0*Y*Y/9.0D0+272.0D0*Y/3.0D0-16.0D0/Y/9.0D0
1599      B3=DL*(-64.0D0*Y)+
1600     1DLM*(32.0D0*Y-32.0D0*Y*Y)-
1601     216.0D0-80.0D0*Y+272.0D0*Y*Y/3.0D0+16.0D0/Y/3.0D0
1602      B4=-88.0D0*Y/3.0D0+88.0D0*Y*Y/3.0D0
1603      FCA=CA*(B1+B2+B3*DL1+B4*EL1)
1604      C1=16.0D0*Y/3.0D0-16.0D0*Y*Y/3.0D0
1605      FNF=NF*C1*EL1
1606      CLG2=NF*(FCF+FCA+FNF)
1607      RETURN
1608      END
1609
1610C THIS IS THE LONGITUDINAL WILSON COEFFICIENT CORRESPONDING TO THE PROCESSES :
1611C V + Q -> Q + G + G ; V + Q -> Q + Q + QB WITH NONIDENTICAL QUARKS
1612C IN THE FINAL STATE (NONSINGLET) , FL(X,Q^2) , DIS SCHEME .
1613
1614
1615
1616      DOUBLE PRECISION FUNCTION DLNSP(Y)
1617      IMPLICIT REAL*8(A-H,O-Z)
1618      COMPLEX*16 WGPLG
1619      DOUBLE PRECISION NF
1620      COMMON/PAR/NF,Q2,FS,RS,DN
1621      COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
1622      FS2=FS*FS
1623      RS2=RS*RS
1624      DL1=DLOG(Q2/FS2)
1625      CL1=DLOG(Q2/RS2)
1626      S121MX=WGPLG(1,2,1.0D0-Y)
1627      S21MX=WGPLG(2,1,-Y)
1628      S111MX=WGPLG(1,1,1.0D0-Y)
1629      S11MX=WGPLG(1,1,-Y)
1630      ZETA2=1.6449340668D0
1631      DL=DLOG(Y)
1632      DLP=DLOG(1.0D0+Y)
1633      DLM=DLOG(1.0D0-Y)
1634      A1=DL**2*(24.0D0*Y*Y*Y/5.0D0+16.0D0*Y*Y-4.0D0*Y)+
1635     1DL*DLM*(-16.0D0*Y)+
1636     2DL*DLP*(32.0D0-48.0D0*Y*Y*Y/5.0D0-32.0D0*Y*Y+16.0D0/Y)+
1637     332.0D0*DL*DLP/Y/Y/5.0D0-32.0D0*DL/Y/5.0D0+
1638     4DL*S11MX*(32.0D0*Y)+
1639     5DL*(16.0D0/5.0D0+48.0D0*Y*Y/5.0D0-4.0D0*Y/5.0D0)+
1640     7ZETA2*DLM*(-32.0D0*Y)
1641      A2=DLM*(48.0D0*Y)+
1642     1ZETA2*(-48.0D0*Y*Y*Y/5.0D0-32.0D0*Y*Y+16.0D0*Y)+
1643     2S111MX*(-16.0D0*Y)+
1644     3S121MX*(32.0D0*Y)+
1645     4S11MX*(32.0D0-48.0D0*Y*Y*Y/5.0D0-32.0D0*Y*Y+16.0D0/Y)+
1646     532.0D0*S11MX/Y/Y/5.0D0+32.0D0/Y/5.0D0+
1647     6S21MX*(-64.0D0*Y)+232.0D0/5.0D0+
1648     748.0D0*Y*Y/5.0D0-302.0D0*Y/5.0D0
1649      A3=DL*(-8.0D0*Y)+
1650     1DLM*(16.0D0*Y)+
1651     28.0D0+4.0D0*Y
1652      FCF=CF2*(A1+A2+A3*DL1)
1653      B1=DL**2*(-12.0D0*Y*Y*Y/5.0D0-8.0D0*Y*Y+4.0D0*Y)+
1654     1DL*DLP*(-16.0D0+24.0D0*Y*Y*Y/5.0D0+16.0D0*Y*Y-8.0D0/Y)
1655     2-16.0D0*DL*DLP/Y/Y/5.0D0+16.0D0*DL/Y/5.0D0+
1656     3DL*S11MX*(-16.0D0*Y)+
1657     4DL*(-8.0D0/5.0D0-24.0D0*Y*Y/5.0D0+416.0D0*Y/15.0D0)+
1658     5ZETA2*DLM*(16.0D0*Y)+
1659     6DLM*(-92.0D0*Y/3.0D0)
1660      B2=ZETA2*(24.0D0*Y*Y*Y/5.0D0+16.0D0*Y*Y-8.0D0*Y)+
1661     1S121MX*(-16.0D0*Y)+
1662     2S11MX*(-16.0D0+24.0D0*Y*Y*Y/5.0D0+16.0D0*Y*Y-8.0D0/Y)
1663     3-16.0D0*S11MX/Y/Y/5.0D0-16.0D0/Y/5.0D0+
1664     4S21MX*(32.0D0*Y)-508.0D0/15.0D0
1665     5-24.0D0*Y*Y/5.0D0+4034.0D0*Y/45.0D0
1666      B3=-44.0D0*Y/3.0D0
1667      FCA=CACF*(B1+B2+B3*CL1)
1668      C1=DL*(-16.0D0*Y/3.0)+
1669     1DLM*(8.0D0*Y/3.0D0)+
1670     28.0D0/3.0D0-100.0D0*Y/9.0D0
1671      C2=8.0D0*Y/3.0D0
1672      FNF=2.0D0*NF*CFTF*(C1+C2*CL1)
1673      DLNSP=FCF+FCA+FNF
1674      RETURN
1675      END
1676
1677
1678C THIS IS THE LONGITUDINAL WILSON COEFFICIENT CORRESPONDING TO THE PROCESS :
1679C V + Q -> Q + QB + Q (SINGLET) , FL(X,Q^2)  , DIS SCHEME .
1680
1681
1682
1683      DOUBLE PRECISION FUNCTION DLS(Y)
1684      IMPLICIT REAL*8(A-H,O-Z)
1685      COMPLEX*16 WGPLG
1686      DOUBLE PRECISION NF
1687      COMMON/PAR/NF,Q2,FS,RS,DN
1688      COMMON/DIS/NDIS
1689      COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
1690      S111MX=WGPLG(1,1,1.0D0-Y)
1691      ZETA2=1.6449340668D0
1692      DL=DLOG(Y)
1693      DLM=DLOG(1.0D0-Y)
1694      FS2=FS*FS
1695      DL1=DLOG(Q2/FS2)
1696      IF (NDIS.EQ.1) THEN
1697      A1=DL**2*(-16.0D0*Y*Y+24.0D0*Y)+
1698     1DL*DLM*(32.0D0*Y*Y-32.0D0*Y)+
1699     2DL*(8.0D0-80.0D0*Y*Y+8.0D0*Y)+
1700     3DLM**2*(-16.0D0*Y*Y+16.0D0*Y)+
1701     4DLM*(-8.0D0+176.0D0*Y*Y/3.0D0-56.0D0*Y+16.0D0/Y/3.0D0)+
1702     5ZETA2*(32.0D0*Y*Y-32.0D0*Y)+16.0D0/3.0D0+
1703     5736.0D0*Y*Y/9.0D0-256.0D0*Y/3.0D0-16.0D0/Y/9.0D0
1704      A2=DL*(-16.0D0*Y)-
1705     116.0D0+32.0D0*Y*Y/3.0D0+16.0D0/Y/3.0D0
1706      ELSE
1707      A1=DL**2*(16.0D0*Y)+
1708     1DL*DLM*(-16.0D0*Y)+
1709     2DL*(16.0D0-32.0D0*Y*Y-16.0D0*Y)+
1710     3DLM*(-16.0D0+32.0D0*Y*Y/3.0D0+16.0D0/Y/3.0D0)+
1711     4S111MX*(-16.0D0*Y)+16.0D0/3.0D0+
1712     5160.0D0*Y*Y/9.0D0-64.0D0*Y/3.0D0-16.0D0/Y/9.0D0
1713      A2=DL*(-16.0D0*Y)-
1714     116.0D0+32.0D0*Y*Y/3.0D0+16.0D0/Y/3.0D0
1715      ENDIF
1716      DLS=2.0D0*NF*CFTF*(A1+A2*DL1)
1717      RETURN
1718      END
1719
1720
1721
1722
1723C THIS IS THE LONGITUDINAL WILSON COEFFICIENT CORRESPONDING TO THE PROCESS :
1724C V + G -> Q + QB + G , FL(X,Q^2) , DIS SCHEME .
1725
1726
1727      DOUBLE PRECISION FUNCTION DLG2(Y)
1728      IMPLICIT REAL*8(A-H,O-Z)
1729      COMPLEX*16 WGPLG
1730      DOUBLE PRECISION NF
1731      COMMON/PAR/NF,Q2,FS,RS,DN
1732      COMMON/DIS/NDIS
1733      COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
1734      FS2=FS*FS
1735      RS2=RS*RS
1736      DL1=DLOG(Q2/FS2)
1737      EL1=DLOG(FS2/RS2)
1738      S111MX=WGPLG(1,1,1.0D0-Y)
1739      S11MX=WGPLG(1,1,-Y)
1740      ZETA2=1.6449340668D0
1741      DL=DLOG(Y)
1742      DLP=DLOG(1.0D0+Y)
1743      DLM=DLOG(1.0D0-Y)
1744      A1=DL**2*(-32.0D0*Y*Y*Y/5.0D0-8.0D0*Y/3.0D0)+
1745     1DL*DLP*(64.0D0*Y*Y*Y/5.0D0-32.0D0*Y/3.0D0)+
1746     232.0D0*DL*DLP/Y/Y/15.0D0-32.0D0*DL/Y/15.0D0+
1747     3DL*(16.0D0/15.0D0+16.0D0*Y*Y/5.0D0+72.0D0*Y/5.0D0)+
1748     4DLM*(-16.0D0*Y*Y+16.0D0*Y)
1749      A2=ZETA2*(64.0D0*Y*Y*Y/5.0D0-32.0D0*Y/3.0D0)+
1750     1S11MX*(64.0D0*Y*Y*Y/5.0D0-32.0D0*Y/3.0D0)+
1751     232.0D0*S11MX/Y/Y/15.0D0+32.0D0/Y/15.0D0
1752     3+112.0D0/15.0D0+16.0D0*Y*Y/5.0D0-64.0D0*Y/5.0D0
1753      A3=DL*(16.0D0*Y)+
1754     18.0D0+8.0D0*Y-16.0D0*Y*Y
1755      FCF=CF*(A1+A2+A3*DL1)
1756      IF (NDIS.EQ.1) THEN
1757      B1=DL**2*(48.0D0*Y)+
1758     1DL*DLM*(32.0D0*Y*Y-96.0D0*Y)+
1759     2DL*DLP*(32.0D0*Y*Y+32.0D0*Y)+
1760     3DL*(16.0D0-208.0D0*Y*Y+128.0D0*Y)+
1761     4DLM**2*(-16.0D0*Y*Y+16.0D0*Y)+
1762     5DLM*(-16.0D0+464.0D0*Y*Y/3.0D0-144.0D0*Y+16.0D0/Y/3.0D0)
1763      B2=ZETA2*(32.0D0*Y*Y)+
1764     1S111MX*(-64.0D0*Y)+
1765     2S11MX*(32.0D0*Y*Y+32.0D0*Y)+16.0D0/3.0D0
1766     3-848.0D0*Y*Y/9.0D0+272.0D0*Y/3.0D0-16.0D0/Y/9.0D0
1767      B3=DL*(-64.0D0*Y)+
1768     1DLM*(32.0D0*Y-32.0D0*Y*Y)-
1769     216.0D0-80.0D0*Y+272.0D0*Y*Y/3.0D0+16.0D0/Y/3.0D0
1770      B4=-88.0D0*Y/3.0D0+88.0D0*Y*Y/3.0D0
1771      FCA=CA*(B1+B2+B3*DL1+B4*EL1)
1772      C1=DL**2*(-16.0D0*Y-16.0D0*Y*Y)+
1773     1DL*DLM*(32.0D0*Y+32.0D0*Y*Y)+
1774     2DL*(-8.0D0-144.0D0*Y-72.0D0*Y*Y)+
1775     3DLM*(8.0D0+48.0D0*Y-56.0D0*Y*Y)+
1776     4S111MX*(32.0D0*Y+32.0D0*Y*Y)-
1777     520.0D0-248.0D0*Y+268.0D0*Y*Y
1778      C2=16.0D0*Y/3.0D0-16.0D0*Y*Y/3.0D0
1779      FNF=NF*(C1+C2*EL1)
1780      ELSE
1781      B1=DL**2*(48.0D0*Y)+
1782     1DL*DLM*(32.0D0*Y*Y-96.0D0*Y)+
1783     2DL*DLP*(32.0D0*Y*Y+32.0D0*Y)+
1784     3DL*(16.0D0-208.0D0*Y*Y+128.0D0*Y)+
1785     4DLM**2*(-16.0D0*Y*Y+16.0D0*Y)+
1786     5DLM*(-16.0D0+464.0D0*Y*Y/3.0D0-144.0D0*Y+16.0D0/Y/3.0D0)
1787      B2=ZETA2*(32.0D0*Y*Y)+
1788     1S111MX*(-64.0D0*Y)+
1789     2S11MX*(32.0D0*Y*Y+32.0D0*Y)+16.0D0/3.0D0
1790     3-848.0D0*Y*Y/9.0D0+272.0D0*Y/3.0D0-16.0D0/Y/9.0D0
1791      B3=DL*(-64.0D0*Y)+
1792     1DLM*(32.0D0*Y-32.0D0*Y*Y)-
1793     216.0D0-80.0D0*Y+272.0D0*Y*Y/3.0D0+16.0D0/Y/3.0D0
1794      B4=-88.0D0*Y/3.0D0+88.0D0*Y*Y/3.0D0
1795      CA=3.0D0
1796      FCA=CA*(B1+B2+B3*DL1+B4*EL1)
1797      C1=16.0D0*Y/3.0D0-16.0D0*Y*Y/3.0D0
1798      FNF=NF*C1*EL1
1799      ENDIF
1800      DLG2=NF*(FCF+FCA+FNF)
1801      RETURN
1802      END
1803
1804
1805C THIS IS THE POLARIZED WILSON COEFFICIENT CORRESPONDING TO
1806C THE PROCESS : V + G -> Q + QB  , G1(X,Q^2) , MSBAR SCHEME .
1807
1808
1809      REAL*8 FUNCTION POLG1(Y)
1810      IMPLICIT REAL*8(A-H,O-Z)
1811      DOUBLE PRECISION NF
1812      COMMON/PAR/NF,Q2,FS,RS,DN
1813      COMMON/DIS/NDIS
1814      COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
1815      FS2=FS*FS
1816      DL1=DLOG(Q2/FS2)
1817      DL=DLOG(Y)
1818      DLM=DLOG(1.0D0-Y)
1819      A1=(8.0D0*Y-4.0D0)*(DL1+DLM-DL)+4.0D0*(3.0D0-4.0D0*Y)
1820C THE 0.5 BELOW IS DUE TO TF=1/2
1821      POLG1=0.5D0*NF*A1
1822      RETURN
1823      END
1824
1825C THIS IS THE POLARIZED WILSON COEFFICIENT CORRESPONDING TO THE PROCESS :
1826C V + Q -> Q + QB + Q (SINGLET) , G1(X,Q^2) , MSBAR SCHEME .
1827
1828
1829      REAL*8 FUNCTION POLS(Y)
1830      IMPLICIT REAL*8(A-H,O-Z)
1831      COMPLEX*16 WGPLG
1832      DOUBLE PRECISION NF
1833      COMMON/PAR/NF,Q2,FS,RS,DN
1834      COMMON/DIS/NDIS
1835      COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
1836      FS2=FS*FS
1837      DL1=DLOG(Q2/FS2)
1838      DL2=DL1*DL1
1839      S211MX=WGPLG(2,1,1.0D0-Y)
1840      S111MX=WGPLG(1,1,1.0D0-Y)
1841      S11MX=WGPLG(1,1,-Y)
1842      ZETA2=1.6449340668D0
1843      DL=DLOG(Y)
1844      DLP=DLOG(1.0D0+Y)
1845      DLM=DLOG(1.0D0-Y)
1846      A1=DL**3*(20.0D0/3.0D0+20.0D0*Y/3.0D0)+
1847     1DL**2*DLM*(-16.0D0-16.0D0*Y)+
1848     2DL**2*(50.0D0+16.0D0*Y*Y/3.0D0-10.0D0*Y)+
1849     3DL*DLM**2*(8.0D0+8.0D0*Y)+
1850     4DL*DLM*(-64.0D0+32.0D0*Y)+
1851     5DL*DLP*(-32.0D0-32.0D0*Y*Y/3.0D0-32.0D0*Y-32.0D0/Y/3.0D0)
1852      A2=DL*ZETA2*(-16.0D0-16.0D0*Y)+
1853     1DL*S111MX*(-16.0D0-16.0D0*Y)+
1854     2DL*(476.0D0/3.0D0-52.0D0*Y/3.0D0)+
1855     3DLM**2*(20.0D0-20.0D0*Y)+
1856     5DLM*S111MX*(16.0D0+16.0D0*Y)+
1857     6DLM*(-88.0D0+88.0D0*Y)
1858      A3=ZETA2*(-72.0D0-32.0D0*Y*Y/3.0D0+40.0D0*Y)+
1859     1S111MX*(-24.0D0-8.0D0*Y)+
1860     2S211MX*(-16.0D0-16.0D0*Y)+
1861     3S11MX*(-32.0D0-32.0D0*Y*Y/3.0D0-32.0D0*Y-32.0D0/Y/3.0D0)+
1862     4760.0D0/3.0D0-760.0D0*Y/3.0D0
1863      A4=DL*(8.0D0+8.0D0*Y)+
1864     120.0D0-20.0D0*Y
1865      A5=DL**2*(-16.0D0-16.0D0*Y)+
1866     1DL*DLM*(16.0D0+16.0D0*Y)+
1867     2DL*(-64.0D0+32.0D0*Y)+
1868     3DLM*(40.0D0-40.0D0*Y)+
1869     4S111MX*(16.0D0+16.0D0*Y)-
1870     588.0D0+88.0D0*Y
1871      POLS=NF*CFTF*(A1+A2+A3+A4*DL2+A5*DL1)
1872      RETURN
1873      END
1874
1875C THIS IS THE POLARIZED WILSON COEFFICIENT CORRESPONDING TO
1876C THE PROCESS : V + G -> Q + QB + G , G1(X,Q^2) , MSBAR SCHEME .
1877
1878
1879      REAL*8 FUNCTION POLG2(Y)
1880      IMPLICIT REAL*8(A-H,O-Z)
1881      COMPLEX*16 WGPLG
1882      DOUBLE PRECISION NF
1883      COMMON/PAR/NF,Q2,FS,RS,DN
1884      COMMON/DIS/NDIS
1885      COMMON/COLOUR/CA,CF,CF2,CACF,CFTF
1886      FS2=FS*FS
1887      RS2=RS*RS
1888      DL1=DLOG(Q2/FS2)
1889      DL2=DL1*DL1
1890      EL1=DLOG(FS2/RS2)
1891      S121MX=WGPLG(1,2,1.0D0-Y)
1892      S12MX=WGPLG(1,2,-Y)
1893      S211MX=WGPLG(2,1,1.0D0-Y)
1894      S21MX=WGPLG(2,1,-Y)
1895      Z=(1.0D0-Y)/(1.0D0+Y)
1896      S21Z=WGPLG(2,1,Z)
1897      S21MZ=WGPLG(2,1,-Z)
1898      S111MX=WGPLG(1,1,1.0D0-Y)
1899      S11MX=WGPLG(1,1,-Y)
1900      ZETA2=1.6449340668D0
1901      ZETA3=1.2020569031D0
1902      DL=DLOG(Y)
1903      DLP=DLOG(1.0D0+Y)
1904      DLM=DLOG(1.0D0-Y)
1905      A1=DL**3*(10.0D0/3.0D0-20.0D0*Y/3.0D0)+
1906     1DL**2*DLM*(-16.0D0+32.0D0*Y)+
1907     2DL**2*DLP*(16.0D0+16.0D0*Y*Y+32.0D0*Y)+
1908     3DL**2*(41.0D0-104.0D0*Y*Y/3.0D0-16.0D0*Y)+
1909     4DL*DLM**2*(16.0D0-32.0D0*Y)+
1910     5DL*DLM*(-88.0D0+96.0D0*Y)+
1911     6DL*DLP**2*(-32.0D0-32.0D0*Y*Y-64.0D0*Y)+
1912     7DL*DLP*(208.0D0*Y*Y/3.0D0+64.0D0*Y+64.0D0/Y/3.0D0)
1913      A2=DL*ZETA2*(-24.0D0+48.0D0*Y)+
1914     2DL*S111MX*(-8.0D0+16.0D0*Y)+
1915     3DL*S11MX*(32.0D0+32.0D0*Y*Y-64.0D0*Y)+
1916     4DL*(416.0D0/3.0D0-16.0D0*Y*Y-274.0D0*Y/3.0D0)+
1917     5DLM**3*(-20.0D0/3.0D0+40.0D0*Y/3.0D0)+
1918     6DLM**2*(54.0D0-72.0D0*Y)+
1919     7ZETA2*DLM*(-32.0D0*Y*Y)
1920      A3=DLM*S111MX*(-16.0D0+32.0D0*Y)+
1921     1DLM*(-124.0D0+16.0D0*Y*Y+104.0D0*Y)+
1922     2ZETA2*DLP*(-32.0D0-32.0D0*Y*Y-64.0D0*Y)+
1923     3DLP*S11MX*(-64.0D0-64.0D0*Y*Y-128.0D0*Y)+
1924     4ZETA2*(-108.0D0+208.0D0*Y*Y/3.0D0+144.0D0*Y)+
1925     5ZETA3*(-40.0D0+32.0D0*Y*Y+208.0D0*Y)
1926      A4=S111MX*(20.0D0-48.0D0*Y)+
1927     1S211MX*(32.0D0-64.0D0*Y)+
1928     2S121MX*(8.0D0+32.0D0*Y*Y-16.0D0*Y)+
1929     3S11MX*(208.0D0*Y*Y/3.0D0+64.0D0*Y+64.0D0/Y/3.0D0)+
1930     4S21MX*(-32.0D0-32.0D0*Y*Y+192.0D0*Y)+
1931     5S12MX*(-64.0D0-64.0D0*Y*Y-128.0D0*Y)+
1932     6466.0D0/3.0D0-478.0D0*Y/3.0D0
1933      A5=DL*(4.0D0-8.0D0*Y)+
1934     1DLM*(-8.0D0+16.0D0*Y)+
1935     26.0D0
1936      A6=DL**2*(-8.0D0+16.0D0*Y)+
1937     1DL*DLM*(24.0D0-48.0D0*Y)+
1938     2DL*(-48.0D0+32.0D0*Y)+
1939     3DLM**2*(-16.0D0+32.0D0*Y)+
1940     4DLM*(68.0D0-80.0D0*Y)+
1941     5ZETA2*(32.0D0-64.0D0*Y)+
1942     6S111MX*(-8.0D0+16.0D0*Y)-
1943     768.0D0+52.0D0*Y
1944      FCF=CF*(A1+A2+A3+A4+A5*DL2+A6*DL1)
1945      B1=DL**3*(28.0D0/3.0D0+40.0D0*Y/3.0D0)+
1946     1DL**2*DLM*(-32.0D0-16.0D0*Y)+
1947     2DL**2*DLP*(8.0D0-8.0D0*Y*Y+16.0D0*Y)+
1948     3DL**2*(90.0D0+88.0D0*Y*Y/3.0D0-32.0D0*Y)+
1949     4DL*DLM**2*(24.0D0)+
1950     5DL*DLM*DLP*(-16.0D0-32.0D0*Y)+
1951     6DL*DLM*(-136.0D0+128.0D0*Y)+
1952     8DL*DLP**2*(16.0D0+16.0D0*Y*Y+32.0D0*Y)
1953      B2=DL*DLP*(-32.0D0-176.0D0*Y*Y/3.0D0-64.0D0*Y-32.0D0/Y/3.0D0)+
1954     1DL*ZETA2*(-48.0D0)+
1955     2DL*S111MX*(-16.0D0-32.0D0*Y)+
1956     3DL*S11MX*(16.0D0-16.0D0*Y*Y+32.0D0*Y)+
1957     4DL*(944.0D0/3.0D0+8.0D0*Y*Y-208.0D0*Y/3.0D0)+
1958     5DLM**3*(-4.0D0/3.0D0+8.0D0*Y/3.0D0)+
1959     6DLM**2*(48.0D0-56.0D0*Y)+
1960     7ZETA2*DLM*(24.0D0+16.0D0*Y*Y-80.0D0*Y)
1961      B3=DLM*S111MX*(40.0D0+16.0D0*Y)+
1962     1DLM*S11MX*(-16.0D0-32.0D0*Y)+
1963     2DLM*(-176.0D0-8.0D0*Y*Y+212.0D0*Y)+
1964     3ZETA2*DLP*(16.0D0+16.0D0*Y*Y+32.0D0*Y)+
1965     4DLP*S11MX*(32.0D0+32.0D0*Y*Y+64.0D0*Y)+
1966     5S21MZ*(-16.0D0-32.0D0*Y)+
1967     6S21Z*(16.0D0+32.0D0*Y)
1968      B4=ZETA2*(-144.0D0-176.0D0*Y*Y/3.0D0+128.0D0*Y)+
1969     1ZETA3*(-12.0D0-16.0D0*Y*Y-40.0D0*Y)+
1970     2S111MX*(-24.0D0)+
1971     3S211MX*(-72.0D0-16.0D0*Y)+
1972     4S121MX*(8.0D0-16.0D0*Y*Y+16.0D0*Y)+
1973     5S11MX*(-32.0D0-176.0D0*Y*Y/3.0D0-64.0D0*Y-32.0D0/Y/3.0D0)+
1974     6S21MX*(-16.0D0+16.0D0*Y*Y-32.0D0*Y)+
1975     7S12MX*(32.0D0+32.0D0*Y*Y+64.0D0*Y)+
1976     81420.0D0/3.0D0-1468.0D0*Y/3.0D0
1977      B5=DL*(16.0D0+16.0D0*Y)+
1978     1DLM*(-8.0D0+16.0D0*Y)+
1979     248.0D0-48.0D0*Y
1980      B6=44.0D0/3.0D0-88.0D0*Y/3.0D0
1981      B7=DL**2*(-24.0D0-32.0D0*Y)+
1982     1DL*DLM*(48.0D0)+
1983     2DL*DLP*(-16.0D0-32.0D0*Y)+
1984     3DL*(-120.0D0+96.0D0*Y)+
1985     4DLM**2*(-8.0D0+16.0D0*Y)+
1986     5DLM*(112.0D0-128.0D0*Y)+
1987     6ZETA2*(-32.0D0*Y)+
1988     7S111MX*(32.0D0+32.0D0*Y)+
1989     8S11MX*(-16.0D0-32.0D0*Y)-
1990     9160.0D0+168.0D0*Y
1991      B8=44.0D0*(2.0D0*Y-1.0D0)*(DL-DLM)/3.0D0+
1992     1176.0D0*Y/3.0D0-44.0D0
1993      FCA=CA*(B1+B2+B3+B4+B5*DL2+B6*DL1*EL1+B7*DL1+B8*EL1)
1994      C1=16.0D0*Y/3.0D0-8.0D0/3.0D0
1995      C2=8.0D0*(2.0D0*Y-1.0D0)*(DLM-DL)/3.0D0-
1996     132.0D0*Y/3.0D0+8.0D0
1997      FNF=NF*(C1*DL1*EL1+C2*EL1)
1998C THE 0.5 BELOW IS DUE TO TF=1/2
1999      POLG2=0.5D0*NF*(FCF+FCA+FNF)
2000      RETURN
2001      END
2002
2003C ZEROTH ORDER ALPHAS CONTRIBUTION TO F3, F2 : V + Q -> Q + G
2004
2005
2006      DOUBLE PRECISION FUNCTION F30R(X)
2007      IMPLICIT REAL*8(A-H,O-Z)
2008      DOUBLE PRECISION NF
2009      COMMON/PAR/NF,Q2,FS,RS,DN
2010      COMMON/CHOICE/MODE
2011      F30R=1.0D0
2012      RETURN
2013      END
2014
2015
2016      DOUBLE PRECISION FUNCTION F20R(X)
2017      IMPLICIT REAL*8(A-H,O-Z)
2018      DOUBLE PRECISION NF
2019      COMMON/PAR/NF,Q2,FS,RS,DN
2020      COMMON/CHOICE/MODE
2021      F20R=1.0D0
2022      RETURN
2023      END
2024
2025
2026C ORDER ALPHAS CONTRIBUTION TO F3, F2 : V + Q -> Q + G
2027C REGULAR PART , MSBAR SCHEME .
2028
2029
2030      DOUBLE PRECISION FUNCTION F31R(X)
2031      IMPLICIT REAL*8(A-H,O-Z)
2032      DOUBLE PRECISION NF
2033      COMMON/PAR/NF,Q2,FS,RS,DN
2034      COMMON/CHOICE/MODE
2035      FA=X**(DN-1.d0)
2036      F31R=FA*C3NS1(X)
2037      RETURN
2038      END
2039
2040      DOUBLE PRECISION FUNCTION F21R(X)
2041      IMPLICIT REAL*8(A-H,O-Z)
2042      DOUBLE PRECISION NF
2043      COMMON/PAR/NF,Q2,FS,RS,DN
2044      COMMON/CHOICE/MODE
2045      FA=X**DN
2046      F21R=FA*C2NS1(X)
2047      RETURN
2048      END
2049
2050      DOUBLE PRECISION FUNCTION FG21R(X)
2051      IMPLICIT REAL*8(A-H,O-Z)
2052      DOUBLE PRECISION NF
2053      COMMON/PAR/NF,Q2,FS,RS,DN
2054      COMMON/CHOICE/MODE
2055      FA=X**DN
2056csm gluon coeff function
2057      FG21R=FA*C2G1(X)
2058      RETURN
2059      END
2060
2061      DOUBLE PRECISION FUNCTION FG22R(X)
2062      IMPLICIT REAL*8(A-H,O-Z)
2063      DOUBLE PRECISION NF
2064      COMMON/PAR/NF,Q2,FS,RS,DN
2065      COMMON/CHOICE/MODE
2066      FA=X**DN
2067csm gluon coeff function
2068      FG22R=FA*C2G2(X)
2069      RETURN
2070      END
2071
2072
2073      DOUBLE PRECISION FUNCTION FPS22R(X)
2074      IMPLICIT REAL*8(A-H,O-Z)
2075      DOUBLE PRECISION NF
2076      COMMON/PAR/NF,Q2,FS,RS,DN
2077      COMMON/CHOICE/MODE
2078      FA=X**DN
2079csm pure singlet coeff function
2080      FPS22R=FA*C2S(X)
2081      RETURN
2082      END
2083
2084      DOUBLE PRECISION FUNCTION FL1R(X)
2085      IMPLICIT REAL*8(A-H,O-Z)
2086      DOUBLE PRECISION NF
2087      COMMON/PAR/NF,Q2,FS,RS,DN
2088      COMMON/CHOICE/MODE
2089      FA=X**DN
2090      FL1R=FA*CLNS1(X)
2091      RETURN
2092      END
2093
2094      DOUBLE PRECISION FUNCTION FL2R(X)
2095      IMPLICIT REAL*8(A-H,O-Z)
2096      DOUBLE PRECISION NF
2097      COMMON/PAR/NF,Q2,FS,RS,DN
2098      COMMON/CHOICE/MODE
2099      FA=X**DN
2100      FL2R=FA*(CLNSP(X)+CLNSM(X))
2101      RETURN
2102      END
2103
2104      DOUBLE PRECISION FUNCTION FPSL2R(X)
2105      IMPLICIT REAL*8(A-H,O-Z)
2106      DOUBLE PRECISION NF
2107      COMMON/PAR/NF,Q2,FS,RS,DN
2108      COMMON/CHOICE/MODE
2109      FA=X**DN
2110csm pure singlet coeff function
2111      FPSL2R=FA*CLS(X)
2112      RETURN
2113      END
2114
2115
2116      DOUBLE PRECISION FUNCTION FLG1R(X)
2117      IMPLICIT REAL*8(A-H,O-Z)
2118      DOUBLE PRECISION NF
2119      COMMON/PAR/NF,Q2,FS,RS,DN
2120      COMMON/CHOICE/MODE
2121      FA=X**DN
2122      FLG1R=FA*CLG1(X)
2123      RETURN
2124      END
2125
2126
2127      DOUBLE PRECISION FUNCTION FLG2R(X)
2128      IMPLICIT REAL*8(A-H,O-Z)
2129      DOUBLE PRECISION NF
2130      COMMON/PAR/NF,Q2,FS,RS,DN
2131      COMMON/CHOICE/MODE
2132      FA=X**DN
2133      FLG2R=FA*CLG2(X)
2134      RETURN
2135      END
2136
2137
2138C ORDER ALPHAS CONTRIBUTION TO F3, F2 : V + Q -> Q + G
2139C SINGULAR PART , MSBAR SCHEME .
2140
2141      DOUBLE PRECISION FUNCTION F31S(X)
2142      IMPLICIT REAL*8(A-H,O-Z)
2143      DOUBLE PRECISION NF
2144      COMMON/PAR/NF,Q2,FS,RS,DN
2145      COMMON/CHOICE/MODE
2146      FA=X**(DN-1.d0)
2147      F31S=(FA-1.0D0)*SOFT1(X)
2148      RETURN
2149      END
2150
2151      DOUBLE PRECISION FUNCTION F21S(X)
2152      IMPLICIT REAL*8(A-H,O-Z)
2153      DOUBLE PRECISION NF
2154      COMMON/PAR/NF,Q2,FS,RS,DN
2155      COMMON/CHOICE/MODE
2156      FA=X**DN
2157      F21S=(FA-1.0D0)*SOFT1(X)
2158      RETURN
2159      END
2160
2161
2162
2163C ORDER ALPHAS**2 CORRECTION TO F3, F2 : V + Q -> Q + G + G ;
2164C V + Q -> Q + Q + QB (ID + NONID)
2165C REGULAR PART , MSBAR SCHEME .
2166
2167
2168      DOUBLE PRECISION FUNCTION F32R(X)
2169      IMPLICIT REAL*8(A-H,O-Z)
2170      DOUBLE PRECISION NF
2171      COMMON/PAR/NF,Q2,FS,RS,DN
2172      COMMON/CHOICE/MODE
2173      FA=X**(DN-1.d0)
2174      F32R=FA*(C3NSP(X)-C3NSM(X))
2175      RETURN
2176      END
2177
2178      DOUBLE PRECISION FUNCTION F22R(X)
2179      IMPLICIT REAL*8(A-H,O-Z)
2180      DOUBLE PRECISION NF
2181      COMMON/PAR/NF,Q2,FS,RS,DN
2182      COMMON/CHOICE/MODE
2183      FA=X**DN
2184C THIS IS THE COMBINATION FOR F2 IN CHARGED CURRENT INTERACTIONS
2185C      F22R=FA*(C2NSP(X)-C2NSM(X))
2186C THIS IS THE COMBINATION FOR F2 IN NEUTRAL CURRENT INTERACTIONS
2187      F22R=FA*(C2NSP(X)+C2NSM(X))
2188      RETURN
2189      END
2190
2191
2192C ORDER ALPHAS**2 CORRECTION TO F3, F2 : V + Q -> Q + G
2193C SINGULAR PART , MSBAR SCHEME .
2194
2195      DOUBLE PRECISION FUNCTION F32S(X)
2196      IMPLICIT REAL*8(A-H,O-Z)
2197      DOUBLE PRECISION NF
2198      COMMON/PAR/NF,Q2,FS,RS,DN
2199      COMMON/CHOICE/MODE
2200      FA=X**(DN-1.d0)
2201      F32S=(FA-1.0D0)*SOFT2(X)
2202      RETURN
2203      END
2204
2205
2206      DOUBLE PRECISION FUNCTION F22S(X)
2207      IMPLICIT REAL*8(A-H,O-Z)
2208      DOUBLE PRECISION NF
2209      COMMON/PAR/NF,Q2,FS,RS,DN
2210      COMMON/CHOICE/MODE
2211      FA=X**DN
2212      F22S=(FA-1.0D0)*SOFT2(X)
2213      RETURN
2214      END
2215
2216
2217C ORDER ALPHAS CONTRIBUTION TO G1 : V + G -> Q + QB
2218C MSBAR SCHEME .
2219
2220
2221      DOUBLE PRECISION FUNCTION F1G1(X)
2222      IMPLICIT REAL*8(A-H,O-Z)
2223      DOUBLE PRECISION NF
2224      COMMON/PAR/NF,Q2,FS,RS,DN
2225      COMMON/CHOICE/MODE
2226      FA=X**DN
2227      F1G1=FA*POLG1(X)
2228      RETURN
2229      END
2230
2231
2232C ORDER ALPHAS**2 CORRECTION TO G1, : V + G -> G + Q + QB ;
2233C MSBAR SCHEME .
2234
2235
2236      DOUBLE PRECISION FUNCTION F1G2(X)
2237      IMPLICIT REAL*8(A-H,O-Z)
2238      DOUBLE PRECISION NF
2239      COMMON/PAR/NF,Q2,FS,RS,DN
2240      COMMON/CHOICE/MODE
2241      FA=X**DN
2242      F1G2=FA*POLG2(X)
2243      RETURN
2244      END
2245
2246
2247C ORDER ALPHAS**2 CORRECTION TO G1, : V + Q -> Q + Q + QB ;
2248C MSBAR SCHEME .
2249
2250      DOUBLE PRECISION FUNCTION F1Q2(X)
2251      IMPLICIT REAL*8(A-H,O-Z)
2252      DOUBLE PRECISION NF
2253      COMMON/PAR/NF,Q2,FS,RS,DN
2254      COMMON/CHOICE/MODE
2255      FA=X**DN
2256      F1Q2=FA*POLS(X)
2257      RETURN
2258      END
2259
2260
2261
2262
2263C HERE FOLLOWS THE TWO LOOP CORRECTED RUNNING COUPLING CONSTANT
2264
2265      DOUBLE PRECISION FUNCTION ALFAS1(RS,LAMBDA)
2266      IMPLICIT REAL*8(A-H,O-Z)
2267      DOUBLE PRECISION LAMBDA,NF
2268      COMMON/PAR/NF,Q2,FSS,RSS,DN
2269      ALFAS1=ALF(RS,LAMBDA,NF)
2270      RETURN
2271      END
2272
2273
2274      DOUBLE PRECISION FUNCTION ALFAS2(RS,LAMBDA)
2275      IMPLICIT REAL*8(A-H,O-Z)
2276      DOUBLE PRECISION LAMBDA,MB,MC
2277      MB=4.75D0
2278      MC=1.5D0
2279      A=1.0D0
2280      AMB=A*MB
2281      AMC=A*MC
2282      IF (RS.LT.AMB) GO TO 1
2283      ALF1=ALF(RS,LAMBDA,5.0D0)
2284      GO TO 3
2285    1 IF (RS.LT.AMC) GO TO 2
2286      ALF1I=1.0D0/ALF(RS,LAMBDA,4.0D0)
2287      ALF2I=1.0D0/ALF(AMB,LAMBDA,5.0D0)
2288      ALF3I=1.0D0/ALF(AMB,LAMBDA,4.0D0)
2289      ALF1=1.0D0/(ALF1I+ALF2I-ALF3I)
2290      GO TO 3
2291    2 ALF1I=1.0D0/ALF(RS,LAMBDA,3.0D0)
2292      ALF2I=1.0D0/ALF(AMC,LAMBDA,4.0D0)
2293      ALF3I=1.0D0/ALF(AMB,LAMBDA,5.0D0)
2294      ALF4I=1.0D0/ALF(AMB,LAMBDA,4.0D0)
2295      ALF5I=1.0D0/ALF(AMC,LAMBDA,3.0D0)
2296      ALF1=1.0D0/(ALF1I+ALF2I+ALF3I-ALF4I-ALF5I)
2297    3 ALFAS2=ALF1
2298      RETURN
2299      END
2300
2301
2302      DOUBLE PRECISION FUNCTION ALF(RS,LAMBDA,NF)
2303      IMPLICIT REAL*8(A-H,O-Z)
2304      DOUBLE PRECISION LAMBDA,LAMBD2,NF
2305      PI=3.14159265359D0
2306      LAMBD2=LAMBDA*LAMBDA
2307      BF=(33.0D0-2.0D0*NF)/12.0D0/PI
2308      BFP=(153.0D0-19.0D0*NF)/(66.0D0-4.0D0*NF)/PI
2309      RS2=RS*RS
2310      DL1=DLOG(RS2/LAMBD2)
2311      DL2=DLOG(DL1)
2312      ALF=(1.0D0-(BFP*DL2)/(BF*DL1))/DL1/BF
2313C      ALF=1.0D0/DL1/BF
2314      RETURN
2315      END
2316      SUBROUTINE STRUCT(Z,FS,MODE,UV,DV,SS,ST,CH,BOT,GL)
2317      IMPLICIT REAL*8 (A-H,O-Z)
2318      UV=8*Z*(1-Z)**3
2319      DV=5*Z*(1-Z)**4
2320      GL=21*Z*(1-Z)**5
2321      SS=9*Z*(1-Z)**8
2322      ST=SS
2323      CH=0.0D0
2324      BOT=0.0D0
2325      RETURN
2326      END
2327      INCLUDE 'daind.f'
2328      INCLUDE 'daind1.f'
2329      INCLUDE 'wgplg.f'