Analysis Software
Documentation for sPHENIX simulation software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
pytbhg.f
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file pytbhg.f
1 
2 C AMPLITUDE SQUARED (MATRIX ELEMENTS) FOR THE PROCESSES:
3 C GG->TBH^+, QQBAR->TBH^+
4 C AS A FUNCTION OF 4-MOMENTA FOR SUITABLE INTERFACE
5 C (FOR INSTANCE WITH PYTHIA)
6 C------------------------------------------------------------
7 C BASED ON F. BORZUMATI, J.-L. KNEUR, N. POLONSKY HEP-PH/9905443,
8 C PHYS REV. D 60 (1999) 115011
9 C (THESE FILES PREPARED BY J.-L. KNEUR)
10 C------------------------------------------------------------
11 C 1) GG->TBH^+
12  SUBROUTINE pytbhg(Q1,Q2,P1,P2,P3,MT,MB,RMB,MHP,AMP2)
13 C
14 C CONVENTIONS AND INPUT/OUTPUT DEFINITIONS:
15 C
16 C INPUT: Q1,Q2 ARE ENTERING 4-MOMENTA OF INITIAL GLUONS OR QUARKS;
17 C P1, P2 ARE THE TOP AND BOTTOM OUTGOING 4-MOMENTA;
18 C P3 IS OUTGOING CHARGED HIGGS 4-MOMENTA.
19 C (NB FOR ALL 4-MOMENTA P(4) IS TIME-COMPONENT)
20 C "PHYSICAL PARAMETERS" INPUT:
21 C MT,MB TOP AND BOTTOM MASSES;
22 C MHP CHARGED HIGGS MASS
23 C FURTHER PARAMETERS INPUT IS NEEDED FROM COMMON/PARAM/ (SEE BELOW)
24 C
25 C OUTPUT: AMP2 IS MATRIX ELEMENT (AMPLITUDE**2) FOR GG->TB H^+
26 C (NB AMP2 IS TRULY AMPLITUDE SQUARRED, I.E. WITHOUT ANY
27 C PHASE SPACE FACTORS INCLUDED. IT INCLUDES COLOUR AND COUPLING
28 C FACTORS, AS EXPLICIT BELOW. ACCORDINGLY, FOR EXAMPLE THE TOTAL
29 C CROSS-SECTION SHOULD BE (SYMBOLICALLY):
30 C SIGMA = INTEGRATE [PARTON DENSITY FUNCTIONS * 3-PARTICLE FINAL
31 C STATE PHASE-SPACE (STANDARDLY NORMALIZED) * AMP2 ]
32 C
33  IMPLICIT DOUBLE PRECISION(a-h, o-z)
34  IMPLICIT INTEGER(i-n)
35  DOUBLE PRECISION mw2,mt,mb,mhp,mw
36  dimension q1(4),q2(4),p1(4),p2(4),p3(4)
37  common/pydat1/mstu(200),paru(200),mstj(200),parj(200)
38  common/pydat2/kchg(500,4),pmas(500,4),parf(2000),vckm(4,4)
39  common/pymssm/imss(0:99),rmss(0:99)
40 
41  common/pyctbh/ alpha,alphas,sw2,mw2,tanb,vtb,v,a
42  SAVE /pydat1/,/pydat2/,/pymssm/,/pyctbh/
43 C !THE RELEVANT INPUT PARAMETERS ABOVE ARE NEEDED FOR CALCULATION
44 C BUT ARE NOT DEFINED HERE SO THAT ONE MAY CHOOSE/VARY THEIR VALUES:
45 C ACCORDINGLY, WHEN CALLING THESE SUBROUTINES, PLEASE SUPPLY VIA
46 C THIS COMMON/PARAM/ YOUR PREFERRED ALPHA, ALPHAS,..AND TANB
47 C (TAN BETA) VALUES
48 C
49 C THE NORMALIZED V,A COUPLINGS ARE DEFINED BELOW AND USED BOTH
50 C IN THIS ROUTINE AND IN THE TOP WIDTH CALCULATION PYTBHB(..).
51 
52  pi = 4*datan(1.d0)
53  mw = dsqrt(mw2)
54 C
55 C COLLECTING THE RELEVANT OVERALL FACTORS:
56 C 8X8 INITIAL GLUON COLOR AVERAGE, 2X2 GLUON SPIN AVERAGE
57  ps=1.d0/(8.d0*8.d0 *2.d0*2.d0)
58 C COUPLING CONSTANT (OVERALL NORMALIZATION)
59  fact=(4.d0*pi*alpha)*(4.d0*pi*alphas)**2/sw2/2.d0
60 C NB ALPHA IS E^2/4/PI, BUT BETTER DEFINED IN TERMS OF G_FERMI:
61 C ALPHA= DSQRT(2.D0)*GF*SW2*MW**2/PI
62 C ALPHAS IS ALPHA_STRONG;
63 C SW2 IS SIN(THETA_W)**2.
64 C
65 C VTB=.998D0
66 C VTB IS TOP-BOTTOM CKM MATRIX ELEMENT (APPROXIMATE VALUE HERE)
67 C
68  v = ( mt/mw/tanb +rmb/mw*tanb)/2.d0
69  a = (-mt/mw/tanb +rmb/mw*tanb)/2.d0
70 C V AND A ARE (NORMALIZED) VECTOR AND AXIAL TBH^+ COUPLINGS
71 C
72 C REDEFINING P2 INGOING FROM OVERALL MOMENTUM CONSERVATION
73 C (BECAUSE P2 INGOING WAS USED IN OUR GRAPH CALCULATION CONVENTIONS)
74  DO 100 kk=1,4
75  p2(kk)=p3(kk)-q1(kk)-q2(kk)+p1(kk)
76  100 CONTINUE
77 C DEFINING VARIOUS RELEVANT 4-SCALAR PRODUCTS:
78  s = 2*pytbhs(q1,q2)
79  p1q1=pytbhs(q1,p1)
80  p1q2=pytbhs(p1,q2)
81  p2q1=pytbhs(p2,q1)
82  p2q2=pytbhs(p2,q2)
83  p1p2=pytbhs(p1,p2)
84 C
85 C TOP WIDTH CALCULATION
86  CALL pytbhb(mt,mb,mhp,br,gamt)
87 C GAMT IS THE TOP WIDTH: T->BH^+ AND/OR T->B W^+
88 C THEN DEFINE TOP (RESONANT) PROPAGATOR:
89  a1inv= s -2*p1q1 -2*p1q2
90  a1 =a1inv/(a1inv**2+ (gamt*mt)**2)
91 C (I.E. INTRODUCE THE TOP WIDTH IN A1 TO REGULARISE THE POLE)
92 C NB: A12 = A1*A1 BUT CORRECT EXPRESSION BELOW BECAUSE OF
93 C THE TOP WIDTH
94  a12 = 1.d0/(a1inv**2+ (gamt*mt)**2)
95  a2 =1.d0/(s +2*p2q1 +2*p2q2)
96 C NOTE A2 IS B PROPAGATOR, DOES NOT NEED A WIDTH
97 C NOW COMES THE AMP**2:
98 C NB COLOR FACTOR (COMING FROM GRAPHS) ALREADY INCLUDED IN
99 C THE EXPRESSIONS BELOW
100  v18=0.d0
101  a18=0.d0
102  v18= 640*a1/3+640*a2/3+32*a1*a2*mb**2-368*a12*mb*mt-
103  &512*a1*a2*mb*mt/3-
104  &368*a2**2*mb*mt+32*a1*a2*mt**2+496*a12*p1p2/3+
105  &320*a1*a2*p1p2+496*a2**2*p1p2/3+128*a1*mb*mt**3/(3*p1q1**2)+
106  &128*a1*mt**4/(3*p1q1**2)-256*a12*mb*mt**5/(3*p1q1**2)+
107  &256*a1*mt**2*p1p2/(3*p1q1**2)-256*a12*mt**4*p1p2/(3*p1q1**2)+
108  &8/(3*p1q1)-32*a1*mb*mt/p1q1-56*a2*mb*mt/(3*p1q1)+
109  &88*a1*mt**2/(3*p1q1)+72*a2*mt**2/p1q1+
110  &704*a12*mb*mt**3/(3*p1q1)-224*a1*a2*mb*mt**3/(3*p1q1)+
111  &104*a1*p1p2/(3*p1q1)+48*a2*p1p2/p1q1+
112  &128*a1*a2*mb*mt*p1p2/(3*p1q1)+512*a12*mt**2*p1p2/(3*p1q1)-
113  &448*a1*a2*mt**2*p1p2/(3*p1q1)-32*a1*a2*p1p2**2/p1q1-
114  &656*a1*a2*p1q1/3-224*a2**2*p1q1+128*a1*mb*mt**3/(3*p1q2**2)+
115  &128*a1*mt**4/(3*p1q2**2)-256*a12*mb*mt**5/(3*p1q2**2)+
116  &256*a1*mt**2*p1p2/(3*p1q2**2)-256*a12*mt**4*p1p2/(3*p1q2**2)+
117  &256*a1*mt**2*p1q1/(3*p1q2**2)+256*a12*mb*mt**3*p1q1/(3*p1q2**2)+
118  &8/(3*p1q2)-32*a1*mb*mt/p1q2-56*a2*mb*mt/(3*p1q2)
119  v18=v18+88*a1*mt**2/(3*p1q2)+72*a2*mt**2/p1q2+
120  &704*a12*mb*mt**3/(3*p1q2)-224*a1*a2*mb*mt**3/(3*p1q2)+
121  &104*a1*p1p2/(3*p1q2)+48*a2*p1p2/p1q2+
122  &128*a1*a2*mb*mt*p1p2/(3*p1q2)+512*a12*mt**2*p1p2/(3*p1q2)-
123  &448*a1*a2*mt**2*p1p2/(3*p1q2)-32*a1*a2*p1p2**2/p1q2-
124  &32*a1*mb*mt**3/(3*p1q1*p1q2)-32*a1*mt**4/(3*p1q1*p1q2)+
125  &64*a12*mb*mt**5/(3*p1q1*p1q2)+16*p1p2/(3*p1q1*p1q2)-
126  &64*a1*mt**2*p1p2/(3*p1q1*p1q2)+64*a12*mt**4*p1p2/(3*p1q1*p1q2)+
127  &112*a1*p1q1/p1q2+272*a2*p1q1/(3*p1q2)-
128  &272*a1*a2*mb**2*p1q1/(3*p1q2)+208*a12*mb*mt*p1q1/(3*p1q2)-
129  &400*a1*a2*mb*mt*p1q1/(3*p1q2)-80*a1*a2*mt**2*p1q1/p1q2+
130  &96*a12*p1p2*p1q1/p1q2-320*a1*a2*p1p2*p1q1/p1q2-
131  &544*a1*a2*p1q1**2/(3*p1q2)-656*a1*a2*p1q2/3-224*a2**2*p1q2+
132  &256*a1*mt**2*p1q2/(3*p1q1**2)+256*a12*mb*mt**3*p1q2/(3*p1q1**2)+
133  &112*a1*p1q2/p1q1+272*a2*p1q2/(3*p1q1)-
134  &272*a1*a2*mb**2*p1q2/(3*p1q1)+208*a12*mb*mt*p1q2/(3*p1q1)-
135  &400*a1*a2*mb*mt*p1q2/(3*p1q1)-80*a1*a2*mt**2*p1q2/p1q1
136  v18=v18+96*a12*p1p2*p1q2/p1q1-320*a1*a2*p1p2*p1q2/p1q1-
137  &544*a1*a2*p1q2**2/(3*p1q1)+128*a2*mb**4/(3*p2q1**2)+
138  &128*a2*mb**3*mt/(3*p2q1**2)-256*a2**2*mb**5*mt/(3*p2q1**2)+
139  &256*a2*mb**2*p1p2/(3*p2q1**2)-256*a2**2*mb**4*p1p2/(3*p2q1**2)+
140  &256*a2*mb**2*p1q1/(3*p2q1**2)-256*a2**2*mb**4*p1q1/(3*p2q1**2)-
141  &64*mb**3*mt**3/(3*p1q2**2*p2q1**2)-
142  &64*mb**2*mt**2*p1p2/(3*p1q2**2*p2q1**2)-
143  &64*mb**2*mt**2*p1q1/(3*p1q2**2*p2q1**2)+
144  &64*mb**3*mt/(3*p1q2*p2q1**2)+
145  &256*a2*mb**3*mt*p1p2/(3*p1q2*p2q1**2)+
146  &256*a2*mb**2*p1p2**2/(3*p1q2*p2q1**2)+
147  &256*a2*mb**3*mt*p1q1/(3*p1q2*p2q1**2)+
148  &512*a2*mb**2*p1p2*p1q1/(3*p1q2*p2q1**2)+
149  &256*a2*mb**2*p1q1**2/(3*p1q2*p2q1**2)-
150  &256*a2**2*mb**4*p1q2/(3*p2q1**2)-8/(3*p2q1)-72*a1*mb**2/p2q1-
151  &88*a2*mb**2/(3*p2q1)+56*a1*mb*mt/(3*p2q1)+32*a2*mb*mt/p2q1+
152  &224*a1*a2*mb**3*mt/(3*p2q1)-704*a2**2*mb**3*mt/(3*p2q1)
153  v18=v18-48*a1*p1p2/p2q1-104*a2*p1p2/(3*p2q1)+
154  &448*a1*a2*mb**2*p1p2/(3*p2q1)-512*a2**2*mb**2*p1p2/(3*p2q1)-
155  &128*a1*a2*mb*mt*p1p2/(3*p2q1)+32*a1*a2*p1p2**2/p2q1-
156  &16*p1p2/(3*p1q1*p2q1)-32*a1*mb*mt*p1p2/(3*p1q1*p2q1)-
157  &32*a2*mb*mt*p1p2/(3*p1q1*p2q1)-
158  &64*a1*a2*mb*mt*p1p2**2/(3*p1q1*p2q1)-
159  &64*a1*a2*p1p2**3/(3*p1q1*p2q1)-256*a2*p1q1/(3*p2q1)+
160  &448*a1*a2*mb**2*p1q1/(3*p2q1)-368*a2**2*mb**2*p1q1/(3*p2q1)+
161  &224*a1*a2*mb*mt*p1q1/(3*p2q1)+304*a1*a2*p1p2*p1q1/(3*p2q1)-
162  &64*mb*mt**3/(3*p1q2**2*p2q1)-
163  &256*a1*mb*mt**3*p1p2/(3*p1q2**2*p2q1)-
164  &256*a1*mt**2*p1p2**2/(3*p1q2**2*p2q1)+
165  &64*mt**2*p1q1/(3*p1q2**2*p2q1)-
166  &128*a1*mb**2*mt**2*p1q1/(3*p1q2**2*p2q1)-
167  &128*a1*mb*mt**3*p1q1/(3*p1q2**2*p2q1)-
168  &256*a1*mt**2*p1p2*p1q1/(3*p1q2**2*p2q1)-4*mb**2/(3*p1q2*p2q1)+
169  &64*mb*mt/(3*p1q2*p2q1)-128*a2*mb**3*mt/(3*p1q2*p2q1)
170  v18=v18-4*mt**2/(3*p1q2*p2q1)-128*a1*mb**2*mt**2/(3*p1q2*p2q1)-
171  &128*a2*mb**2*mt**2/(3*p1q2*p2q1)-128*a1*mb*mt**3/(3*p1q2*p2q1)-
172  &112*a2*mb**2*p1p2/(3*p1q2*p2q1)-32*a1*mb*mt*p1p2/(3*p1q2*p2q1)-
173  &32*a2*mb*mt*p1p2/(3*p1q2*p2q1)-112*a1*mt**2*p1p2/(3*p1q2*p2q1)-
174  &48*a1*p1p2**2/(p1q2*p2q1)-48*a2*p1p2**2/(p1q2*p2q1)+
175  &512*a1*a2*mb*mt*p1p2**2/(3*p1q2*p2q1)+
176  &512*a1*a2*p1p2**3/(3*p1q2*p2q1)-8*mb*mt*p1p2/(3*p1q1*p1q2*p2q1)-
177  &8*mt**2*p1p2/(3*p1q1*p1q2*p2q1)+
178  &32*a1*mb*mt**3*p1p2/(3*p1q1*p1q2*p2q1)-
179  &16*p1p2**2/(3*p1q1*p1q2*p2q1)+
180  &32*a1*mt**2*p1p2**2/(3*p1q1*p1q2*p2q1)+8*p1q1/(3*p1q2*p2q1)-
181  &160*a1*mb**2*p1q1/(3*p1q2*p2q1)-272*a2*mb**2*p1q1/(3*p1q2*p2q1)+
182  &56*a1*mb*mt*p1q1/(3*p1q2*p2q1)+200*a2*mb*mt*p1q1/(3*p1q2*p2q1)-
183  &48*a1*p1p2*p1q1/(p1q2*p2q1)-256*a2*p1p2*p1q1/(3*p1q2*p2q1)+
184  &256*a1*a2*mb**2*p1p2*p1q1/(3*p1q2*p2q1)+
185  &256*a1*a2*mb*mt*p1p2*p1q1/(p1q2*p2q1)+
186  &1024*a1*a2*p1p2**2*p1q1/(3*p1q2*p2q1)
187  v18=v18-272*a2*p1q1**2/(3*p1q2*p2q1)+
188  &256*a1*a2*mb**2*p1q1**2/(3*p1q2*p2q1)+
189  &256*a1*a2*mb*mt*p1q1**2/(3*p1q2*p2q1)+
190  &512*a1*a2*p1p2*p1q1**2/(3*p1q2*p2q1)+16*a2*p1q2/(3*p2q1)+
191  &64*a1*a2*mb**2*p1q2/p2q1+32*a2**2*mb**2*p1q2/(3*p2q1)+
192  &112*a1*a2*mb*mt*p1q2/(3*p2q1)+368*a1*a2*p1p2*p1q2/(3*p2q1)+
193  &32*a2*p1p2*p1q2/(3*p1q1*p2q1)-
194  &32*a1*a2*mb**2*p1p2*p1q2/(3*p1q1*p2q1)-
195  &32*a1*a2*mb*mt*p1p2*p1q2/(3*p1q1*p2q1)-
196  &64*a1*a2*p1p2**2*p1q2/(3*p1q1*p2q1)+224*a12*p2q1+
197  &656*a1*a2*p2q1/3-256*a1*mt**2*p2q1/(3*p1q1**2)+
198  &256*a12*mt**4*p2q1/(3*p1q1**2)-256*a1*p2q1/(3*p1q1)+
199  &224*a1*a2*mb*mt*p2q1/(3*p1q1)-368*a12*mt**2*p2q1/(3*p1q1)+
200  &448*a1*a2*mt**2*p2q1/(3*p1q1)+304*a1*a2*p1p2*p2q1/(3*p1q1)+
201  &256*a12*mt**4*p2q1/(3*p1q2**2)+
202  &256*a12*mt**2*p1q1*p2q1/(3*p1q2**2)+16*a1*p2q1/(3*p1q2)+
203  &112*a1*a2*mb*mt*p2q1/(3*p1q2)+32*a12*mt**2*p2q1/(3*p1q2)
204  v18=v18+64*a1*a2*mt**2*p2q1/p1q2+368*a1*a2*p1p2*p2q1/(3*p1q2)+
205  &16*a1*mt**2*p2q1/(3*p1q1*p1q2)-64*a12*mt**4*p2q1/(3*p1q1*p1q2)+
206  &640*a12*p1q1*p2q1/(3*p1q2)+544*a1*a2*p1q1*p2q1/(3*p1q2)+
207  &32*a12*p1q2*p2q1/p1q1+944*a1*a2*p1q2*p2q1/(3*p1q1)+
208  &128*a2*mb**4/(3*p2q2**2)+128*a2*mb**3*mt/(3*p2q2**2)-
209  &256*a2**2*mb**5*mt/(3*p2q2**2)+256*a2*mb**2*p1p2/(3*p2q2**2)-
210  &256*a2**2*mb**4*p1p2/(3*p2q2**2)-
211  &64*mb**3*mt**3/(3*p1q1**2*p2q2**2)-
212  &64*mb**2*mt**2*p1p2/(3*p1q1**2*p2q2**2)+
213  &64*mb**3*mt/(3*p1q1*p2q2**2)+
214  &256*a2*mb**3*mt*p1p2/(3*p1q1*p2q2**2)+
215  &256*a2*mb**2*p1p2**2/(3*p1q1*p2q2**2)-
216  &256*a2**2*mb**4*p1q1/(3*p2q2**2)+256*a2*mb**2*p1q2/(3*p2q2**2)-
217  &256*a2**2*mb**4*p1q2/(3*p2q2**2)-
218  &64*mb**2*mt**2*p1q2/(3*p1q1**2*p2q2**2)+
219  &256*a2*mb**3*mt*p1q2/(3*p1q1*p2q2**2)+
220  &512*a2*mb**2*p1p2*p1q2/(3*p1q1*p2q2**2)
221  v18=v18+256*a2*mb**2*p1q2**2/(3*p1q1*p2q2**2)-
222  &256*a2*mb**2*p2q1/(3*p2q2**2)-256*a2**2*mb**3*mt*p2q1/(3*p2q2**2)+
223  &64*mb**2*mt**2*p2q1/(3*p1q1**2*p2q2**2)+
224  &64*mb**2*p2q1/(3*p1q1*p2q2**2)-
225  &128*a2*mb**3*mt*p2q1/(3*p1q1*p2q2**2)-
226  &128*a2*mb**2*mt**2*p2q1/(3*p1q1*p2q2**2)-
227  &256*a2*mb**2*p1p2*p2q1/(3*p1q1*p2q2**2)+
228  &256*a2**2*mb**2*p1q1*p2q1/(3*p2q2**2)-
229  &256*a2*mb**2*p1q2*p2q1/(3*p1q1*p2q2**2)-8/(3*p2q2)-
230  &72*a1*mb**2/p2q2-88*a2*mb**2/(3*p2q2)+56*a1*mb*mt/(3*p2q2)+
231  &32*a2*mb*mt/p2q2+224*a1*a2*mb**3*mt/(3*p2q2)-
232  &704*a2**2*mb**3*mt/(3*p2q2)-48*a1*p1p2/p2q2-
233  &104*a2*p1p2/(3*p2q2)+448*a1*a2*mb**2*p1p2/(3*p2q2)-
234  &512*a2**2*mb**2*p1p2/(3*p2q2)-128*a1*a2*mb*mt*p1p2/(3*p2q2)+
235  &32*a1*a2*p1p2**2/p2q2-64*mb*mt**3/(3*p1q1**2*p2q2)-
236  &256*a1*mb*mt**3*p1p2/(3*p1q1**2*p2q2)-
237  &256*a1*mt**2*p1p2**2/(3*p1q1**2*p2q2)-4*mb**2/(3*p1q1*p2q2)
238  v18=v18+64*mb*mt/(3*p1q1*p2q2)-128*a2*mb**3*mt/(3*p1q1*p2q2)-
239  &4*mt**2/(3*p1q1*p2q2)-128*a1*mb**2*mt**2/(3*p1q1*p2q2)-
240  &128*a2*mb**2*mt**2/(3*p1q1*p2q2)-128*a1*mb*mt**3/(3*p1q1*p2q2)-
241  &112*a2*mb**2*p1p2/(3*p1q1*p2q2)-32*a1*mb*mt*p1p2/(3*p1q1*p2q2)-
242  &32*a2*mb*mt*p1p2/(3*p1q1*p2q2)-112*a1*mt**2*p1p2/(3*p1q1*p2q2)-
243  &48*a1*p1p2**2/(p1q1*p2q2)-48*a2*p1p2**2/(p1q1*p2q2)+
244  &512*a1*a2*mb*mt*p1p2**2/(3*p1q1*p2q2)+
245  &512*a1*a2*p1p2**3/(3*p1q1*p2q2)+16*a2*p1q1/(3*p2q2)+
246  &64*a1*a2*mb**2*p1q1/p2q2+32*a2**2*mb**2*p1q1/(3*p2q2)+
247  &112*a1*a2*mb*mt*p1q1/(3*p2q2)+368*a1*a2*p1p2*p1q1/(3*p2q2)-
248  &16*p1p2/(3*p1q2*p2q2)-32*a1*mb*mt*p1p2/(3*p1q2*p2q2)-
249  &32*a2*mb*mt*p1p2/(3*p1q2*p2q2)-
250  &64*a1*a2*mb*mt*p1p2**2/(3*p1q2*p2q2)-
251  &64*a1*a2*p1p2**3/(3*p1q2*p2q2)-8*mb*mt*p1p2/(3*p1q1*p1q2*p2q2)-
252  &8*mt**2*p1p2/(3*p1q1*p1q2*p2q2)+
253  &32*a1*mb*mt**3*p1p2/(3*p1q1*p1q2*p2q2)-
254  &16*p1p2**2/(3*p1q1*p1q2*p2q2)
255  v18=v18+32*a1*mt**2*p1p2**2/(3*p1q1*p1q2*p2q2)+
256  &32*a2*p1p2*p1q1/(3*p1q2*p2q2)-
257  &32*a1*a2*mb**2*p1p2*p1q1/(3*p1q2*p2q2)-
258  &32*a1*a2*mb*mt*p1p2*p1q1/(3*p1q2*p2q2)-
259  &64*a1*a2*p1p2**2*p1q1/(3*p1q2*p2q2)-256*a2*p1q2/(3*p2q2)+
260  &448*a1*a2*mb**2*p1q2/(3*p2q2)-368*a2**2*mb**2*p1q2/(3*p2q2)+
261  &224*a1*a2*mb*mt*p1q2/(3*p2q2)+304*a1*a2*p1p2*p1q2/(3*p2q2)+
262  &64*mt**2*p1q2/(3*p1q1**2*p2q2)-
263  &128*a1*mb**2*mt**2*p1q2/(3*p1q1**2*p2q2)-
264  &128*a1*mb*mt**3*p1q2/(3*p1q1**2*p2q2)-
265  &256*a1*mt**2*p1p2*p1q2/(3*p1q1**2*p2q2)+8*p1q2/(3*p1q1*p2q2)-
266  &160*a1*mb**2*p1q2/(3*p1q1*p2q2)-272*a2*mb**2*p1q2/(3*p1q1*p2q2)+
267  &56*a1*mb*mt*p1q2/(3*p1q1*p2q2)+200*a2*mb*mt*p1q2/(3*p1q1*p2q2)-
268  &48*a1*p1p2*p1q2/(p1q1*p2q2)-256*a2*p1p2*p1q2/(3*p1q1*p2q2)+
269  &256*a1*a2*mb**2*p1p2*p1q2/(3*p1q1*p2q2)+
270  &256*a1*a2*mb*mt*p1p2*p1q2/(p1q1*p2q2)+
271  &1024*a1*a2*p1p2**2*p1q2/(3*p1q1*p2q2)
272  v18=v18-272*a2*p1q2**2/(3*p1q1*p2q2)+
273  &256*a1*a2*mb**2*p1q2**2/(3*p1q1*p2q2)+
274  &256*a1*a2*mb*mt*p1q2**2/(3*p1q1*p2q2)+
275  &512*a1*a2*p1p2*p1q2**2/(3*p1q1*p2q2)-32*a2*mb**4/(3*p2q1*p2q2)-
276  &32*a2*mb**3*mt/(3*p2q1*p2q2)+64*a2**2*mb**5*mt/(3*p2q1*p2q2)+
277  &16*p1p2/(3*p2q1*p2q2)-64*a2*mb**2*p1p2/(3*p2q1*p2q2)+
278  &64*a2**2*mb**4*p1p2/(3*p2q1*p2q2)+8*mb**2*p1p2/(3*p1q1*p2q1*p2q2)+
279  &8*mb*mt*p1p2/(3*p1q1*p2q1*p2q2)-
280  &32*a2*mb**3*mt*p1p2/(3*p1q1*p2q1*p2q2)+
281  &16*p1p2**2/(3*p1q1*p2q1*p2q2)-
282  &32*a2*mb**2*p1p2**2/(3*p1q1*p2q1*p2q2)-
283  &16*a2*mb**2*p1q1/(3*p2q1*p2q2)+64*a2**2*mb**4*p1q1/(3*p2q1*p2q2)+
284  &8*mb**2*p1p2/(3*p1q2*p2q1*p2q2)+8*mb*mt*p1p2/(3*p1q2*p2q1*p2q2)-
285  &32*a2*mb**3*mt*p1p2/(3*p1q2*p2q1*p2q2)+
286  &16*p1p2**2/(3*p1q2*p2q1*p2q2)-
287  &32*a2*mb**2*p1p2**2/(3*p1q2*p2q1*p2q2)+
288  &16*mb*mt*p1p2**2/(3*p1q1*p1q2*p2q1*p2q2)
289  v18=v18+16*p1p2**3/(3*p1q1*p1q2*p2q1*p2q2)-
290  &32*a2*mb**2*p1p2*p1q1/(3*p1q2*p2q1*p2q2)-
291  &16*a2*mb**2*p1q2/(3*p2q1*p2q2)+64*a2**2*mb**4*p1q2/(3*p2q1*p2q2)-
292  &32*a2*mb**2*p1p2*p1q2/(3*p1q1*p2q1*p2q2)+272*a1*p2q1/(3*p2q2)+
293  &112*a2*p2q1/p2q2-80*a1*a2*mb**2*p2q1/p2q2-
294  &400*a1*a2*mb*mt*p2q1/(3*p2q2)+208*a2**2*mb*mt*p2q1/(3*p2q2)-
295  &272*a1*a2*mt**2*p2q1/(3*p2q2)-320*a1*a2*p1p2*p2q1/p2q2+
296  &96*a2**2*p1p2*p2q1/p2q2+256*a1*mb*mt**3*p2q1/(3*p1q1**2*p2q2)+
297  &512*a1*mt**2*p1p2*p2q1/(3*p1q1**2*p2q2)-8*p2q1/(3*p1q1*p2q2)-
298  &200*a1*mb*mt*p2q1/(3*p1q1*p2q2)-56*a2*mb*mt*p2q1/(3*p1q1*p2q2)+
299  &272*a1*mt**2*p2q1/(3*p1q1*p2q2)+160*a2*mt**2*p2q1/(3*p1q1*p2q2)+
300  &256*a1*p1p2*p2q1/(3*p1q1*p2q2)+48*a2*p1p2*p2q1/(p1q1*p2q2)-
301  &256*a1*a2*mb*mt*p1p2*p2q1/(p1q1*p2q2)-
302  &256*a1*a2*mt**2*p1p2*p2q1/(3*p1q1*p2q2)-
303  &1024*a1*a2*p1p2**2*p2q1/(3*p1q1*p2q2)-
304  &544*a1*a2*p1q1*p2q1/(3*p2q2)-640*a2**2*p1q1*p2q1/(3*p2q2)-
305  &32*a1*p1p2*p2q1/(3*p1q2*p2q2)
306  v18=v18+32*a1*a2*mb*mt*p1p2*p2q1/(3*p1q2*p2q2)+
307  &32*a1*a2*mt**2*p1p2*p2q1/(3*p1q2*p2q2)+
308  &64*a1*a2*p1p2**2*p2q1/(3*p1q2*p2q2)-
309  &32*a1*mt**2*p1p2*p2q1/(3*p1q1*p1q2*p2q2)+
310  &64*a1*a2*p1p2*p1q1*p2q1/(3*p1q2*p2q2)-
311  &944*a1*a2*p1q2*p2q1/(3*p2q2)-32*a2**2*p1q2*p2q1/p2q2+
312  &256*a1*mt**2*p1q2*p2q1/(3*p1q1**2*p2q2)+
313  &96*a1*p1q2*p2q1/(p1q1*p2q2)+96*a2*p1q2*p2q1/(p1q1*p2q2)-
314  &128*a1*a2*mb**2*p1q2*p2q1/(3*p1q1*p2q2)-
315  &256*a1*a2*mb*mt*p1q2*p2q1/(p1q1*p2q2)-
316  &128*a1*a2*mt**2*p1q2*p2q1/(3*p1q1*p2q2)-
317  &512*a1*a2*p1p2*p1q2*p2q1/(p1q1*p2q2)-
318  &512*a1*a2*p1q2**2*p2q1/(3*p1q1*p2q2)+544*a1*a2*p2q1**2/(3*p2q2)-
319  &256*a1*mt**2*p2q1**2/(3*p1q1**2*p2q2)-
320  &272*a1*p2q1**2/(3*p1q1*p2q2)+
321  &256*a1*a2*mb*mt*p2q1**2/(3*p1q1*p2q2)+
322  &256*a1*a2*mt**2*p2q1**2/(3*p1q1*p2q2)
323  v18=v18+512*a1*a2*p1p2*p2q1**2/(3*p1q1*p2q2)+
324  &512*a1*a2*p1q2*p2q1**2/(3*p1q1*p2q2)+224*a12*p2q2+
325  &656*a1*a2*p2q2/3+256*a12*mt**4*p2q2/(3*p1q1**2)+
326  &16*a1*p2q2/(3*p1q1)+112*a1*a2*mb*mt*p2q2/(3*p1q1)+
327  &32*a12*mt**2*p2q2/(3*p1q1)+64*a1*a2*mt**2*p2q2/p1q1+
328  &368*a1*a2*p1p2*p2q2/(3*p1q1)-256*a1*mt**2*p2q2/(3*p1q2**2)+
329  &256*a12*mt**4*p2q2/(3*p1q2**2)-256*a1*p2q2/(3*p1q2)+
330  &224*a1*a2*mb*mt*p2q2/(3*p1q2)-368*a12*mt**2*p2q2/(3*p1q2)+
331  &448*a1*a2*mt**2*p2q2/(3*p1q2)+304*a1*a2*p1p2*p2q2/(3*p1q2)+
332  &16*a1*mt**2*p2q2/(3*p1q1*p1q2)-64*a12*mt**4*p2q2/(3*p1q1*p1q2)+
333  &32*a12*p1q1*p2q2/p1q2+944*a1*a2*p1q1*p2q2/(3*p1q2)+
334  &256*a12*mt**2*p1q2*p2q2/(3*p1q1**2)+
335  &640*a12*p1q2*p2q2/(3*p1q1)+544*a1*a2*p1q2*p2q2/(3*p1q1)-
336  &256*a2*mb**2*p2q2/(3*p2q1**2)-256*a2**2*mb**3*mt*p2q2/(3*p2q1**2)+
337  &64*mb**2*mt**2*p2q2/(3*p1q2**2*p2q1**2)+
338  &64*mb**2*p2q2/(3*p1q2*p2q1**2)-
339  &128*a2*mb**3*mt*p2q2/(3*p1q2*p2q1**2)
340  v18=v18-128*a2*mb**2*mt**2*p2q2/(3*p1q2*p2q1**2)-
341  &256*a2*mb**2*p1p2*p2q2/(3*p1q2*p2q1**2)-
342  &256*a2*mb**2*p1q1*p2q2/(3*p1q2*p2q1**2)+
343  &256*a2**2*mb**2*p1q2*p2q2/(3*p2q1**2)+272*a1*p2q2/(3*p2q1)+
344  &112*a2*p2q2/p2q1-80*a1*a2*mb**2*p2q2/p2q1-
345  &400*a1*a2*mb*mt*p2q2/(3*p2q1)+208*a2**2*mb*mt*p2q2/(3*p2q1)-
346  &272*a1*a2*mt**2*p2q2/(3*p2q1)-320*a1*a2*p1p2*p2q2/p2q1+
347  &96*a2**2*p1p2*p2q2/p2q1-32*a1*p1p2*p2q2/(3*p1q1*p2q1)+
348  &32*a1*a2*mb*mt*p1p2*p2q2/(3*p1q1*p2q1)+
349  &32*a1*a2*mt**2*p1p2*p2q2/(3*p1q1*p2q1)+
350  &64*a1*a2*p1p2**2*p2q2/(3*p1q1*p2q1)-944*a1*a2*p1q1*p2q2/(3*p2q1)-
351  &32*a2**2*p1q1*p2q2/p2q1+256*a1*mb*mt**3*p2q2/(3*p1q2**2*p2q1)+
352  &512*a1*mt**2*p1p2*p2q2/(3*p1q2**2*p2q1)+
353  &256*a1*mt**2*p1q1*p2q2/(3*p1q2**2*p2q1)-8*p2q2/(3*p1q2*p2q1)-
354  &200*a1*mb*mt*p2q2/(3*p1q2*p2q1)-56*a2*mb*mt*p2q2/(3*p1q2*p2q1)+
355  &272*a1*mt**2*p2q2/(3*p1q2*p2q1)+160*a2*mt**2*p2q2/(3*p1q2*p2q1)+
356  &256*a1*p1p2*p2q2/(3*p1q2*p2q1)+48*a2*p1p2*p2q2/(p1q2*p2q1)
357  v18=v18-256*a1*a2*mb*mt*p1p2*p2q2/(p1q2*p2q1)-
358  &256*a1*a2*mt**2*p1p2*p2q2/(3*p1q2*p2q1)-
359  &1024*a1*a2*p1p2**2*p2q2/(3*p1q2*p2q1)-
360  &32*a1*mt**2*p1p2*p2q2/(3*p1q1*p1q2*p2q1)+
361  &96*a1*p1q1*p2q2/(p1q2*p2q1)+96*a2*p1q1*p2q2/(p1q2*p2q1)-
362  &128*a1*a2*mb**2*p1q1*p2q2/(3*p1q2*p2q1)-
363  &256*a1*a2*mb*mt*p1q1*p2q2/(p1q2*p2q1)-
364  &128*a1*a2*mt**2*p1q1*p2q2/(3*p1q2*p2q1)-
365  &512*a1*a2*p1p2*p1q1*p2q2/(p1q2*p2q1)-
366  &512*a1*a2*p1q1**2*p2q2/(3*p1q2*p2q1)-544*a1*a2*p1q2*p2q2/(3*p2q1)-
367  &640*a2**2*p1q2*p2q2/(3*p2q1)+
368  &64*a1*a2*p1p2*p1q2*p2q2/(3*p1q1*p2q1)+544*a1*a2*p2q2**2/(3*p2q1)-
369  &256*a1*mt**2*p2q2**2/(3*p1q2**2*p2q1)-
370  &272*a1*p2q2**2/(3*p1q2*p2q1)+
371  &256*a1*a2*mb*mt*p2q2**2/(3*p1q2*p2q1)+
372  &256*a1*a2*mt**2*p2q2**2/(3*p1q2*p2q1)+
373  &512*a1*a2*p1p2*p2q2**2/(3*p1q2*p2q1)
374  v18=v18+512*a1*a2*p1q1*p2q2**2/(3*p1q2*p2q1)+
375  &384*a12*mb*mt*p1q1**2/s**2+
376  &384*a12*p1p2*p1q1**2/s**2+2688*a12*mb*mt*p1q1*p1q2/s**2+
377  &2688*a12*p1p2*p1q1*p1q2/s**2+384*a12*mb*mt*p1q2**2/s**2+
378  &384*a12*p1p2*p1q2**2/s**2+768*a1*a2*mb*mt*p1q1*p2q1/s**2+
379  &768*a1*a2*p1p2*p1q1*p2q1/s**2+2688*a1*a2*mb*mt*p1q2*p2q1/s**2+
380  &2688*a1*a2*p1p2*p1q2*p2q1/s**2-960*a12*p1q1*p1q2*p2q1/s**2-
381  &960*a1*a2*p1q1*p1q2*p2q1/s**2+960*a12*p1q2**2*p2q1/s**2+
382  &960*a1*a2*p1q2**2*p2q1/s**2+384*a2**2*mb*mt*p2q1**2/s**2+
383  &384*a2**2*p1p2*p2q1**2/s**2-960*a1*a2*p1q2*p2q1**2/s**2-
384  &960*a2**2*p1q2*p2q1**2/s**2+2688*a1*a2*mb*mt*p1q1*p2q2/s**2+
385  &2688*a1*a2*p1p2*p1q1*p2q2/s**2+960*a12*p1q1**2*p2q2/s**2+
386  &960*a1*a2*p1q1**2*p2q2/s**2+768*a1*a2*mb*mt*p1q2*p2q2/s**2+
387  &768*a1*a2*p1p2*p1q2*p2q2/s**2-960*a12*p1q1*p1q2*p2q2/s**2-
388  &960*a1*a2*p1q1*p1q2*p2q2/s**2+2688*a2**2*mb*mt*p2q1*p2q2/s**2+
389  &2688*a2**2*p1p2*p2q1*p2q2/s**2+960*a1*a2*p1q1*p2q1*p2q2/s**2+
390  &960*a2**2*p1q1*p2q1*p2q2/s**2+960*a1*a2*p1q2*p2q1*p2q2/s**2+
391  &960*a2**2*p1q2*p2q1*p2q2/s**2+384*a2**2*mb*mt*p2q2**2/s**2
392  v18=v18+384*a2**2*p1p2*p2q2**2/s**2-960*a1*a2*p1q1*p2q2**2/s**2-
393  &960*a2**2*p1q1*p2q2**2/s**2+96*a1*mb*mt/s+96*a2*mb*mt/s-
394  &768*a2**2*mb**3*mt/s-768*a12*mb*mt**3/s-192*a1*p1p2/s-
395  &192*a2*p1p2/s-768*a2**2*mb**2*p1p2/s-2304*a1*a2*mb*mt*p1p2/s-
396  &768*a12*mt**2*p1p2/s-2304*a1*a2*p1p2**2/s-
397  &96*a1*mb*mt**3/(p1q1*s)-192*a2*mb*mt*p1p2/(p1q1*s)-
398  &96*a1*mt**2*p1p2/(p1q1*s)-192*a2*p1p2**2/(p1q1*s)-192*a1*p1q1/s-
399  &144*a2*p1q1/s-384*a1*a2*mb**2*p1q1/s-480*a2**2*mb**2*p1q1/s-
400  &480*a12*mb*mt*p1q1/s+96*a1*a2*mb*mt*p1q1/s-
401  &864*a12*p1p2*p1q1/s-672*a1*a2*p1p2*p1q1/s-96*a1*a2*p1q1**2/s-
402  &96*a1*mb*mt**3/(p1q2*s)-192*a2*mb*mt*p1p2/(p1q2*s)-
403  &96*a1*mt**2*p1p2/(p1q2*s)-192*a2*p1p2**2/(p1q2*s)-
404  &48*a1*mb*mt*p1q1/(p1q2*s)+96*a2*mb*mt*p1q1/(p1q2*s)-
405  &48*a1*mt**2*p1q1/(p1q2*s)-192*a1*p1p2*p1q1/(p1q2*s)-
406  &192*a2*p1p2*p1q1/(p1q2*s)+192*a1*a2*mb*mt*p1p2*p1q1/(p1q2*s)+
407  &192*a1*a2*p1p2**2*p1q1/(p1q2*s)-192*a1*p1q1**2/(p1q2*s)-
408  &192*a2*p1q1**2/(p1q2*s)+192*a1*a2*mb**2*p1q1**2/(p1q2*s)
409  v18=v18-192*a12*mb*mt*p1q1**2/(p1q2*s)+
410  &96*a1*a2*mb*mt*p1q1**2/(p1q2*s)+
411  &192*a1*a2*p1p2*p1q1**2/(p1q2*s)-192*a1*p1q2/s-144*a2*p1q2/s-
412  &384*a1*a2*mb**2*p1q2/s-480*a2**2*mb**2*p1q2/s-
413  &480*a12*mb*mt*p1q2/s+96*a1*a2*mb*mt*p1q2/s-
414  &864*a12*p1p2*p1q2/s-672*a1*a2*p1p2*p1q2/s-
415  &48*a1*mb*mt*p1q2/(p1q1*s)+96*a2*mb*mt*p1q2/(p1q1*s)-
416  &48*a1*mt**2*p1q2/(p1q1*s)-192*a1*p1p2*p1q2/(p1q1*s)-
417  &192*a2*p1p2*p1q2/(p1q1*s)+192*a1*a2*mb*mt*p1p2*p1q2/(p1q1*s)+
418  &192*a1*a2*p1p2**2*p1q2/(p1q1*s)-576*a1*a2*p1q1*p1q2/s-
419  &96*a1*a2*p1q2**2/s-192*a1*p1q2**2/(p1q1*s)-
420  &192*a2*p1q2**2/(p1q1*s)+192*a1*a2*mb**2*p1q2**2/(p1q1*s)-
421  &192*a12*mb*mt*p1q2**2/(p1q1*s)+96*a1*a2*mb*mt*p1q2**2/(p1q1*s)+
422  &192*a1*a2*p1p2*p1q2**2/(p1q1*s)+96*a2*mb**3*mt/(p2q1*s)+
423  &96*a2*mb**2*p1p2/(p2q1*s)+192*a1*mb*mt*p1p2/(p2q1*s)+
424  &192*a1*p1p2**2/(p2q1*s)+96*a1*mb**2*p1q1/(p2q1*s)+
425  &192*a2*mb**2*p1q1/(p2q1*s)+96*a1*mb*mt*p1q1/(p2q1*s)+
426  &192*a1*a2*mb**3*mt*p1q1/(p2q1*s)+192*a1*p1p2*p1q1/(p2q1*s)
427  v18=v18+192*a1*a2*mb**2*p1p2*p1q1/(p2q1*s)+
428  &96*a1*a2*mb**2*p1q1**2/(p2q1*s)+
429  &192*a2*mb**3*mt*p1q1/(p1q2*p2q1*s)+
430  &192*a2*mb**2*p1p2*p1q1/(p1q2*p2q1*s)+
431  &96*a1*mb*mt*p1p2*p1q1/(p1q2*p2q1*s)+
432  &96*a1*p1p2**2*p1q1/(p1q2*p2q1*s)+
433  &96*a1*mb**2*p1q1**2/(p1q2*p2q1*s)+
434  &192*a2*mb**2*p1q1**2/(p1q2*p2q1*s)+
435  &48*a1*mb*mt*p1q1**2/(p1q2*p2q1*s)+
436  &96*a1*p1p2*p1q1**2/(p1q2*p2q1*s)+96*a1*mb**2*p1q2/(p2q1*s)+
437  &48*a2*mb**2*p1q2/(p2q1*s)-192*a1*a2*mb**3*mt*p1q2/(p2q1*s)-
438  &192*a1*a2*mb**2*p1p2*p1q2/(p2q1*s)-
439  &96*a1*a2*mb**2*p1q2**2/(p2q1*s)+144*a1*p2q1/s+192*a2*p2q1/s-
440  &96*a1*a2*mb*mt*p2q1/s+480*a2**2*mb*mt*p2q1/s+
441  &480*a12*mt**2*p2q1/s+384*a1*a2*mt**2*p2q1/s+
442  &672*a1*a2*p1p2*p2q1/s+864*a2**2*p1p2*p2q1/s+
443  &96*a2*mb*mt*p2q1/(p1q1*s)+192*a1*mt**2*p2q1/(p1q1*s)
444  v18=v18+96*a2*mt**2*p2q1/(p1q1*s)+
445  &192*a1*a2*mb*mt**3*p2q1/(p1q1*s)+
446  &192*a2*p1p2*p2q1/(p1q1*s)+192*a1*a2*mt**2*p1p2*p2q1/(p1q1*s)-
447  &192*a12*p1q1*p2q1/s-192*a2**2*p1q1*p2q1/s+
448  &48*a1*mt**2*p2q1/(p1q2*s)+96*a2*mt**2*p2q1/(p1q2*s)-
449  &192*a1*a2*mb*mt**3*p2q1/(p1q2*s)-
450  &192*a1*a2*mt**2*p1p2*p2q1/(p1q2*s)-
451  &96*a1*a2*mb*mt*p1q1*p2q1/(p1q2*s)-
452  &192*a12*mt**2*p1q1*p2q1/(p1q2*s)-
453  &96*a1*a2*mt**2*p1q1*p2q1/(p1q2*s)-
454  &384*a1*a2*p1p2*p1q1*p2q1/(p1q2*s)-384*a12*p1q1**2*p2q1/(p1q2*s)-
455  &384*a1*a2*p1q1**2*p2q1/(p1q2*s)-480*a12*p1q2*p2q1/s-
456  &960*a1*a2*p1q2*p2q1/s-480*a2**2*p1q2*p2q1/s+
457  &144*a1*p1q2*p2q1/(p1q1*s)+96*a2*p1q2*p2q1/(p1q1*s)-
458  &384*a1*a2*mb*mt*p1q2*p2q1/(p1q1*s)-
459  &96*a12*mt**2*p1q2*p2q1/(p1q1*s)+
460  &96*a1*a2*mt**2*p1q2*p2q1/(p1q1*s)-
461  &576*a1*a2*p1p2*p1q2*p2q1/(p1q1*s)-192*a12*p1q2**2*p2q1/(p1q1*s)
462  v18=v18-384*a1*a2*p1q2**2*p2q1/(p1q1*s)-96*a1*a2*p2q1**2/s-
463  &96*a1*a2*mt**2*p2q1**2/(p1q1*s)+96*a1*a2*mt**2*p2q1**2/(p1q2*s)+
464  &288*a1*a2*p1q2*p2q1**2/(p1q1*s)+96*a2*mb**3*mt/(p2q2*s)+
465  &96*a2*mb**2*p1p2/(p2q2*s)+192*a1*mb*mt*p1p2/(p2q2*s)+
466  &192*a1*p1p2**2/(p2q2*s)+96*a1*mb**2*p1q1/(p2q2*s)+
467  &48*a2*mb**2*p1q1/(p2q2*s)-192*a1*a2*mb**3*mt*p1q1/(p2q2*s)-
468  &192*a1*a2*mb**2*p1p2*p1q1/(p2q2*s)-
469  &96*a1*a2*mb**2*p1q1**2/(p2q2*s)+96*a1*mb**2*p1q2/(p2q2*s)+
470  &192*a2*mb**2*p1q2/(p2q2*s)+96*a1*mb*mt*p1q2/(p2q2*s)+
471  &192*a1*a2*mb**3*mt*p1q2/(p2q2*s)+192*a1*p1p2*p1q2/(p2q2*s)+
472  &192*a1*a2*mb**2*p1p2*p1q2/(p2q2*s)+
473  &192*a2*mb**3*mt*p1q2/(p1q1*p2q2*s)+
474  &192*a2*mb**2*p1p2*p1q2/(p1q1*p2q2*s)+
475  &96*a1*mb*mt*p1p2*p1q2/(p1q1*p2q2*s)+
476  &96*a1*p1p2**2*p1q2/(p1q1*p2q2*s)+96*a1*a2*mb**2*p1q2**2/(p2q2*s)+
477  &96*a1*mb**2*p1q2**2/(p1q1*p2q2*s)+
478  &192*a2*mb**2*p1q2**2/(p1q1*p2q2*s)
479  v18=v18+48*a1*mb*mt*p1q2**2/(p1q1*p2q2*s)+
480  &96*a1*p1p2*p1q2**2/(p1q1*p2q2*s)-48*a2*mb**2*p2q1/(p2q2*s)+
481  &96*a1*mb*mt*p2q1/(p2q2*s)-48*a2*mb*mt*p2q1/(p2q2*s)-
482  &192*a1*p1p2*p2q1/(p2q2*s)-192*a2*p1p2*p2q1/(p2q2*s)+
483  &192*a1*a2*mb*mt*p1p2*p2q1/(p2q2*s)+
484  &192*a1*a2*p1p2**2*p2q1/(p2q2*s)-
485  &192*a1*mb*mt**3*p2q1/(p1q1*p2q2*s)-
486  &96*a2*mb*mt*p1p2*p2q1/(p1q1*p2q2*s)-
487  &192*a1*mt**2*p1p2*p2q1/(p1q1*p2q2*s)-
488  &96*a2*p1p2**2*p2q1/(p1q1*p2q2*s)+
489  &96*a1*a2*mb**2*p1q1*p2q1/(p2q2*s)+
490  &192*a2**2*mb**2*p1q1*p2q1/(p2q2*s)+
491  &96*a1*a2*mb*mt*p1q1*p2q1/(p2q2*s)+
492  &384*a1*a2*p1p2*p1q1*p2q1/(p2q2*s)-96*a1*p1q2*p2q1/(p2q2*s)-
493  &144*a2*p1q2*p2q1/(p2q2*s)-96*a1*a2*mb**2*p1q2*p2q1/(p2q2*s)+
494  &96*a2**2*mb**2*p1q2*p2q1/(p2q2*s)+
495  &384*a1*a2*mb*mt*p1q2*p2q1/(p2q2*s)
496  v18=v18+576*a1*a2*p1p2*p1q2*p2q1/(p2q2*s)-
497  &96*a2*mb**2*p1q2*p2q1/(p1q1*p2q2*s)+
498  &48*a1*mb*mt*p1q2*p2q1/(p1q1*p2q2*s)+
499  &48*a2*mb*mt*p1q2*p2q1/(p1q1*p2q2*s)-
500  &96*a1*mt**2*p1q2*p2q1/(p1q1*p2q2*s)-
501  &96*a1*p1p2*p1q2*p2q1/(p1q1*p2q2*s)-
502  &96*a2*p1p2*p1q2*p2q1/(p1q1*p2q2*s)+
503  &96*a1*a2*p1q1*p1q2*p2q1/(p2q2*s)+288*a1*a2*p1q2**2*p2q1/(p2q2*s)-
504  &96*a1*p1q2**2*p2q1/(p1q1*p2q2*s)-96*a2*p1q2**2*p2q1/(p1q1*p2q2*s)+
505  &192*a1*p2q1**2/(p2q2*s)+192*a2*p2q1**2/(p2q2*s)-
506  &96*a1*a2*mb*mt*p2q1**2/(p2q2*s)+192*a2**2*mb*mt*p2q1**2/(p2q2*s)-
507  &192*a1*a2*mt**2*p2q1**2/(p2q2*s)-192*a1*a2*p1p2*p2q1**2/(p2q2*s)+
508  &48*a2*mb*mt*p2q1**2/(p1q1*p2q2*s)+
509  &192*a1*mt**2*p2q1**2/(p1q1*p2q2*s)+
510  &96*a2*mt**2*p2q1**2/(p1q1*p2q2*s)+
511  &96*a2*p1p2*p2q1**2/(p1q1*p2q2*s)-384*a1*a2*p1q1*p2q1**2/(p2q2*s)-
512  &384*a2**2*p1q1*p2q1**2/(p2q2*s)-384*a1*a2*p1q2*p2q1**2/(p2q2*s)
513  v18=v18-192*a2**2*p1q2*p2q1**2/(p2q2*s)+
514  &96*a1*p1q2*p2q1**2/(p1q1*p2q2*s)+
515  &96*a2*p1q2*p2q1**2/(p1q1*p2q2*s)+144*a1*p2q2/s+192*a2*p2q2/s-
516  &96*a1*a2*mb*mt*p2q2/s+480*a2**2*mb*mt*p2q2/s+
517  &480*a12*mt**2*p2q2/s+384*a1*a2*mt**2*p2q2/s+
518  &672*a1*a2*p1p2*p2q2/s+864*a2**2*p1p2*p2q2/s+
519  &48*a1*mt**2*p2q2/(p1q1*s)+96*a2*mt**2*p2q2/(p1q1*s)-
520  &192*a1*a2*mb*mt**3*p2q2/(p1q1*s)-
521  &192*a1*a2*mt**2*p1p2*p2q2/(p1q1*s)-480*a12*p1q1*p2q2/s-
522  &960*a1*a2*p1q1*p2q2/s-480*a2**2*p1q1*p2q2/s+
523  &96*a2*mb*mt*p2q2/(p1q2*s)+192*a1*mt**2*p2q2/(p1q2*s)+
524  &96*a2*mt**2*p2q2/(p1q2*s)+192*a1*a2*mb*mt**3*p2q2/(p1q2*s)+
525  &192*a2*p1p2*p2q2/(p1q2*s)+192*a1*a2*mt**2*p1p2*p2q2/(p1q2*s)+
526  &144*a1*p1q1*p2q2/(p1q2*s)+96*a2*p1q1*p2q2/(p1q2*s)-
527  &384*a1*a2*mb*mt*p1q1*p2q2/(p1q2*s)-
528  &96*a12*mt**2*p1q1*p2q2/(p1q2*s)+
529  &96*a1*a2*mt**2*p1q1*p2q2/(p1q2*s)
530  v18=v18-576*a1*a2*p1p2*p1q1*p2q2/(p1q2*s)-
531  &192*a12*p1q1**2*p2q2/(p1q2*s)-
532  &384*a1*a2*p1q1**2*p2q2/(p1q2*s)-192*a12*p1q2*p2q2/s-
533  &192*a2**2*p1q2*p2q2/s-96*a1*a2*mb*mt*p1q2*p2q2/(p1q1*s)-
534  &192*a12*mt**2*p1q2*p2q2/(p1q1*s)-
535  &96*a1*a2*mt**2*p1q2*p2q2/(p1q1*s)-
536  &384*a1*a2*p1p2*p1q2*p2q2/(p1q1*s)-384*a12*p1q2**2*p2q2/(p1q1*s)-
537  &384*a1*a2*p1q2**2*p2q2/(p1q1*s)-48*a2*mb**2*p2q2/(p2q1*s)+
538  &96*a1*mb*mt*p2q2/(p2q1*s)-48*a2*mb*mt*p2q2/(p2q1*s)-
539  &192*a1*p1p2*p2q2/(p2q1*s)-192*a2*p1p2*p2q2/(p2q1*s)+
540  &192*a1*a2*mb*mt*p1p2*p2q2/(p2q1*s)+
541  &192*a1*a2*p1p2**2*p2q2/(p2q1*s)-96*a1*p1q1*p2q2/(p2q1*s)-
542  &144*a2*p1q1*p2q2/(p2q1*s)-96*a1*a2*mb**2*p1q1*p2q2/(p2q1*s)+
543  &96*a2**2*mb**2*p1q1*p2q2/(p2q1*s)+
544  &384*a1*a2*mb*mt*p1q1*p2q2/(p2q1*s)+
545  &576*a1*a2*p1p2*p1q1*p2q2/(p2q1*s)+288*a1*a2*p1q1**2*p2q2/(p2q1*s)-
546  &192*a1*mb*mt**3*p2q2/(p1q2*p2q1*s)
547  v18=v18-96*a2*mb*mt*p1p2*p2q2/(p1q2*p2q1*s)-
548  &192*a1*mt**2*p1p2*p2q2/(p1q2*p2q1*s)-
549  &96*a2*p1p2**2*p2q2/(p1q2*p2q1*s)-
550  &96*a2*mb**2*p1q1*p2q2/(p1q2*p2q1*s)+
551  &48*a1*mb*mt*p1q1*p2q2/(p1q2*p2q1*s)
552 
553  v18bis=
554  &48*a2*mb*mt*p1q1*p2q2/(p1q2*p2q1*s)-
555  &96*a1*mt**2*p1q1*p2q2/(p1q2*p2q1*s)-
556  &96*a1*p1p2*p1q1*p2q2/(p1q2*p2q1*s)-
557  &96*a2*p1p2*p1q1*p2q2/(p1q2*p2q1*s)-
558  &96*a1*p1q1**2*p2q2/(p1q2*p2q1*s)-96*a2*p1q1**2*p2q2/(p1q2*p2q1*s)+
559  &96*a1*a2*mb**2*p1q2*p2q2/(p2q1*s)+
560  &192*a2**2*mb**2*p1q2*p2q2/(p2q1*s)+
561  &96*a1*a2*mb*mt*p1q2*p2q2/(p2q1*s)+
562  &384*a1*a2*p1p2*p1q2*p2q2/(p2q1*s)+
563  &96*a1*a2*p1q1*p1q2*p2q2/(p2q1*s)-576*a1*a2*p2q1*p2q2/s+
564  &96*a1*a2*p1q1*p2q1*p2q2/(p1q2*s)+96*a1*a2*p1q2*p2q1*p2q2/(p1q1*s)-
565  &96*a1*a2*p2q2**2/s+96*a1*a2*mt**2*p2q2**2/(p1q1*s)-
566  &96*a1*a2*mt**2*p2q2**2/(p1q2*s)+288*a1*a2*p1q1*p2q2**2/(p1q2*s)+
567  &192*a1*p2q2**2/(p2q1*s)+192*a2*p2q2**2/(p2q1*s)-
568  &96*a1*a2*mb*mt*p2q2**2/(p2q1*s)+192*a2**2*mb*mt*p2q2**2/(p2q1*s)-
569  &192*a1*a2*mt**2*p2q2**2/(p2q1*s)-192*a1*a2*p1p2*p2q2**2/(p2q1*s)
570  v18bis=v18bis-384*a1*a2*p1q1*p2q2**2/(p2q1*s)-
571  &192*a2**2*p1q1*p2q2**2/(p2q1*s)+
572  &48*a2*mb*mt*p2q2**2/(p1q2*p2q1*s)+
573  &192*a1*mt**2*p2q2**2/(p1q2*p2q1*s)+
574  &96*a2*mt**2*p2q2**2/(p1q2*p2q1*s)+
575  &96*a2*p1p2*p2q2**2/(p1q2*p2q1*s)+96*a1*p1q1*p2q2**2/(p1q2*p2q1*s)+
576  &96*a2*p1q1*p2q2**2/(p1q2*p2q1*s)-384*a1*a2*p1q2*p2q2**2/(p2q1*s)-
577  &384*a2**2*p1q2*p2q2**2/(p2q1*s)+512*a1*a2*s/3-
578  &128*a1*mt**2*s/(3*p1q1**2)-128*a12*mb*mt**3*s/(3*p1q1**2)-
579  &152*a1*s/(3*p1q1)+152*a12*mb*mt*s/(3*p1q1)+
580  &128*a1*a2*mb*mt*s/(3*p1q1)+112*a1*a2*mt**2*s/(3*p1q1)-
581  &16*a12*p1p2*s/p1q1+152*a1*a2*p1p2*s/(3*p1q1)-
582  &128*a1*mt**2*s/(3*p1q2**2)-128*a12*mb*mt**3*s/(3*p1q2**2)-
583  &152*a1*s/(3*p1q2)+152*a12*mb*mt*s/(3*p1q2)+
584  &128*a1*a2*mb*mt*s/(3*p1q2)+112*a1*a2*mt**2*s/(3*p1q2)-
585  &16*a12*p1p2*s/p1q2+152*a1*a2*p1p2*s/(3*p1q2)-
586  &16*a1*mb*mt*s/(3*p1q1*p1q2)+32*a12*mb*mt**3*s/(3*p1q1*p1q2)
587  v18bis=v18bis-16*a1*p1p2*s/(3*p1q1*p1q2)+
588  &272*a1*a2*p1q1*s/(3*p1q2)+
589  &272*a1*a2*p1q2*s/(3*p1q1)-128*a2*mb**2*s/(3*p2q1**2)-
590  &128*a2**2*mb**3*mt*s/(3*p2q1**2)+
591  &32*mb**2*mt**2*s/(3*p1q2**2*p2q1**2)+32*mb**2*s/(3*p1q2*p2q1**2)-
592  &64*a2*mb**3*mt*s/(3*p1q2*p2q1**2)-
593  &64*a2*mb**2*mt**2*s/(3*p1q2*p2q1**2)-
594  &128*a2*mb**2*p1p2*s/(3*p1q2*p2q1**2)-
595  &128*a2*mb**2*p1q1*s/(3*p1q2*p2q1**2)+
596  &128*a2**2*mb**2*p1q2*s/(3*p2q1**2)+152*a2*s/(3*p2q1)-
597  &112*a1*a2*mb**2*s/(3*p2q1)-128*a1*a2*mb*mt*s/(3*p2q1)-
598  &152*a2**2*mb*mt*s/(3*p2q1)-152*a1*a2*p1p2*s/(3*p2q1)+
599  &16*a2**2*p1p2*s/p2q1+8*a1*a2*mb**3*mt*s/(3*p1q1*p2q1)+
600  &16*a1*a2*mb**2*mt**2*s/(3*p1q1*p2q1)+
601  &8*a1*a2*mb*mt**3*s/(3*p1q1*p2q1)-8*a1*p1p2*s/(3*p1q1*p2q1)-
602  &8*a2*p1p2*s/(3*p1q1*p2q1)+8*a1*a2*mb**2*p1p2*s/(3*p1q1*p2q1)+
603  &16*a1*a2*mb*mt*p1p2*s/(3*p1q1*p2q1)
604  v18bis=v18bis+8*a1*a2*mt**2*p1p2*s/(3*p1q1*p2q1)+
605  &32*a1*a2*p1p2**2*s/(3*p1q1*p2q1)-32*a2**2*p1q1*s/(3*p2q1)-
606  &32*mt**2*s/(3*p1q2**2*p2q1)+64*a1*mb**2*mt**2*s/(3*p1q2**2*p2q1)+
607  &64*a1*mb*mt**3*s/(3*p1q2**2*p2q1)+
608  &128*a1*mt**2*p1p2*s/(3*p1q2**2*p2q1)-12*s/(p1q2*p2q1)+
609  &24*a1*mb**2*s/(p1q2*p2q1)-64*a1*a2*mb**3*mt*s/(3*p1q2*p2q1)+
610  &24*a2*mt**2*s/(p1q2*p2q1)-128*a1*a2*mb**2*mt**2*s/(3*p1q2*p2q1)-
611  &64*a1*a2*mb*mt**3*s/(3*p1q2*p2q1)+56*a1*p1p2*s/(3*p1q2*p2q1)+
612  &56*a2*p1p2*s/(3*p1q2*p2q1)-64*a1*a2*mb**2*p1p2*s/(3*p1q2*p2q1)-
613  &128*a1*a2*mb*mt*p1p2*s/(3*p1q2*p2q1)-
614  &64*a1*a2*mt**2*p1p2*s/(3*p1q2*p2q1)-
615  &256*a1*a2*p1p2**2*s/(3*p1q2*p2q1)+4*p1p2*s/(3*p1q1*p1q2*p2q1)+
616  &8*a1*mb*mt*p1p2*s/(3*p1q1*p1q2*p2q1)-
617  &8*a1*mt**2*p1p2*s/(3*p1q1*p1q2*p2q1)+136*a2*p1q1*s/(3*p1q2*p2q1)-
618  &128*a1*a2*mb**2*p1q1*s/(3*p1q2*p2q1)-
619  &128*a1*a2*mb*mt*p1q1*s/(3*p1q2*p2q1)-
620  &256*a1*a2*p1p2*p1q1*s/(3*p1q2*p2q1)-160*a2**2*p1q2*s/(3*p2q1)
621  v18bis=v18bis+16*a1*a2*p1p2*p1q2*s/(3*p1q1*p2q1)-
622  &32*a12*p2q1*s/(3*p1q1)-
623  &128*a12*mt**2*p2q1*s/(3*p1q2**2)-160*a12*p2q1*s/(3*p1q2)-
624  &128*a2*mb**2*s/(3*p2q2**2)-128*a2**2*mb**3*mt*s/(3*p2q2**2)+
625  &32*mb**2*mt**2*s/(3*p1q1**2*p2q2**2)+32*mb**2*s/(3*p1q1*p2q2**2)-
626  &64*a2*mb**3*mt*s/(3*p1q1*p2q2**2)-
627  &64*a2*mb**2*mt**2*s/(3*p1q1*p2q2**2)-
628  &128*a2*mb**2*p1p2*s/(3*p1q1*p2q2**2)+
629  &128*a2**2*mb**2*p1q1*s/(3*p2q2**2)-
630  &128*a2*mb**2*p1q2*s/(3*p1q1*p2q2**2)+152*a2*s/(3*p2q2)-
631  &112*a1*a2*mb**2*s/(3*p2q2)-128*a1*a2*mb*mt*s/(3*p2q2)-
632  &152*a2**2*mb*mt*s/(3*p2q2)-152*a1*a2*p1p2*s/(3*p2q2)+
633  &16*a2**2*p1p2*s/p2q2-32*mt**2*s/(3*p1q1**2*p2q2)+
634  &64*a1*mb**2*mt**2*s/(3*p1q1**2*p2q2)+
635  &64*a1*mb*mt**3*s/(3*p1q1**2*p2q2)+
636  &128*a1*mt**2*p1p2*s/(3*p1q1**2*p2q2)-12*s/(p1q1*p2q2)+
637  &24*a1*mb**2*s/(p1q1*p2q2)-64*a1*a2*mb**3*mt*s/(3*p1q1*p2q2)
638  v18bis=v18bis+24*a2*mt**2*s/(p1q1*p2q2)-
639  &128*a1*a2*mb**2*mt**2*s/(3*p1q1*p2q2)-
640  &64*a1*a2*mb*mt**3*s/(3*p1q1*p2q2)+56*a1*p1p2*s/(3*p1q1*p2q2)+
641  &56*a2*p1p2*s/(3*p1q1*p2q2)-64*a1*a2*mb**2*p1p2*s/(3*p1q1*p2q2)-
642  &128*a1*a2*mb*mt*p1p2*s/(3*p1q1*p2q2)-
643  &64*a1*a2*mt**2*p1p2*s/(3*p1q1*p2q2)-
644  &256*a1*a2*p1p2**2*s/(3*p1q1*p2q2)-160*a2**2*p1q1*s/(3*p2q2)+
645  &8*a1*a2*mb**3*mt*s/(3*p1q2*p2q2)+
646  &16*a1*a2*mb**2*mt**2*s/(3*p1q2*p2q2)+
647  &8*a1*a2*mb*mt**3*s/(3*p1q2*p2q2)-8*a1*p1p2*s/(3*p1q2*p2q2)-
648  &8*a2*p1p2*s/(3*p1q2*p2q2)+8*a1*a2*mb**2*p1p2*s/(3*p1q2*p2q2)+
649  &16*a1*a2*mb*mt*p1p2*s/(3*p1q2*p2q2)+
650  &8*a1*a2*mt**2*p1p2*s/(3*p1q2*p2q2)+
651  &32*a1*a2*p1p2**2*s/(3*p1q2*p2q2)+4*p1p2*s/(3*p1q1*p1q2*p2q2)+
652  &8*a1*mb*mt*p1p2*s/(3*p1q1*p1q2*p2q2)-
653  &8*a1*mt**2*p1p2*s/(3*p1q1*p1q2*p2q2)+
654  &16*a1*a2*p1p2*p1q1*s/(3*p1q2*p2q2)-32*a2**2*p1q2*s/(3*p2q2)
655  v18bis=v18bis+136*a2*p1q2*s/(3*p1q1*p2q2)-
656  &128*a1*a2*mb**2*p1q2*s/(3*p1q1*p2q2)-
657  &128*a1*a2*mb*mt*p1q2*s/(3*p1q1*p2q2)-
658  &256*a1*a2*p1p2*p1q2*s/(3*p1q1*p2q2)-16*a2*mb*mt*s/(3*p2q1*p2q2)+
659  &32*a2**2*mb**3*mt*s/(3*p2q1*p2q2)-16*a2*p1p2*s/(3*p2q1*p2q2)-
660  &4*p1p2*s/(3*p1q1*p2q1*p2q2)+8*a2*mb**2*p1p2*s/(3*p1q1*p2q1*p2q2)-
661  &8*a2*mb*mt*p1p2*s/(3*p1q1*p2q1*p2q2)-4*p1p2*s/(3*p1q2*p2q1*p2q2)+
662  &8*a2*mb**2*p1p2*s/(3*p1q2*p2q1*p2q2)-
663  &8*a2*mb*mt*p1p2*s/(3*p1q2*p2q1*p2q2)+
664  &2*mb**3*mt*s/(3*p1q1*p1q2*p2q1*p2q2)+
665  &4*mb**2*mt**2*s/(3*p1q1*p1q2*p2q1*p2q2)+
666  &2*mb*mt**3*s/(3*p1q1*p1q2*p2q1*p2q2)-
667  &2*mb**2*p1p2*s/(3*p1q1*p1q2*p2q1*p2q2)-
668  &4*mb*mt*p1p2*s/(3*p1q1*p1q2*p2q1*p2q2)-
669  &2*mt**2*p1p2*s/(3*p1q1*p1q2*p2q1*p2q2)-
670  &8*p1p2**2*s/(3*p1q1*p1q2*p2q1*p2q2)+
671  &8*a2*p1p2*p1q1*s/(3*p1q2*p2q1*p2q2)
672  v18bis=v18bis+8*a2*p1p2*p1q2*s/(3*p1q1*p2q1*p2q2)+
673  &272*a1*a2*p2q1*s/(3*p2q2)-
674  &128*a1*mt**2*p2q1*s/(3*p1q1**2*p2q2)-136*a1*p2q1*s/(3*p1q1*p2q2)+
675  &128*a1*a2*mb*mt*p2q1*s/(3*p1q1*p2q2)+
676  &128*a1*a2*mt**2*p2q1*s/(3*p1q1*p2q2)+
677  &256*a1*a2*p1p2*p2q1*s/(3*p1q1*p2q2)-
678  &16*a1*a2*p1p2*p2q1*s/(3*p1q2*p2q2)+
679  &8*a1*p1p2*p2q1*s/(3*p1q1*p1q2*p2q2)+
680  &256*a1*a2*p1q2*p2q1*s/(3*p1q1*p2q2)-
681  &128*a12*mt**2*p2q2*s/(3*p1q1**2)-160*a12*p2q2*s/(3*p1q1)-
682  &32*a12*p2q2*s/(3*p1q2)+272*a1*a2*p2q2*s/(3*p2q1)-
683  &16*a1*a2*p1p2*p2q2*s/(3*p1q1*p2q1)-
684  &128*a1*mt**2*p2q2*s/(3*p1q2**2*p2q1)-136*a1*p2q2*s/(3*p1q2*p2q1)+
685  &128*a1*a2*mb*mt*p2q2*s/(3*p1q2*p2q1)+
686  &128*a1*a2*mt**2*p2q2*s/(3*p1q2*p2q1)+
687  &256*a1*a2*p1p2*p2q2*s/(3*p1q2*p2q1)+
688  &8*a1*p1p2*p2q2*s/(3*p1q1*p1q2*p2q1)
689  v18bis=v18bis+256*a1*a2*p1q1*p2q2*s/(3*p1q2*p2q1)+
690  &8*a12*mb*mt*s**2/(3*p1q1*p1q2)+16*a12*p1p2*s**2/(3*p1q1*p1q2)-
691  &8*a1*a2*p1p2*s**2/(3*p1q1*p2q1)+4*a1*p1p2*s**2/(3*p1q1*p1q2*p2q1)-
692  &8*a1*a2*p1p2*s**2/(3*p1q2*p2q2)+4*a1*p1p2*s**2/(3*p1q1*p1q2*p2q2)+
693  &8*a2**2*mb*mt*s**2/(3*p2q1*p2q2)+16*a2**2*p1p2*s**2/(3*p2q1*p2q2)-
694  &4*a2*p1p2*s**2/(3*p1q1*p2q1*p2q2)-
695  &4*a2*p1p2*s**2/(3*p1q2*p2q1*p2q2)+
696  &2*p1p2*s**2/(3*p1q1*p1q2*p2q1*p2q2)
697 C
698 
699  a18 = 640*a1/3+640*a2/3+32*a1*a2*mb**2+368*a12*mb*mt+
700  &512*a1*a2*mb*mt/3+
701  &368*a2**2*mb*mt+32*a1*a2*mt**2+496*a12*p1p2/3+
702  &320*a1*a2*p1p2+496*a2**2*p1p2/3-128*a1*mb*mt**3/(3*p1q1**2)+
703  &128*a1*mt**4/(3*p1q1**2)+256*a12*mb*mt**5/(3*p1q1**2)+
704  &256*a1*mt**2*p1p2/(3*p1q1**2)-256*a12*mt**4*p1p2/(3*p1q1**2)+
705  &8/(3*p1q1)+32*a1*mb*mt/p1q1+56*a2*mb*mt/(3*p1q1)+
706  &88*a1*mt**2/(3*p1q1)+72*a2*mt**2/p1q1-
707  &704*a12*mb*mt**3/(3*p1q1)+224*a1*a2*mb*mt**3/(3*p1q1)+
708  &104*a1*p1p2/(3*p1q1)+48*a2*p1p2/p1q1-
709  &128*a1*a2*mb*mt*p1p2/(3*p1q1)+512*a12*mt**2*p1p2/(3*p1q1)-
710  &448*a1*a2*mt**2*p1p2/(3*p1q1)-32*a1*a2*p1p2**2/p1q1-
711  &656*a1*a2*p1q1/3-224*a2**2*p1q1-128*a1*mb*mt**3/(3*p1q2**2)+
712  &128*a1*mt**4/(3*p1q2**2)+256*a12*mb*mt**5/(3*p1q2**2)+
713  &256*a1*mt**2*p1p2/(3*p1q2**2)-256*a12*mt**4*p1p2/(3*p1q2**2)+
714  &256*a1*mt**2*p1q1/(3*p1q2**2)-256*a12*mb*mt**3*p1q1/(3*p1q2**2)+
715  &8/(3*p1q2)+32*a1*mb*mt/p1q2+56*a2*mb*mt/(3*p1q2)
716  a18=a18+88*a1*mt**2/(3*p1q2)+72*a2*mt**2/p1q2-
717  &704*a12*mb*mt**3/(3*p1q2)+224*a1*a2*mb*mt**3/(3*p1q2)+
718  &104*a1*p1p2/(3*p1q2)+48*a2*p1p2/p1q2-
719  &128*a1*a2*mb*mt*p1p2/(3*p1q2)+512*a12*mt**2*p1p2/(3*p1q2)-
720  &448*a1*a2*mt**2*p1p2/(3*p1q2)-32*a1*a2*p1p2**2/p1q2+
721  &32*a1*mb*mt**3/(3*p1q1*p1q2)-32*a1*mt**4/(3*p1q1*p1q2)-
722  &64*a12*mb*mt**5/(3*p1q1*p1q2)+16*p1p2/(3*p1q1*p1q2)-
723  &64*a1*mt**2*p1p2/(3*p1q1*p1q2)+64*a12*mt**4*p1p2/(3*p1q1*p1q2)+
724  &112*a1*p1q1/p1q2+272*a2*p1q1/(3*p1q2)-
725  &272*a1*a2*mb**2*p1q1/(3*p1q2)-208*a12*mb*mt*p1q1/(3*p1q2)+
726  &400*a1*a2*mb*mt*p1q1/(3*p1q2)-80*a1*a2*mt**2*p1q1/p1q2+
727  &96*a12*p1p2*p1q1/p1q2-320*a1*a2*p1p2*p1q1/p1q2-
728  &544*a1*a2*p1q1**2/(3*p1q2)-656*a1*a2*p1q2/3-224*a2**2*p1q2+
729  &256*a1*mt**2*p1q2/(3*p1q1**2)-256*a12*mb*mt**3*p1q2/(3*p1q1**2)+
730  &112*a1*p1q2/p1q1+272*a2*p1q2/(3*p1q1)-
731  &272*a1*a2*mb**2*p1q2/(3*p1q1)-208*a12*mb*mt*p1q2/(3*p1q1)+
732  &400*a1*a2*mb*mt*p1q2/(3*p1q1)-80*a1*a2*mt**2*p1q2/p1q1
733  a18=a18+96*a12*p1p2*p1q2/p1q1-320*a1*a2*p1p2*p1q2/p1q1-
734  &544*a1*a2*p1q2**2/(3*p1q1)+128*a2*mb**4/(3*p2q1**2)-
735  &128*a2*mb**3*mt/(3*p2q1**2)+256*a2**2*mb**5*mt/(3*p2q1**2)+
736  &256*a2*mb**2*p1p2/(3*p2q1**2)-256*a2**2*mb**4*p1p2/(3*p2q1**2)+
737  &256*a2*mb**2*p1q1/(3*p2q1**2)-256*a2**2*mb**4*p1q1/(3*p2q1**2)+
738  &64*mb**3*mt**3/(3*p1q2**2*p2q1**2)-
739  &64*mb**2*mt**2*p1p2/(3*p1q2**2*p2q1**2)-
740  &64*mb**2*mt**2*p1q1/(3*p1q2**2*p2q1**2)-
741  &64*mb**3*mt/(3*p1q2*p2q1**2)-
742  &256*a2*mb**3*mt*p1p2/(3*p1q2*p2q1**2)+
743  &256*a2*mb**2*p1p2**2/(3*p1q2*p2q1**2)-
744  &256*a2*mb**3*mt*p1q1/(3*p1q2*p2q1**2)+
745  &512*a2*mb**2*p1p2*p1q1/(3*p1q2*p2q1**2)+
746  &256*a2*mb**2*p1q1**2/(3*p1q2*p2q1**2)-
747  &256*a2**2*mb**4*p1q2/(3*p2q1**2)-8/(3*p2q1)-72*a1*mb**2/p2q1-
748  &88*a2*mb**2/(3*p2q1)-56*a1*mb*mt/(3*p2q1)-32*a2*mb*mt/p2q1-
749  &224*a1*a2*mb**3*mt/(3*p2q1)+704*a2**2*mb**3*mt/(3*p2q1)
750  a18=a18-48*a1*p1p2/p2q1-104*a2*p1p2/(3*p2q1)+
751  &448*a1*a2*mb**2*p1p2/(3*p2q1)-512*a2**2*mb**2*p1p2/(3*p2q1)+
752  &128*a1*a2*mb*mt*p1p2/(3*p2q1)+32*a1*a2*p1p2**2/p2q1-
753  &16*p1p2/(3*p1q1*p2q1)+32*a1*mb*mt*p1p2/(3*p1q1*p2q1)+
754  &32*a2*mb*mt*p1p2/(3*p1q1*p2q1)+
755  &64*a1*a2*mb*mt*p1p2**2/(3*p1q1*p2q1)-
756  &64*a1*a2*p1p2**3/(3*p1q1*p2q1)-256*a2*p1q1/(3*p2q1)+
757  &448*a1*a2*mb**2*p1q1/(3*p2q1)-368*a2**2*mb**2*p1q1/(3*p2q1)-
758  &224*a1*a2*mb*mt*p1q1/(3*p2q1)+304*a1*a2*p1p2*p1q1/(3*p2q1)+
759  &64*mb*mt**3/(3*p1q2**2*p2q1)+
760  &256*a1*mb*mt**3*p1p2/(3*p1q2**2*p2q1)-
761  &256*a1*mt**2*p1p2**2/(3*p1q2**2*p2q1)+
762  &64*mt**2*p1q1/(3*p1q2**2*p2q1)-
763  &128*a1*mb**2*mt**2*p1q1/(3*p1q2**2*p2q1)+
764  &128*a1*mb*mt**3*p1q1/(3*p1q2**2*p2q1)-
765  &256*a1*mt**2*p1p2*p1q1/(3*p1q2**2*p2q1)-4*mb**2/(3*p1q2*p2q1)-
766  &64*mb*mt/(3*p1q2*p2q1)+128*a2*mb**3*mt/(3*p1q2*p2q1)
767  a18=a18-4*mt**2/(3*p1q2*p2q1)-128*a1*mb**2*mt**2/(3*p1q2*p2q1)-
768  &128*a2*mb**2*mt**2/(3*p1q2*p2q1)+128*a1*mb*mt**3/(3*p1q2*p2q1)-
769  &112*a2*mb**2*p1p2/(3*p1q2*p2q1)+32*a1*mb*mt*p1p2/(3*p1q2*p2q1)+
770  &32*a2*mb*mt*p1p2/(3*p1q2*p2q1)-112*a1*mt**2*p1p2/(3*p1q2*p2q1)-
771  &48*a1*p1p2**2/(p1q2*p2q1)-48*a2*p1p2**2/(p1q2*p2q1)-
772  &512*a1*a2*mb*mt*p1p2**2/(3*p1q2*p2q1)+
773  &512*a1*a2*p1p2**3/(3*p1q2*p2q1)+8*mb*mt*p1p2/(3*p1q1*p1q2*p2q1)-
774  &8*mt**2*p1p2/(3*p1q1*p1q2*p2q1)-
775  &32*a1*mb*mt**3*p1p2/(3*p1q1*p1q2*p2q1)-
776  &16*p1p2**2/(3*p1q1*p1q2*p2q1)+
777  &32*a1*mt**2*p1p2**2/(3*p1q1*p1q2*p2q1)+8*p1q1/(3*p1q2*p2q1)-
778  &160*a1*mb**2*p1q1/(3*p1q2*p2q1)-272*a2*mb**2*p1q1/(3*p1q2*p2q1)-
779  &56*a1*mb*mt*p1q1/(3*p1q2*p2q1)-200*a2*mb*mt*p1q1/(3*p1q2*p2q1)-
780  &48*a1*p1p2*p1q1/(p1q2*p2q1)-256*a2*p1p2*p1q1/(3*p1q2*p2q1)+
781  &256*a1*a2*mb**2*p1p2*p1q1/(3*p1q2*p2q1)-
782  &256*a1*a2*mb*mt*p1p2*p1q1/(p1q2*p2q1)+
783  &1024*a1*a2*p1p2**2*p1q1/(3*p1q2*p2q1)
784  a18=a18-272*a2*p1q1**2/(3*p1q2*p2q1)+
785  &256*a1*a2*mb**2*p1q1**2/(3*p1q2*p2q1)-
786  &256*a1*a2*mb*mt*p1q1**2/(3*p1q2*p2q1)+
787  &512*a1*a2*p1p2*p1q1**2/(3*p1q2*p2q1)+16*a2*p1q2/(3*p2q1)+
788  &64*a1*a2*mb**2*p1q2/p2q1+32*a2**2*mb**2*p1q2/(3*p2q1)-
789  &112*a1*a2*mb*mt*p1q2/(3*p2q1)+368*a1*a2*p1p2*p1q2/(3*p2q1)+
790  &32*a2*p1p2*p1q2/(3*p1q1*p2q1)-
791  &32*a1*a2*mb**2*p1p2*p1q2/(3*p1q1*p2q1)+
792  &32*a1*a2*mb*mt*p1p2*p1q2/(3*p1q1*p2q1)-
793  &64*a1*a2*p1p2**2*p1q2/(3*p1q1*p2q1)+224*a12*p2q1+
794  &656*a1*a2*p2q1/3-256*a1*mt**2*p2q1/(3*p1q1**2)+
795  &256*a12*mt**4*p2q1/(3*p1q1**2)-256*a1*p2q1/(3*p1q1)-
796  &224*a1*a2*mb*mt*p2q1/(3*p1q1)-368*a12*mt**2*p2q1/(3*p1q1)+
797  &448*a1*a2*mt**2*p2q1/(3*p1q1)+304*a1*a2*p1p2*p2q1/(3*p1q1)+
798  &256*a12*mt**4*p2q1/(3*p1q2**2)+
799  &256*a12*mt**2*p1q1*p2q1/(3*p1q2**2)+16*a1*p2q1/(3*p1q2)-
800  &112*a1*a2*mb*mt*p2q1/(3*p1q2)+32*a12*mt**2*p2q1/(3*p1q2)
801  a18=a18+64*a1*a2*mt**2*p2q1/p1q2+368*a1*a2*p1p2*p2q1/(3*p1q2)+
802  &16*a1*mt**2*p2q1/(3*p1q1*p1q2)-64*a12*mt**4*p2q1/(3*p1q1*p1q2)+
803  &640*a12*p1q1*p2q1/(3*p1q2)+544*a1*a2*p1q1*p2q1/(3*p1q2)+
804  &32*a12*p1q2*p2q1/p1q1+944*a1*a2*p1q2*p2q1/(3*p1q1)+
805  &128*a2*mb**4/(3*p2q2**2)-128*a2*mb**3*mt/(3*p2q2**2)+
806  &256*a2**2*mb**5*mt/(3*p2q2**2)+256*a2*mb**2*p1p2/(3*p2q2**2)-
807  &256*a2**2*mb**4*p1p2/(3*p2q2**2)+
808  &64*mb**3*mt**3/(3*p1q1**2*p2q2**2)-
809  &64*mb**2*mt**2*p1p2/(3*p1q1**2*p2q2**2)-
810  &64*mb**3*mt/(3*p1q1*p2q2**2)-
811  &256*a2*mb**3*mt*p1p2/(3*p1q1*p2q2**2)+
812  &256*a2*mb**2*p1p2**2/(3*p1q1*p2q2**2)-
813  &256*a2**2*mb**4*p1q1/(3*p2q2**2)+256*a2*mb**2*p1q2/(3*p2q2**2)-
814  &256*a2**2*mb**4*p1q2/(3*p2q2**2)-
815  &64*mb**2*mt**2*p1q2/(3*p1q1**2*p2q2**2)-
816  &256*a2*mb**3*mt*p1q2/(3*p1q1*p2q2**2)+
817  &512*a2*mb**2*p1p2*p1q2/(3*p1q1*p2q2**2)
818  a18=a18+256*a2*mb**2*p1q2**2/(3*p1q1*p2q2**2)-
819  &256*a2*mb**2*p2q1/(3*p2q2**2)+256*a2**2*mb**3*mt*p2q1/(3*p2q2**2)+
820  &64*mb**2*mt**2*p2q1/(3*p1q1**2*p2q2**2)+
821  &64*mb**2*p2q1/(3*p1q1*p2q2**2)+
822  &128*a2*mb**3*mt*p2q1/(3*p1q1*p2q2**2)-
823  &128*a2*mb**2*mt**2*p2q1/(3*p1q1*p2q2**2)-
824  &256*a2*mb**2*p1p2*p2q1/(3*p1q1*p2q2**2)+
825  &256*a2**2*mb**2*p1q1*p2q1/(3*p2q2**2)-
826  &256*a2*mb**2*p1q2*p2q1/(3*p1q1*p2q2**2)-8/(3*p2q2)-
827  &72*a1*mb**2/p2q2-88*a2*mb**2/(3*p2q2)-56*a1*mb*mt/(3*p2q2)-
828  &32*a2*mb*mt/p2q2-224*a1*a2*mb**3*mt/(3*p2q2)+
829  &704*a2**2*mb**3*mt/(3*p2q2)-48*a1*p1p2/p2q2-
830  &104*a2*p1p2/(3*p2q2)+448*a1*a2*mb**2*p1p2/(3*p2q2)-
831  &512*a2**2*mb**2*p1p2/(3*p2q2)+128*a1*a2*mb*mt*p1p2/(3*p2q2)+
832  &32*a1*a2*p1p2**2/p2q2+64*mb*mt**3/(3*p1q1**2*p2q2)+
833  &256*a1*mb*mt**3*p1p2/(3*p1q1**2*p2q2)-
834  &256*a1*mt**2*p1p2**2/(3*p1q1**2*p2q2)-4*mb**2/(3*p1q1*p2q2)
835  a18=a18-64*mb*mt/(3*p1q1*p2q2)+128*a2*mb**3*mt/(3*p1q1*p2q2)-
836  &4*mt**2/(3*p1q1*p2q2)-128*a1*mb**2*mt**2/(3*p1q1*p2q2)-
837  &128*a2*mb**2*mt**2/(3*p1q1*p2q2)+128*a1*mb*mt**3/(3*p1q1*p2q2)-
838  &112*a2*mb**2*p1p2/(3*p1q1*p2q2)+32*a1*mb*mt*p1p2/(3*p1q1*p2q2)+
839  &32*a2*mb*mt*p1p2/(3*p1q1*p2q2)-112*a1*mt**2*p1p2/(3*p1q1*p2q2)-
840  &48*a1*p1p2**2/(p1q1*p2q2)-48*a2*p1p2**2/(p1q1*p2q2)-
841  &512*a1*a2*mb*mt*p1p2**2/(3*p1q1*p2q2)+
842  &512*a1*a2*p1p2**3/(3*p1q1*p2q2)+16*a2*p1q1/(3*p2q2)+
843  &64*a1*a2*mb**2*p1q1/p2q2+32*a2**2*mb**2*p1q1/(3*p2q2)-
844  &112*a1*a2*mb*mt*p1q1/(3*p2q2)+368*a1*a2*p1p2*p1q1/(3*p2q2)-
845  &16*p1p2/(3*p1q2*p2q2)+32*a1*mb*mt*p1p2/(3*p1q2*p2q2)+
846  &32*a2*mb*mt*p1p2/(3*p1q2*p2q2)+
847  &64*a1*a2*mb*mt*p1p2**2/(3*p1q2*p2q2)-
848  &64*a1*a2*p1p2**3/(3*p1q2*p2q2)+8*mb*mt*p1p2/(3*p1q1*p1q2*p2q2)-
849  &8*mt**2*p1p2/(3*p1q1*p1q2*p2q2)-
850  &32*a1*mb*mt**3*p1p2/(3*p1q1*p1q2*p2q2)-
851  &16*p1p2**2/(3*p1q1*p1q2*p2q2)
852  a18=a18+32*a1*mt**2*p1p2**2/(3*p1q1*p1q2*p2q2)+
853  &32*a2*p1p2*p1q1/(3*p1q2*p2q2)-
854  &32*a1*a2*mb**2*p1p2*p1q1/(3*p1q2*p2q2)+
855  &32*a1*a2*mb*mt*p1p2*p1q1/(3*p1q2*p2q2)-
856  &64*a1*a2*p1p2**2*p1q1/(3*p1q2*p2q2)-256*a2*p1q2/(3*p2q2)+
857  &448*a1*a2*mb**2*p1q2/(3*p2q2)-368*a2**2*mb**2*p1q2/(3*p2q2)-
858  &224*a1*a2*mb*mt*p1q2/(3*p2q2)+304*a1*a2*p1p2*p1q2/(3*p2q2)+
859  &64*mt**2*p1q2/(3*p1q1**2*p2q2)-
860  &128*a1*mb**2*mt**2*p1q2/(3*p1q1**2*p2q2)+
861  &128*a1*mb*mt**3*p1q2/(3*p1q1**2*p2q2)-
862  &256*a1*mt**2*p1p2*p1q2/(3*p1q1**2*p2q2)+8*p1q2/(3*p1q1*p2q2)-
863  &160*a1*mb**2*p1q2/(3*p1q1*p2q2)-272*a2*mb**2*p1q2/(3*p1q1*p2q2)-
864  &56*a1*mb*mt*p1q2/(3*p1q1*p2q2)-200*a2*mb*mt*p1q2/(3*p1q1*p2q2)-
865  &48*a1*p1p2*p1q2/(p1q1*p2q2)-256*a2*p1p2*p1q2/(3*p1q1*p2q2)+
866  &256*a1*a2*mb**2*p1p2*p1q2/(3*p1q1*p2q2)-
867  &256*a1*a2*mb*mt*p1p2*p1q2/(p1q1*p2q2)+
868  &1024*a1*a2*p1p2**2*p1q2/(3*p1q1*p2q2)
869  a18=a18-272*a2*p1q2**2/(3*p1q1*p2q2)+
870  &256*a1*a2*mb**2*p1q2**2/(3*p1q1*p2q2)-
871  &256*a1*a2*mb*mt*p1q2**2/(3*p1q1*p2q2)+
872  &512*a1*a2*p1p2*p1q2**2/(3*p1q1*p2q2)-32*a2*mb**4/(3*p2q1*p2q2)+
873  &32*a2*mb**3*mt/(3*p2q1*p2q2)-64*a2**2*mb**5*mt/(3*p2q1*p2q2)+
874  &16*p1p2/(3*p2q1*p2q2)-64*a2*mb**2*p1p2/(3*p2q1*p2q2)+
875  &64*a2**2*mb**4*p1p2/(3*p2q1*p2q2)+8*mb**2*p1p2/(3*p1q1*p2q1*p2q2)-
876  &8*mb*mt*p1p2/(3*p1q1*p2q1*p2q2)+
877  &32*a2*mb**3*mt*p1p2/(3*p1q1*p2q1*p2q2)+
878  &16*p1p2**2/(3*p1q1*p2q1*p2q2)-
879  &32*a2*mb**2*p1p2**2/(3*p1q1*p2q1*p2q2)-
880  &16*a2*mb**2*p1q1/(3*p2q1*p2q2)+64*a2**2*mb**4*p1q1/(3*p2q1*p2q2)+
881  &8*mb**2*p1p2/(3*p1q2*p2q1*p2q2)-8*mb*mt*p1p2/(3*p1q2*p2q1*p2q2)+
882  &32*a2*mb**3*mt*p1p2/(3*p1q2*p2q1*p2q2)+
883  &16*p1p2**2/(3*p1q2*p2q1*p2q2)-
884  &32*a2*mb**2*p1p2**2/(3*p1q2*p2q1*p2q2)-
885  &16*mb*mt*p1p2**2/(3*p1q1*p1q2*p2q1*p2q2)
886  a18=a18+16*p1p2**3/(3*p1q1*p1q2*p2q1*p2q2)-
887  &32*a2*mb**2*p1p2*p1q1/(3*p1q2*p2q1*p2q2)-
888  &16*a2*mb**2*p1q2/(3*p2q1*p2q2)+64*a2**2*mb**4*p1q2/(3*p2q1*p2q2)-
889  &32*a2*mb**2*p1p2*p1q2/(3*p1q1*p2q1*p2q2)+272*a1*p2q1/(3*p2q2)+
890  &112*a2*p2q1/p2q2-80*a1*a2*mb**2*p2q1/p2q2+
891  &400*a1*a2*mb*mt*p2q1/(3*p2q2)-208*a2**2*mb*mt*p2q1/(3*p2q2)-
892  &272*a1*a2*mt**2*p2q1/(3*p2q2)-320*a1*a2*p1p2*p2q1/p2q2+
893  &96*a2**2*p1p2*p2q1/p2q2-256*a1*mb*mt**3*p2q1/(3*p1q1**2*p2q2)+
894  &512*a1*mt**2*p1p2*p2q1/(3*p1q1**2*p2q2)-8*p2q1/(3*p1q1*p2q2)+
895  &200*a1*mb*mt*p2q1/(3*p1q1*p2q2)+56*a2*mb*mt*p2q1/(3*p1q1*p2q2)+
896  &272*a1*mt**2*p2q1/(3*p1q1*p2q2)+160*a2*mt**2*p2q1/(3*p1q1*p2q2)+
897  &256*a1*p1p2*p2q1/(3*p1q1*p2q2)+48*a2*p1p2*p2q1/(p1q1*p2q2)+
898  &256*a1*a2*mb*mt*p1p2*p2q1/(p1q1*p2q2)-
899  &256*a1*a2*mt**2*p1p2*p2q1/(3*p1q1*p2q2)-
900  &1024*a1*a2*p1p2**2*p2q1/(3*p1q1*p2q2)-
901  &544*a1*a2*p1q1*p2q1/(3*p2q2)-640*a2**2*p1q1*p2q1/(3*p2q2)-
902  &32*a1*p1p2*p2q1/(3*p1q2*p2q2)
903  a18=a18-32*a1*a2*mb*mt*p1p2*p2q1/(3*p1q2*p2q2)+
904  &32*a1*a2*mt**2*p1p2*p2q1/(3*p1q2*p2q2)+
905  &64*a1*a2*p1p2**2*p2q1/(3*p1q2*p2q2)-
906  &32*a1*mt**2*p1p2*p2q1/(3*p1q1*p1q2*p2q2)+
907  &64*a1*a2*p1p2*p1q1*p2q1/(3*p1q2*p2q2)-
908  &944*a1*a2*p1q2*p2q1/(3*p2q2)-32*a2**2*p1q2*p2q1/p2q2+
909  &256*a1*mt**2*p1q2*p2q1/(3*p1q1**2*p2q2)+
910  &96*a1*p1q2*p2q1/(p1q1*p2q2)+96*a2*p1q2*p2q1/(p1q1*p2q2)-
911  &128*a1*a2*mb**2*p1q2*p2q1/(3*p1q1*p2q2)+
912  &256*a1*a2*mb*mt*p1q2*p2q1/(p1q1*p2q2)-
913  &128*a1*a2*mt**2*p1q2*p2q1/(3*p1q1*p2q2)-
914  &512*a1*a2*p1p2*p1q2*p2q1/(p1q1*p2q2)-
915  &512*a1*a2*p1q2**2*p2q1/(3*p1q1*p2q2)+544*a1*a2*p2q1**2/(3*p2q2)-
916  &256*a1*mt**2*p2q1**2/(3*p1q1**2*p2q2)-
917  &272*a1*p2q1**2/(3*p1q1*p2q2)-
918  &256*a1*a2*mb*mt*p2q1**2/(3*p1q1*p2q2)+
919  &256*a1*a2*mt**2*p2q1**2/(3*p1q1*p2q2)
920  a18=a18+512*a1*a2*p1p2*p2q1**2/(3*p1q1*p2q2)+
921  &512*a1*a2*p1q2*p2q1**2/(3*p1q1*p2q2)+224*a12*p2q2+
922  &656*a1*a2*p2q2/3+256*a12*mt**4*p2q2/(3*p1q1**2)+
923  &16*a1*p2q2/(3*p1q1)-112*a1*a2*mb*mt*p2q2/(3*p1q1)+
924  &32*a12*mt**2*p2q2/(3*p1q1)+64*a1*a2*mt**2*p2q2/p1q1+
925  &368*a1*a2*p1p2*p2q2/(3*p1q1)-256*a1*mt**2*p2q2/(3*p1q2**2)+
926  &256*a12*mt**4*p2q2/(3*p1q2**2)-256*a1*p2q2/(3*p1q2)-
927  &224*a1*a2*mb*mt*p2q2/(3*p1q2)-368*a12*mt**2*p2q2/(3*p1q2)+
928  &448*a1*a2*mt**2*p2q2/(3*p1q2)+304*a1*a2*p1p2*p2q2/(3*p1q2)+
929  &16*a1*mt**2*p2q2/(3*p1q1*p1q2)-64*a12*mt**4*p2q2/(3*p1q1*p1q2)+
930  &32*a12*p1q1*p2q2/p1q2+944*a1*a2*p1q1*p2q2/(3*p1q2)+
931  &256*a12*mt**2*p1q2*p2q2/(3*p1q1**2)+
932  &640*a12*p1q2*p2q2/(3*p1q1)+544*a1*a2*p1q2*p2q2/(3*p1q1)-
933  &256*a2*mb**2*p2q2/(3*p2q1**2)+256*a2**2*mb**3*mt*p2q2/(3*p2q1**2)+
934  &64*mb**2*mt**2*p2q2/(3*p1q2**2*p2q1**2)+
935  &64*mb**2*p2q2/(3*p1q2*p2q1**2)+
936  &128*a2*mb**3*mt*p2q2/(3*p1q2*p2q1**2)
937  a18=a18-128*a2*mb**2*mt**2*p2q2/(3*p1q2*p2q1**2)-
938  &256*a2*mb**2*p1p2*p2q2/(3*p1q2*p2q1**2)-
939  &256*a2*mb**2*p1q1*p2q2/(3*p1q2*p2q1**2)+
940  &256*a2**2*mb**2*p1q2*p2q2/(3*p2q1**2)+272*a1*p2q2/(3*p2q1)+
941  &112*a2*p2q2/p2q1-80*a1*a2*mb**2*p2q2/p2q1+
942  &400*a1*a2*mb*mt*p2q2/(3*p2q1)-208*a2**2*mb*mt*p2q2/(3*p2q1)-
943  &272*a1*a2*mt**2*p2q2/(3*p2q1)-320*a1*a2*p1p2*p2q2/p2q1+
944  &96*a2**2*p1p2*p2q2/p2q1-32*a1*p1p2*p2q2/(3*p1q1*p2q1)-
945  &32*a1*a2*mb*mt*p1p2*p2q2/(3*p1q1*p2q1)+
946  &32*a1*a2*mt**2*p1p2*p2q2/(3*p1q1*p2q1)+
947  &64*a1*a2*p1p2**2*p2q2/(3*p1q1*p2q1)-944*a1*a2*p1q1*p2q2/(3*p2q1)-
948  &32*a2**2*p1q1*p2q2/p2q1-256*a1*mb*mt**3*p2q2/(3*p1q2**2*p2q1)+
949  &512*a1*mt**2*p1p2*p2q2/(3*p1q2**2*p2q1)+
950  &256*a1*mt**2*p1q1*p2q2/(3*p1q2**2*p2q1)-8*p2q2/(3*p1q2*p2q1)+
951  &200*a1*mb*mt*p2q2/(3*p1q2*p2q1)+56*a2*mb*mt*p2q2/(3*p1q2*p2q1)+
952  &272*a1*mt**2*p2q2/(3*p1q2*p2q1)+160*a2*mt**2*p2q2/(3*p1q2*p2q1)+
953  &256*a1*p1p2*p2q2/(3*p1q2*p2q1)+48*a2*p1p2*p2q2/(p1q2*p2q1)
954  a18=a18+256*a1*a2*mb*mt*p1p2*p2q2/(p1q2*p2q1)-
955  &256*a1*a2*mt**2*p1p2*p2q2/(3*p1q2*p2q1)-
956  &1024*a1*a2*p1p2**2*p2q2/(3*p1q2*p2q1)-
957  &32*a1*mt**2*p1p2*p2q2/(3*p1q1*p1q2*p2q1)+
958  &96*a1*p1q1*p2q2/(p1q2*p2q1)+96*a2*p1q1*p2q2/(p1q2*p2q1)-
959  &128*a1*a2*mb**2*p1q1*p2q2/(3*p1q2*p2q1)+
960  &256*a1*a2*mb*mt*p1q1*p2q2/(p1q2*p2q1)-
961  &128*a1*a2*mt**2*p1q1*p2q2/(3*p1q2*p2q1)-
962  &512*a1*a2*p1p2*p1q1*p2q2/(p1q2*p2q1)-
963  &512*a1*a2*p1q1**2*p2q2/(3*p1q2*p2q1)-544*a1*a2*p1q2*p2q2/(3*p2q1)-
964  &640*a2**2*p1q2*p2q2/(3*p2q1)+
965  &64*a1*a2*p1p2*p1q2*p2q2/(3*p1q1*p2q1)+544*a1*a2*p2q2**2/(3*p2q1)-
966  &256*a1*mt**2*p2q2**2/(3*p1q2**2*p2q1)-
967  &272*a1*p2q2**2/(3*p1q2*p2q1)-
968  &256*a1*a2*mb*mt*p2q2**2/(3*p1q2*p2q1)+
969  &256*a1*a2*mt**2*p2q2**2/(3*p1q2*p2q1)+
970  &512*a1*a2*p1p2*p2q2**2/(3*p1q2*p2q1)
971  a18=a18+512*a1*a2*p1q1*p2q2**2/(3*p1q2*p2q1)-
972  &384*a12*mb*mt*p1q1**2/s**2+
973  &384*a12*p1p2*p1q1**2/s**2-2688*a12*mb*mt*p1q1*p1q2/s**2+
974  &2688*a12*p1p2*p1q1*p1q2/s**2-384*a12*mb*mt*p1q2**2/s**2+
975  &384*a12*p1p2*p1q2**2/s**2-768*a1*a2*mb*mt*p1q1*p2q1/s**2+
976  &768*a1*a2*p1p2*p1q1*p2q1/s**2-2688*a1*a2*mb*mt*p1q2*p2q1/s**2+
977  &2688*a1*a2*p1p2*p1q2*p2q1/s**2-960*a12*p1q1*p1q2*p2q1/s**2-
978  &960*a1*a2*p1q1*p1q2*p2q1/s**2+960*a12*p1q2**2*p2q1/s**2+
979  &960*a1*a2*p1q2**2*p2q1/s**2-384*a2**2*mb*mt*p2q1**2/s**2+
980  &384*a2**2*p1p2*p2q1**2/s**2-960*a1*a2*p1q2*p2q1**2/s**2-
981  &960*a2**2*p1q2*p2q1**2/s**2-2688*a1*a2*mb*mt*p1q1*p2q2/s**2+
982  &2688*a1*a2*p1p2*p1q1*p2q2/s**2+960*a12*p1q1**2*p2q2/s**2+
983  &960*a1*a2*p1q1**2*p2q2/s**2-768*a1*a2*mb*mt*p1q2*p2q2/s**2+
984  &768*a1*a2*p1p2*p1q2*p2q2/s**2-960*a12*p1q1*p1q2*p2q2/s**2-
985  &960*a1*a2*p1q1*p1q2*p2q2/s**2-2688*a2**2*mb*mt*p2q1*p2q2/s**2+
986  &2688*a2**2*p1p2*p2q1*p2q2/s**2+960*a1*a2*p1q1*p2q1*p2q2/s**2+
987  &960*a2**2*p1q1*p2q1*p2q2/s**2+960*a1*a2*p1q2*p2q1*p2q2/s**2
988  a18=a18+960*a2**2*p1q2*p2q1*p2q2/s**2-
989  &384*a2**2*mb*mt*p2q2**2/s**2+
990  &384*a2**2*p1p2*p2q2**2/s**2-960*a1*a2*p1q1*p2q2**2/s**2-
991  &960*a2**2*p1q1*p2q2**2/s**2-96*a1*mb*mt/s-96*a2*mb*mt/s+
992  &768*a2**2*mb**3*mt/s+768*a12*mb*mt**3/s-192*a1*p1p2/s-
993  &192*a2*p1p2/s-768*a2**2*mb**2*p1p2/s+2304*a1*a2*mb*mt*p1p2/s-
994  &768*a12*mt**2*p1p2/s-2304*a1*a2*p1p2**2/s+
995  &96*a1*mb*mt**3/(p1q1*s)+192*a2*mb*mt*p1p2/(p1q1*s)-
996  &96*a1*mt**2*p1p2/(p1q1*s)-192*a2*p1p2**2/(p1q1*s)-192*a1*p1q1/s-
997  &144*a2*p1q1/s-384*a1*a2*mb**2*p1q1/s-480*a2**2*mb**2*p1q1/s+
998  &480*a12*mb*mt*p1q1/s-96*a1*a2*mb*mt*p1q1/s-
999  &864*a12*p1p2*p1q1/s-672*a1*a2*p1p2*p1q1/s-96*a1*a2*p1q1**2/s+
1000  &96*a1*mb*mt**3/(p1q2*s)+192*a2*mb*mt*p1p2/(p1q2*s)-
1001  &96*a1*mt**2*p1p2/(p1q2*s)-192*a2*p1p2**2/(p1q2*s)+
1002  &48*a1*mb*mt*p1q1/(p1q2*s)-96*a2*mb*mt*p1q1/(p1q2*s)-
1003  &48*a1*mt**2*p1q1/(p1q2*s)-192*a1*p1p2*p1q1/(p1q2*s)-
1004  &192*a2*p1p2*p1q1/(p1q2*s)-192*a1*a2*mb*mt*p1p2*p1q1/(p1q2*s)
1005  a18=a18+192*a1*a2*p1p2**2*p1q1/(p1q2*s)-192*a1*p1q1**2/(p1q2*s)-
1006  &192*a2*p1q1**2/(p1q2*s)+192*a1*a2*mb**2*p1q1**2/(p1q2*s)+
1007  &192*a12*mb*mt*p1q1**2/(p1q2*s)-96*a1*a2*mb*mt*p1q1**2/(p1q2*s)+
1008  &192*a1*a2*p1p2*p1q1**2/(p1q2*s)-192*a1*p1q2/s-144*a2*p1q2/s-
1009  &384*a1*a2*mb**2*p1q2/s-480*a2**2*mb**2*p1q2/s+
1010  &480*a12*mb*mt*p1q2/s-96*a1*a2*mb*mt*p1q2/s-
1011  &864*a12*p1p2*p1q2/s-672*a1*a2*p1p2*p1q2/s+
1012  &48*a1*mb*mt*p1q2/(p1q1*s)-96*a2*mb*mt*p1q2/(p1q1*s)-
1013  &48*a1*mt**2*p1q2/(p1q1*s)-192*a1*p1p2*p1q2/(p1q1*s)-
1014  &192*a2*p1p2*p1q2/(p1q1*s)-192*a1*a2*mb*mt*p1p2*p1q2/(p1q1*s)+
1015  &192*a1*a2*p1p2**2*p1q2/(p1q1*s)-576*a1*a2*p1q1*p1q2/s-
1016  &96*a1*a2*p1q2**2/s-192*a1*p1q2**2/(p1q1*s)-
1017  &192*a2*p1q2**2/(p1q1*s)+192*a1*a2*mb**2*p1q2**2/(p1q1*s)+
1018  &192*a12*mb*mt*p1q2**2/(p1q1*s)-96*a1*a2*mb*mt*p1q2**2/(p1q1*s)+
1019  &192*a1*a2*p1p2*p1q2**2/(p1q1*s)-96*a2*mb**3*mt/(p2q1*s)+
1020  &96*a2*mb**2*p1p2/(p2q1*s)-192*a1*mb*mt*p1p2/(p2q1*s)+
1021  &192*a1*p1p2**2/(p2q1*s)+96*a1*mb**2*p1q1/(p2q1*s)
1022  a18=a18+192*a2*mb**2*p1q1/(p2q1*s)-96*a1*mb*mt*p1q1/(p2q1*s)-
1023  &192*a1*a2*mb**3*mt*p1q1/(p2q1*s)+192*a1*p1p2*p1q1/(p2q1*s)+
1024  &192*a1*a2*mb**2*p1p2*p1q1/(p2q1*s)+
1025  &96*a1*a2*mb**2*p1q1**2/(p2q1*s)-
1026  &192*a2*mb**3*mt*p1q1/(p1q2*p2q1*s)+
1027  &192*a2*mb**2*p1p2*p1q1/(p1q2*p2q1*s)-
1028  &96*a1*mb*mt*p1p2*p1q1/(p1q2*p2q1*s)+
1029  &96*a1*p1p2**2*p1q1/(p1q2*p2q1*s)+
1030  &96*a1*mb**2*p1q1**2/(p1q2*p2q1*s)+
1031  &192*a2*mb**2*p1q1**2/(p1q2*p2q1*s)-
1032  &48*a1*mb*mt*p1q1**2/(p1q2*p2q1*s)+
1033  &96*a1*p1p2*p1q1**2/(p1q2*p2q1*s)+96*a1*mb**2*p1q2/(p2q1*s)+
1034  &48*a2*mb**2*p1q2/(p2q1*s)+192*a1*a2*mb**3*mt*p1q2/(p2q1*s)-
1035  &192*a1*a2*mb**2*p1p2*p1q2/(p2q1*s)-
1036  &96*a1*a2*mb**2*p1q2**2/(p2q1*s)+144*a1*p2q1/s+192*a2*p2q1/s+
1037  &96*a1*a2*mb*mt*p2q1/s-480*a2**2*mb*mt*p2q1/s+
1038  &480*a12*mt**2*p2q1/s+384*a1*a2*mt**2*p2q1/s
1039  a18=a18+672*a1*a2*p1p2*p2q1/s+864*a2**2*p1p2*p2q1/s-
1040  &96*a2*mb*mt*p2q1/(p1q1*s)+192*a1*mt**2*p2q1/(p1q1*s)+
1041  &96*a2*mt**2*p2q1/(p1q1*s)-192*a1*a2*mb*mt**3*p2q1/(p1q1*s)+
1042  &192*a2*p1p2*p2q1/(p1q1*s)+192*a1*a2*mt**2*p1p2*p2q1/(p1q1*s)-
1043  &192*a12*p1q1*p2q1/s-192*a2**2*p1q1*p2q1/s+
1044  &48*a1*mt**2*p2q1/(p1q2*s)+96*a2*mt**2*p2q1/(p1q2*s)+
1045  &192*a1*a2*mb*mt**3*p2q1/(p1q2*s)-
1046  &192*a1*a2*mt**2*p1p2*p2q1/(p1q2*s)+
1047  &96*a1*a2*mb*mt*p1q1*p2q1/(p1q2*s)-
1048  &192*a12*mt**2*p1q1*p2q1/(p1q2*s)-
1049  &96*a1*a2*mt**2*p1q1*p2q1/(p1q2*s)-
1050  &384*a1*a2*p1p2*p1q1*p2q1/(p1q2*s)-384*a12*p1q1**2*p2q1/(p1q2*s)-
1051  &384*a1*a2*p1q1**2*p2q1/(p1q2*s)-480*a12*p1q2*p2q1/s-
1052  &960*a1*a2*p1q2*p2q1/s-480*a2**2*p1q2*p2q1/s+
1053  &144*a1*p1q2*p2q1/(p1q1*s)+96*a2*p1q2*p2q1/(p1q1*s)+
1054  &384*a1*a2*mb*mt*p1q2*p2q1/(p1q1*s)-
1055  &96*a12*mt**2*p1q2*p2q1/(p1q1*s)
1056  a18=a18+96*a1*a2*mt**2*p1q2*p2q1/(p1q1*s)-
1057  &576*a1*a2*p1p2*p1q2*p2q1/(p1q1*s)-192*a12*p1q2**2*p2q1/(p1q1*s)-
1058  &384*a1*a2*p1q2**2*p2q1/(p1q1*s)-96*a1*a2*p2q1**2/s-
1059  &96*a1*a2*mt**2*p2q1**2/(p1q1*s)+96*a1*a2*mt**2*p2q1**2/(p1q2*s)+
1060  &288*a1*a2*p1q2*p2q1**2/(p1q1*s)-96*a2*mb**3*mt/(p2q2*s)+
1061  &96*a2*mb**2*p1p2/(p2q2*s)-192*a1*mb*mt*p1p2/(p2q2*s)+
1062  &192*a1*p1p2**2/(p2q2*s)+96*a1*mb**2*p1q1/(p2q2*s)+
1063  &48*a2*mb**2*p1q1/(p2q2*s)+192*a1*a2*mb**3*mt*p1q1/(p2q2*s)-
1064  &192*a1*a2*mb**2*p1p2*p1q1/(p2q2*s)-
1065  &96*a1*a2*mb**2*p1q1**2/(p2q2*s)+96*a1*mb**2*p1q2/(p2q2*s)+
1066  &192*a2*mb**2*p1q2/(p2q2*s)-96*a1*mb*mt*p1q2/(p2q2*s)-
1067  &192*a1*a2*mb**3*mt*p1q2/(p2q2*s)+192*a1*p1p2*p1q2/(p2q2*s)+
1068  &192*a1*a2*mb**2*p1p2*p1q2/(p2q2*s)-
1069  &192*a2*mb**3*mt*p1q2/(p1q1*p2q2*s)+
1070  &192*a2*mb**2*p1p2*p1q2/(p1q1*p2q2*s)-
1071  &96*a1*mb*mt*p1p2*p1q2/(p1q1*p2q2*s)+
1072  &96*a1*p1p2**2*p1q2/(p1q1*p2q2*s)+96*a1*a2*mb**2*p1q2**2/(p2q2*s)
1073  a18=a18+96*a1*mb**2*p1q2**2/(p1q1*p2q2*s)+
1074  &192*a2*mb**2*p1q2**2/(p1q1*p2q2*s)-
1075  &48*a1*mb*mt*p1q2**2/(p1q1*p2q2*s)+
1076  &96*a1*p1p2*p1q2**2/(p1q1*p2q2*s)-48*a2*mb**2*p2q1/(p2q2*s)-
1077  &96*a1*mb*mt*p2q1/(p2q2*s)+48*a2*mb*mt*p2q1/(p2q2*s)-
1078  &192*a1*p1p2*p2q1/(p2q2*s)-192*a2*p1p2*p2q1/(p2q2*s)-
1079  &192*a1*a2*mb*mt*p1p2*p2q1/(p2q2*s)+
1080  &192*a1*a2*p1p2**2*p2q1/(p2q2*s)+
1081  &192*a1*mb*mt**3*p2q1/(p1q1*p2q2*s)+
1082  &96*a2*mb*mt*p1p2*p2q1/(p1q1*p2q2*s)-
1083  &192*a1*mt**2*p1p2*p2q1/(p1q1*p2q2*s)-
1084  &96*a2*p1p2**2*p2q1/(p1q1*p2q2*s)+
1085  &96*a1*a2*mb**2*p1q1*p2q1/(p2q2*s)+
1086  &192*a2**2*mb**2*p1q1*p2q1/(p2q2*s)-
1087  &96*a1*a2*mb*mt*p1q1*p2q1/(p2q2*s)+
1088  &384*a1*a2*p1p2*p1q1*p2q1/(p2q2*s)-96*a1*p1q2*p2q1/(p2q2*s)-
1089  &144*a2*p1q2*p2q1/(p2q2*s)-96*a1*a2*mb**2*p1q2*p2q1/(p2q2*s)
1090  a18=a18+96*a2**2*mb**2*p1q2*p2q1/(p2q2*s)-
1091  &384*a1*a2*mb*mt*p1q2*p2q1/(p2q2*s)+
1092  &576*a1*a2*p1p2*p1q2*p2q1/(p2q2*s)-
1093  &96*a2*mb**2*p1q2*p2q1/(p1q1*p2q2*s)-
1094  &48*a1*mb*mt*p1q2*p2q1/(p1q1*p2q2*s)-
1095  &48*a2*mb*mt*p1q2*p2q1/(p1q1*p2q2*s)-
1096  &96*a1*mt**2*p1q2*p2q1/(p1q1*p2q2*s)-
1097  &96*a1*p1p2*p1q2*p2q1/(p1q1*p2q2*s)-
1098  &96*a2*p1p2*p1q2*p2q1/(p1q1*p2q2*s)+
1099  &96*a1*a2*p1q1*p1q2*p2q1/(p2q2*s)+288*a1*a2*p1q2**2*p2q1/(p2q2*s)-
1100  &96*a1*p1q2**2*p2q1/(p1q1*p2q2*s)-96*a2*p1q2**2*p2q1/(p1q1*p2q2*s)+
1101  &192*a1*p2q1**2/(p2q2*s)+192*a2*p2q1**2/(p2q2*s)+
1102  &96*a1*a2*mb*mt*p2q1**2/(p2q2*s)-192*a2**2*mb*mt*p2q1**2/(p2q2*s)-
1103  &192*a1*a2*mt**2*p2q1**2/(p2q2*s)-192*a1*a2*p1p2*p2q1**2/(p2q2*s)-
1104  &48*a2*mb*mt*p2q1**2/(p1q1*p2q2*s)+
1105  &192*a1*mt**2*p2q1**2/(p1q1*p2q2*s)+
1106  &96*a2*mt**2*p2q1**2/(p1q1*p2q2*s)
1107  a18=a18+96*a2*p1p2*p2q1**2/(p1q1*p2q2*s)-
1108  &384*a1*a2*p1q1*p2q1**2/(p2q2*s)-
1109  &384*a2**2*p1q1*p2q1**2/(p2q2*s)-384*a1*a2*p1q2*p2q1**2/(p2q2*s)-
1110  &192*a2**2*p1q2*p2q1**2/(p2q2*s)+96*a1*p1q2*p2q1**2/(p1q1*p2q2*s)+
1111  &96*a2*p1q2*p2q1**2/(p1q1*p2q2*s)+144*a1*p2q2/s+192*a2*p2q2/s+
1112  &96*a1*a2*mb*mt*p2q2/s-480*a2**2*mb*mt*p2q2/s+
1113  &480*a12*mt**2*p2q2/s+384*a1*a2*mt**2*p2q2/s+
1114  &672*a1*a2*p1p2*p2q2/s+864*a2**2*p1p2*p2q2/s+
1115  &48*a1*mt**2*p2q2/(p1q1*s)+96*a2*mt**2*p2q2/(p1q1*s)+
1116  &192*a1*a2*mb*mt**3*p2q2/(p1q1*s)-
1117  &192*a1*a2*mt**2*p1p2*p2q2/(p1q1*s)-480*a12*p1q1*p2q2/s-
1118  &960*a1*a2*p1q1*p2q2/s-480*a2**2*p1q1*p2q2/s-
1119  &96*a2*mb*mt*p2q2/(p1q2*s)+192*a1*mt**2*p2q2/(p1q2*s)+
1120  &96*a2*mt**2*p2q2/(p1q2*s)-192*a1*a2*mb*mt**3*p2q2/(p1q2*s)+
1121  &192*a2*p1p2*p2q2/(p1q2*s)+192*a1*a2*mt**2*p1p2*p2q2/(p1q2*s)+
1122  &144*a1*p1q1*p2q2/(p1q2*s)+96*a2*p1q1*p2q2/(p1q2*s)+
1123  &384*a1*a2*mb*mt*p1q1*p2q2/(p1q2*s)
1124  a18=a18-96*a12*mt**2*p1q1*p2q2/(p1q2*s)+
1125  &96*a1*a2*mt**2*p1q1*p2q2/(p1q2*s)-
1126  &576*a1*a2*p1p2*p1q1*p2q2/(p1q2*s)-192*a12*p1q1**2*p2q2/(p1q2*s)-
1127  &384*a1*a2*p1q1**2*p2q2/(p1q2*s)-192*a12*p1q2*p2q2/s-
1128  &192*a2**2*p1q2*p2q2/s+96*a1*a2*mb*mt*p1q2*p2q2/(p1q1*s)-
1129  &192*a12*mt**2*p1q2*p2q2/(p1q1*s)-
1130  &96*a1*a2*mt**2*p1q2*p2q2/(p1q1*s)-
1131  &384*a1*a2*p1p2*p1q2*p2q2/(p1q1*s)-384*a12*p1q2**2*p2q2/(p1q1*s)-
1132  &384*a1*a2*p1q2**2*p2q2/(p1q1*s)-48*a2*mb**2*p2q2/(p2q1*s)-
1133  &96*a1*mb*mt*p2q2/(p2q1*s)+48*a2*mb*mt*p2q2/(p2q1*s)-
1134  &192*a1*p1p2*p2q2/(p2q1*s)-192*a2*p1p2*p2q2/(p2q1*s)-
1135  &192*a1*a2*mb*mt*p1p2*p2q2/(p2q1*s)+
1136  &192*a1*a2*p1p2**2*p2q2/(p2q1*s)-96*a1*p1q1*p2q2/(p2q1*s)-
1137  &144*a2*p1q1*p2q2/(p2q1*s)-96*a1*a2*mb**2*p1q1*p2q2/(p2q1*s)+
1138  &96*a2**2*mb**2*p1q1*p2q2/(p2q1*s)-
1139  &384*a1*a2*mb*mt*p1q1*p2q2/(p2q1*s)+
1140  &576*a1*a2*p1p2*p1q1*p2q2/(p2q1*s)+288*a1*a2*p1q1**2*p2q2/(p2q1*s)
1141  a18=a18+192*a1*mb*mt**3*p2q2/(p1q2*p2q1*s)+
1142  &96*a2*mb*mt*p1p2*p2q2/(p1q2*p2q1*s)-
1143  &192*a1*mt**2*p1p2*p2q2/(p1q2*p2q1*s)-
1144  &96*a2*p1p2**2*p2q2/(p1q2*p2q1*s)-
1145  &96*a2*mb**2*p1q1*p2q2/(p1q2*p2q1*s)-
1146  &48*a1*mb*mt*p1q1*p2q2/(p1q2*p2q1*s)-
1147  &48*a2*mb*mt*p1q1*p2q2/(p1q2*p2q1*s)-
1148  &96*a1*mt**2*p1q1*p2q2/(p1q2*p2q1*s)-
1149  &96*a1*p1p2*p1q1*p2q2/(p1q2*p2q1*s)-
1150  &96*a2*p1p2*p1q1*p2q2/(p1q2*p2q1*s)-
1151  &96*a1*p1q1**2*p2q2/(p1q2*p2q1*s)-96*a2*p1q1**2*p2q2/(p1q2*p2q1*s)+
1152  &96*a1*a2*mb**2*p1q2*p2q2/(p2q1*s)+
1153  &192*a2**2*mb**2*p1q2*p2q2/(p2q1*s)-
1154  &96*a1*a2*mb*mt*p1q2*p2q2/(p2q1*s)+
1155  &384*a1*a2*p1p2*p1q2*p2q2/(p2q1*s)+
1156  &96*a1*a2*p1q1*p1q2*p2q2/(p2q1*s)-576*a1*a2*p2q1*p2q2/s+
1157  &96*a1*a2*p1q1*p2q1*p2q2/(p1q2*s)+96*a1*a2*p1q2*p2q1*p2q2/(p1q1*s)
1158  a18=a18-96*a1*a2*p2q2**2/s+96*a1*a2*mt**2*p2q2**2/(p1q1*s)-
1159  &96*a1*a2*mt**2*p2q2**2/(p1q2*s)+288*a1*a2*p1q1*p2q2**2/(p1q2*s)+
1160  &192*a1*p2q2**2/(p2q1*s)+192*a2*p2q2**2/(p2q1*s)+
1161  &96*a1*a2*mb*mt*p2q2**2/(p2q1*s)-192*a2**2*mb*mt*p2q2**2/(p2q1*s)-
1162  &192*a1*a2*mt**2*p2q2**2/(p2q1*s)-192*a1*a2*p1p2*p2q2**2/(p2q1*s)-
1163  &384*a1*a2*p1q1*p2q2**2/(p2q1*s)-192*a2**2*p1q1*p2q2**2/(p2q1*s)-
1164  &48*a2*mb*mt*p2q2**2/(p1q2*p2q1*s)+
1165  &192*a1*mt**2*p2q2**2/(p1q2*p2q1*s)+
1166  &96*a2*mt**2*p2q2**2/(p1q2*p2q1*s)+
1167  &96*a2*p1p2*p2q2**2/(p1q2*p2q1*s)+96*a1*p1q1*p2q2**2/(p1q2*p2q1*s)+
1168  &96*a2*p1q1*p2q2**2/(p1q2*p2q1*s)-384*a1*a2*p1q2*p2q2**2/(p2q1*s)-
1169  &384*a2**2*p1q2*p2q2**2/(p2q1*s)+512*a1*a2*s/3-
1170  &128*a1*mt**2*s/(3*p1q1**2)+128*a12*mb*mt**3*s/(3*p1q1**2)-
1171  &152*a1*s/(3*p1q1)-152*a12*mb*mt*s/(3*p1q1)-
1172  &128*a1*a2*mb*mt*s/(3*p1q1)+112*a1*a2*mt**2*s/(3*p1q1)-
1173  &16*a12*p1p2*s/p1q1+152*a1*a2*p1p2*s/(3*p1q1)-
1174  &128*a1*mt**2*s/(3*p1q2**2)+128*a12*mb*mt**3*s/(3*p1q2**2)
1175  a18=a18-152*a1*s/(3*p1q2)-152*a12*mb*mt*s/(3*p1q2)-
1176  &128*a1*a2*mb*mt*s/(3*p1q2)+112*a1*a2*mt**2*s/(3*p1q2)-
1177  &16*a12*p1p2*s/p1q2+152*a1*a2*p1p2*s/(3*p1q2)+
1178  &16*a1*mb*mt*s/(3*p1q1*p1q2)-32*a12*mb*mt**3*s/(3*p1q1*p1q2)-
1179  &16*a1*p1p2*s/(3*p1q1*p1q2)+272*a1*a2*p1q1*s/(3*p1q2)+
1180  &272*a1*a2*p1q2*s/(3*p1q1)-128*a2*mb**2*s/(3*p2q1**2)+
1181  &128*a2**2*mb**3*mt*s/(3*p2q1**2)+
1182  &32*mb**2*mt**2*s/(3*p1q2**2*p2q1**2)+32*mb**2*s/(3*p1q2*p2q1**2)
1183 
1184  a18bis=
1185  &64*a2*mb**3*mt*s/(3*p1q2*p2q1**2)-
1186  &64*a2*mb**2*mt**2*s/(3*p1q2*p2q1**2)-
1187  &128*a2*mb**2*p1p2*s/(3*p1q2*p2q1**2)-
1188  &128*a2*mb**2*p1q1*s/(3*p1q2*p2q1**2)+
1189  &128*a2**2*mb**2*p1q2*s/(3*p2q1**2)+152*a2*s/(3*p2q1)-
1190  &112*a1*a2*mb**2*s/(3*p2q1)+128*a1*a2*mb*mt*s/(3*p2q1)+
1191  &152*a2**2*mb*mt*s/(3*p2q1)-152*a1*a2*p1p2*s/(3*p2q1)+
1192  &16*a2**2*p1p2*s/p2q1-8*a1*a2*mb**3*mt*s/(3*p1q1*p2q1)+
1193  &16*a1*a2*mb**2*mt**2*s/(3*p1q1*p2q1)-
1194  &8*a1*a2*mb*mt**3*s/(3*p1q1*p2q1)-8*a1*p1p2*s/(3*p1q1*p2q1)-
1195  &8*a2*p1p2*s/(3*p1q1*p2q1)+8*a1*a2*mb**2*p1p2*s/(3*p1q1*p2q1)-
1196  &16*a1*a2*mb*mt*p1p2*s/(3*p1q1*p2q1)+
1197  &8*a1*a2*mt**2*p1p2*s/(3*p1q1*p2q1)+
1198  &32*a1*a2*p1p2**2*s/(3*p1q1*p2q1)-32*a2**2*p1q1*s/(3*p2q1)-
1199  &32*mt**2*s/(3*p1q2**2*p2q1)+64*a1*mb**2*mt**2*s/(3*p1q2**2*p2q1)-
1200  &64*a1*mb*mt**3*s/(3*p1q2**2*p2q1)
1201  a18bis=a18bis+128*a1*mt**2*p1p2*s/(3*p1q2**2*p2q1)-
1202  &12*s/(p1q2*p2q1)+
1203  &24*a1*mb**2*s/(p1q2*p2q1)+64*a1*a2*mb**3*mt*s/(3*p1q2*p2q1)+
1204  &24*a2*mt**2*s/(p1q2*p2q1)-128*a1*a2*mb**2*mt**2*s/(3*p1q2*p2q1)+
1205  &64*a1*a2*mb*mt**3*s/(3*p1q2*p2q1)+56*a1*p1p2*s/(3*p1q2*p2q1)+
1206  &56*a2*p1p2*s/(3*p1q2*p2q1)-64*a1*a2*mb**2*p1p2*s/(3*p1q2*p2q1)+
1207  &128*a1*a2*mb*mt*p1p2*s/(3*p1q2*p2q1)-
1208  &64*a1*a2*mt**2*p1p2*s/(3*p1q2*p2q1)-
1209  &256*a1*a2*p1p2**2*s/(3*p1q2*p2q1)+4*p1p2*s/(3*p1q1*p1q2*p2q1)-
1210  &8*a1*mb*mt*p1p2*s/(3*p1q1*p1q2*p2q1)-
1211  &8*a1*mt**2*p1p2*s/(3*p1q1*p1q2*p2q1)+136*a2*p1q1*s/(3*p1q2*p2q1)-
1212  &128*a1*a2*mb**2*p1q1*s/(3*p1q2*p2q1)+
1213  &128*a1*a2*mb*mt*p1q1*s/(3*p1q2*p2q1)-
1214  &256*a1*a2*p1p2*p1q1*s/(3*p1q2*p2q1)-160*a2**2*p1q2*s/(3*p2q1)+
1215  &16*a1*a2*p1p2*p1q2*s/(3*p1q1*p2q1)-32*a12*p2q1*s/(3*p1q1)-
1216  &128*a12*mt**2*p2q1*s/(3*p1q2**2)-160*a12*p2q1*s/(3*p1q2)-
1217  &128*a2*mb**2*s/(3*p2q2**2)+128*a2**2*mb**3*mt*s/(3*p2q2**2)
1218  a18bis=a18bis+32*mb**2*mt**2*s/(3*p1q1**2*p2q2**2)+
1219  &32*mb**2*s/(3*p1q1*p2q2**2)+
1220  &64*a2*mb**3*mt*s/(3*p1q1*p2q2**2)-
1221  &64*a2*mb**2*mt**2*s/(3*p1q1*p2q2**2)-
1222  &128*a2*mb**2*p1p2*s/(3*p1q1*p2q2**2)+
1223  &128*a2**2*mb**2*p1q1*s/(3*p2q2**2)-
1224  &128*a2*mb**2*p1q2*s/(3*p1q1*p2q2**2)+152*a2*s/(3*p2q2)-
1225  &112*a1*a2*mb**2*s/(3*p2q2)+128*a1*a2*mb*mt*s/(3*p2q2)+
1226  &152*a2**2*mb*mt*s/(3*p2q2)-152*a1*a2*p1p2*s/(3*p2q2)+
1227  &16*a2**2*p1p2*s/p2q2-32*mt**2*s/(3*p1q1**2*p2q2)+
1228  &64*a1*mb**2*mt**2*s/(3*p1q1**2*p2q2)-
1229  &64*a1*mb*mt**3*s/(3*p1q1**2*p2q2)+
1230  &128*a1*mt**2*p1p2*s/(3*p1q1**2*p2q2)-12*s/(p1q1*p2q2)+
1231  &24*a1*mb**2*s/(p1q1*p2q2)+64*a1*a2*mb**3*mt*s/(3*p1q1*p2q2)+
1232  &24*a2*mt**2*s/(p1q1*p2q2)-128*a1*a2*mb**2*mt**2*s/(3*p1q1*p2q2)+
1233  &64*a1*a2*mb*mt**3*s/(3*p1q1*p2q2)+56*a1*p1p2*s/(3*p1q1*p2q2)+
1234  &56*a2*p1p2*s/(3*p1q1*p2q2)-64*a1*a2*mb**2*p1p2*s/(3*p1q1*p2q2)
1235  a18bis=a18bis+128*a1*a2*mb*mt*p1p2*s/(3*p1q1*p2q2)-
1236  &64*a1*a2*mt**2*p1p2*s/(3*p1q1*p2q2)-
1237  &256*a1*a2*p1p2**2*s/(3*p1q1*p2q2)-160*a2**2*p1q1*s/(3*p2q2)-
1238  &8*a1*a2*mb**3*mt*s/(3*p1q2*p2q2)+
1239  &16*a1*a2*mb**2*mt**2*s/(3*p1q2*p2q2)-
1240  &8*a1*a2*mb*mt**3*s/(3*p1q2*p2q2)-8*a1*p1p2*s/(3*p1q2*p2q2)-
1241  &8*a2*p1p2*s/(3*p1q2*p2q2)+8*a1*a2*mb**2*p1p2*s/(3*p1q2*p2q2)-
1242  &16*a1*a2*mb*mt*p1p2*s/(3*p1q2*p2q2)+
1243  &8*a1*a2*mt**2*p1p2*s/(3*p1q2*p2q2)+
1244  &32*a1*a2*p1p2**2*s/(3*p1q2*p2q2)+4*p1p2*s/(3*p1q1*p1q2*p2q2)-
1245  &8*a1*mb*mt*p1p2*s/(3*p1q1*p1q2*p2q2)-
1246  &8*a1*mt**2*p1p2*s/(3*p1q1*p1q2*p2q2)+
1247  &16*a1*a2*p1p2*p1q1*s/(3*p1q2*p2q2)-32*a2**2*p1q2*s/(3*p2q2)+
1248  &136*a2*p1q2*s/(3*p1q1*p2q2)-128*a1*a2*mb**2*p1q2*s/(3*p1q1*p2q2)+
1249  &128*a1*a2*mb*mt*p1q2*s/(3*p1q1*p2q2)-
1250  &256*a1*a2*p1p2*p1q2*s/(3*p1q1*p2q2)+16*a2*mb*mt*s/(3*p2q1*p2q2)-
1251  &32*a2**2*mb**3*mt*s/(3*p2q1*p2q2)-16*a2*p1p2*s/(3*p2q1*p2q2)
1252  a18bis=a18bis-4*p1p2*s/(3*p1q1*p2q1*p2q2)+
1253  &8*a2*mb**2*p1p2*s/(3*p1q1*p2q1*p2q2)+
1254  &8*a2*mb*mt*p1p2*s/(3*p1q1*p2q1*p2q2)-4*p1p2*s/(3*p1q2*p2q1*p2q2)+
1255  &8*a2*mb**2*p1p2*s/(3*p1q2*p2q1*p2q2)+
1256  &8*a2*mb*mt*p1p2*s/(3*p1q2*p2q1*p2q2)-
1257  &2*mb**3*mt*s/(3*p1q1*p1q2*p2q1*p2q2)+
1258  &4*mb**2*mt**2*s/(3*p1q1*p1q2*p2q1*p2q2)-
1259  &2*mb*mt**3*s/(3*p1q1*p1q2*p2q1*p2q2)-
1260  &2*mb**2*p1p2*s/(3*p1q1*p1q2*p2q1*p2q2)+
1261  &4*mb*mt*p1p2*s/(3*p1q1*p1q2*p2q1*p2q2)-
1262  &2*mt**2*p1p2*s/(3*p1q1*p1q2*p2q1*p2q2)-
1263  &8*p1p2**2*s/(3*p1q1*p1q2*p2q1*p2q2)+
1264  &8*a2*p1p2*p1q1*s/(3*p1q2*p2q1*p2q2)+
1265  &8*a2*p1p2*p1q2*s/(3*p1q1*p2q1*p2q2)+272*a1*a2*p2q1*s/(3*p2q2)-
1266  &128*a1*mt**2*p2q1*s/(3*p1q1**2*p2q2)-136*a1*p2q1*s/(3*p1q1*p2q2)-
1267  &128*a1*a2*mb*mt*p2q1*s/(3*p1q1*p2q2)+
1268  &128*a1*a2*mt**2*p2q1*s/(3*p1q1*p2q2)
1269  a18bis=a18bis+256*a1*a2*p1p2*p2q1*s/(3*p1q1*p2q2)-
1270  &16*a1*a2*p1p2*p2q1*s/(3*p1q2*p2q2)+
1271  &8*a1*p1p2*p2q1*s/(3*p1q1*p1q2*p2q2)+
1272  &256*a1*a2*p1q2*p2q1*s/(3*p1q1*p2q2)-
1273  &128*a12*mt**2*p2q2*s/(3*p1q1**2)-160*a12*p2q2*s/(3*p1q1)-
1274  &32*a12*p2q2*s/(3*p1q2)+272*a1*a2*p2q2*s/(3*p2q1)-
1275  &16*a1*a2*p1p2*p2q2*s/(3*p1q1*p2q1)-
1276  &128*a1*mt**2*p2q2*s/(3*p1q2**2*p2q1)-136*a1*p2q2*s/(3*p1q2*p2q1)-
1277  &128*a1*a2*mb*mt*p2q2*s/(3*p1q2*p2q1)+
1278  &128*a1*a2*mt**2*p2q2*s/(3*p1q2*p2q1)+
1279  &256*a1*a2*p1p2*p2q2*s/(3*p1q2*p2q1)+
1280  &8*a1*p1p2*p2q2*s/(3*p1q1*p1q2*p2q1)+
1281  &256*a1*a2*p1q1*p2q2*s/(3*p1q2*p2q1)-
1282  &8*a12*mb*mt*s**2/(3*p1q1*p1q2)+16*a12*p1p2*s**2/(3*p1q1*p1q2)-
1283  &8*a1*a2*p1p2*s**2/(3*p1q1*p2q1)+4*a1*p1p2*s**2/(3*p1q1*p1q2*p2q1)-
1284  &8*a1*a2*p1p2*s**2/(3*p1q2*p2q2)+4*a1*p1p2*s**2/(3*p1q1*p1q2*p2q2)-
1285  &8*a2**2*mb*mt*s**2/(3*p2q1*p2q2)+16*a2**2*p1p2*s**2/(3*p2q1*p2q2)
1286  a18bis=a18bis-4*a2*p1p2*s**2/(3*p1q1*p2q1*p2q2)-
1287  &4*a2*p1p2*s**2/(3*p1q2*p2q1*p2q2)+
1288  &2*p1p2*s**2/(3*p1q1*p1q2*p2q1*p2q2)
1289 C
1290  v18=v18+v18bis
1291  a18=a18+a18bis
1292  v910 =-48*a12*mb*mt-48*a2**2*mb*mt-48*a12*p1p2-48*a2**2*p1p2-
1293  &384*a12*mb*mt*p1q1*p1q2/s**2-384*a12*p1p2*p1q1*p1q2/s**2-
1294  &384*a1*a2*mb*mt*p1q2*p2q1/s**2-384*a1*a2*p1p2*p1q2*p2q1/s**2+
1295  &192*a12*p1q1*p1q2*p2q1/s**2+192*a1*a2*p1q1*p1q2*p2q1/s**2-
1296  &192*a12*p1q2**2*p2q1/s**2-192*a1*a2*p1q2**2*p2q1/s**2+
1297  &192*a1*a2*p1q2*p2q1**2/s**2+192*a2**2*p1q2*p2q1**2/s**2-
1298  &384*a1*a2*mb*mt*p1q1*p2q2/s**2-384*a1*a2*p1p2*p1q1*p2q2/s**2-
1299  &192*a12*p1q1**2*p2q2/s**2-192*a1*a2*p1q1**2*p2q2/s**2+
1300  &192*a12*p1q1*p1q2*p2q2/s**2+192*a1*a2*p1q1*p1q2*p2q2/s**2-
1301  &384*a2**2*mb*mt*p2q1*p2q2/s**2-384*a2**2*p1p2*p2q1*p2q2/s**2-
1302  &192*a1*a2*p1q1*p2q1*p2q2/s**2-192*a2**2*p1q1*p2q1*p2q2/s**2-
1303  &192*a1*a2*p1q2*p2q1*p2q2/s**2-192*a2**2*p1q2*p2q1*p2q2/s**2+
1304  &192*a1*a2*p1q1*p2q2**2/s**2+192*a2**2*p1q1*p2q2**2/s**2+
1305  &96*a12*mb*mt*p1q1/s-96*a1*a2*mb*mt*p1q1/s+
1306  &96*a12*p1p2*p1q1/s-96*a1*a2*p1p2*p1q1/s+96*a12*mb*mt*p1q2/s-
1307  &96*a1*a2*mb*mt*p1q2/s+96*a12*p1p2*p1q2/s-96*a1*a2*p1p2*p1q2/s+
1308  &96*a1*a2*mb*mt*p2q1/s-96*a2**2*mb*mt*p2q1/s
1309  v910=v910+96*a1*a2*p1p2*p2q1/s-
1310  &96*a2**2*p1p2*p2q1/s+96*a12*p1q2*p2q1/s+
1311  &192*a1*a2*p1q2*p2q1/s+96*a2**2*p1q2*p2q1/s+
1312  &96*a1*a2*mb*mt*p2q2/s-96*a2**2*mb*mt*p2q2/s+
1313  &96*a1*a2*p1p2*p2q2/s-96*a2**2*p1p2*p2q2/s+96*a12*p1q1*p2q2/s+
1314  &192*a1*a2*p1q1*p2q2/s+96*a2**2*p1q1*p2q2/s
1315 C
1316  a910 = 48*a12*mb*mt+48*a2**2*mb*mt-48*a12*p1p2-48*a2**2*p1p2+
1317  &384*a12*mb*mt*p1q1*p1q2/s**2-384*a12*p1p2*p1q1*p1q2/s**2+
1318  &384*a1*a2*mb*mt*p1q2*p2q1/s**2-384*a1*a2*p1p2*p1q2*p2q1/s**2+
1319  &192*a12*p1q1*p1q2*p2q1/s**2+192*a1*a2*p1q1*p1q2*p2q1/s**2-
1320  &192*a12*p1q2**2*p2q1/s**2-192*a1*a2*p1q2**2*p2q1/s**2+
1321  &192*a1*a2*p1q2*p2q1**2/s**2+192*a2**2*p1q2*p2q1**2/s**2+
1322  &384*a1*a2*mb*mt*p1q1*p2q2/s**2-384*a1*a2*p1p2*p1q1*p2q2/s**2-
1323  &192*a12*p1q1**2*p2q2/s**2-192*a1*a2*p1q1**2*p2q2/s**2+
1324  &192*a12*p1q1*p1q2*p2q2/s**2+192*a1*a2*p1q1*p1q2*p2q2/s**2+
1325  &384*a2**2*mb*mt*p2q1*p2q2/s**2-384*a2**2*p1p2*p2q1*p2q2/s**2-
1326  &192*a1*a2*p1q1*p2q1*p2q2/s**2-192*a2**2*p1q1*p2q1*p2q2/s**2-
1327  &192*a1*a2*p1q2*p2q1*p2q2/s**2-192*a2**2*p1q2*p2q1*p2q2/s**2+
1328  &192*a1*a2*p1q1*p2q2**2/s**2+192*a2**2*p1q1*p2q2**2/s**2-
1329  &96*a12*mb*mt*p1q1/s+96*a1*a2*mb*mt*p1q1/s+
1330  &96*a12*p1p2*p1q1/s-96*a1*a2*p1p2*p1q1/s-96*a12*mb*mt*p1q2/s+
1331  &96*a1*a2*mb*mt*p1q2/s+96*a12*p1p2*p1q2/s-96*a1*a2*p1p2*p1q2/s-
1332  &96*a1*a2*mb*mt*p2q1/s+96*a2**2*mb*mt*p2q1/s
1333  a910=a910+96*a1*a2*p1p2*p2q1/s-
1334  &96*a2**2*p1p2*p2q1/s+96*a12*p1q2*p2q1/s+
1335  &192*a1*a2*p1q2*p2q1/s+96*a2**2*p1q2*p2q1/s-
1336  &96*a1*a2*mb*mt*p2q2/s+96*a2**2*mb*mt*p2q2/s+
1337  &96*a1*a2*p1p2*p2q2/s-96*a2**2*p1p2*p2q2/s+96*a12*p1q1*p2q2/s+
1338  &192*a1*a2*p1q1*p2q2/s+96*a2**2*p1q1*p2q2/s
1339 C
1340 C FINAL RESULT;
1341 C
1342  amp2= fact*ps*vtb**2*(v**2 *(v18 +v910)+a**2 *(a18+a910) )
1343 
1344  END