Back to home page

sPhenix code displayed by LXR
 
 

    

File indexing completed on 2025-08-05 08:21:19 

0001  
0002 C AMPLITUDE SQUARED (MATRIX ELEMENTS) FOR THE PROCESSES:
0003 C GG->TBH^+, QQBAR->TBH^+
0004 C AS A FUNCTION OF 4-MOMENTA FOR SUITABLE INTERFACE
0005 C (FOR INSTANCE WITH PYTHIA)
0006 C------------------------------------------------------------
0007 C BASED ON F. BORZUMATI, J.-L. KNEUR, N. POLONSKY  HEP-PH/9905443,
0008 C PHYS REV. D 60 (1999) 115011
0009 C (THESE FILES PREPARED BY J.-L. KNEUR)
0010 C------------------------------------------------------------
0011 C 1)  GG->TBH^+
0012        SUBROUTINE PYTBHG(Q1,Q2,P1,P2,P3,MT,MB,RMB,MHP,AMP2)
0013 C
0014 C CONVENTIONS AND INPUT/OUTPUT DEFINITIONS:
0015 C
0016 C INPUT: Q1,Q2 ARE ENTERING 4-MOMENTA OF INITIAL GLUONS OR QUARKS;
0017 C        P1, P2 ARE THE TOP AND BOTTOM OUTGOING 4-MOMENTA;
0018 C        P3 IS OUTGOING CHARGED HIGGS 4-MOMENTA.
0019 C  (NB FOR ALL 4-MOMENTA P(4) IS TIME-COMPONENT)
0020 C "PHYSICAL PARAMETERS" INPUT:
0021 C        MT,MB TOP AND BOTTOM MASSES;
0022 C        MHP CHARGED HIGGS MASS
0023 C   FURTHER PARAMETERS INPUT IS NEEDED FROM COMMON/PARAM/ (SEE BELOW)
0024 C
0025 C OUTPUT: AMP2  IS MATRIX ELEMENT (AMPLITUDE**2) FOR GG->TB H^+
0026 C (NB AMP2 IS TRULY AMPLITUDE SQUARRED, I.E. WITHOUT ANY
0027 C PHASE SPACE FACTORS INCLUDED. IT INCLUDES COLOUR AND COUPLING
0028 C FACTORS, AS EXPLICIT BELOW. ACCORDINGLY, FOR EXAMPLE THE TOTAL
0029 C CROSS-SECTION SHOULD BE (SYMBOLICALLY):
0030 C   SIGMA = INTEGRATE [PARTON DENSITY FUNCTIONS * 3-PARTICLE FINAL
0031 C           STATE PHASE-SPACE (STANDARDLY NORMALIZED) * AMP2 ]
0032 C
0033       IMPLICIT DOUBLE PRECISION(A-H, O-Z)
0034       IMPLICIT INTEGER(I-N)
0035       DOUBLE PRECISION MW2,MT,MB,MHP,MW
0036       DIMENSION Q1(4),Q2(4),P1(4),P2(4),P3(4)
0037       COMMON/PYDAT1/MSTU(200),PARU(200),MSTJ(200),PARJ(200)
0038       COMMON/PYDAT2/KCHG(500,4),PMAS(500,4),PARF(2000),VCKM(4,4)
0039       COMMON/PYMSSM/IMSS(0:99),RMSS(0:99)
0040  
0041       COMMON/PYCTBH/ ALPHA,ALPHAS,SW2,MW2,TANB,VTB,V,A
0042       SAVE /PYDAT1/,/PYDAT2/,/PYMSSM/,/PYCTBH/
0043 C !THE RELEVANT INPUT PARAMETERS ABOVE ARE NEEDED FOR CALCULATION
0044 C BUT ARE NOT DEFINED HERE SO THAT ONE MAY CHOOSE/VARY THEIR VALUES:
0045 C ACCORDINGLY, WHEN CALLING THESE SUBROUTINES, PLEASE SUPPLY VIA
0046 C THIS COMMON/PARAM/ YOUR PREFERRED ALPHA, ALPHAS,..AND TANB
0047 C (TAN BETA) VALUES
0048 C
0049 C THE NORMALIZED V,A COUPLINGS ARE DEFINED BELOW AND USED BOTH
0050 C IN THIS ROUTINE AND IN THE TOP WIDTH CALCULATION PYTBHB(..).
0051  
0052       PI = 4*DATAN(1.D0)
0053       MW = DSQRT(MW2)
0054 C
0055 C COLLECTING THE RELEVANT OVERALL FACTORS:
0056 C 8X8 INITIAL GLUON COLOR AVERAGE, 2X2 GLUON SPIN AVERAGE
0057       PS=1.D0/(8.D0*8.D0 *2.D0*2.D0)
0058 C COUPLING CONSTANT (OVERALL NORMALIZATION)
0059       FACT=(4.D0*PI*ALPHA)*(4.D0*PI*ALPHAS)**2/SW2/2.D0
0060 C NB ALPHA IS E^2/4/PI, BUT BETTER DEFINED IN TERMS OF G_FERMI:
0061 C ALPHA= DSQRT(2.D0)*GF*SW2*MW**2/PI
0062 C ALPHAS IS ALPHA_STRONG;
0063 C SW2 IS SIN(THETA_W)**2.
0064 C
0065 C      VTB=.998D0
0066 C VTB IS TOP-BOTTOM CKM MATRIX ELEMENT (APPROXIMATE VALUE HERE)
0067 C
0068       V = ( MT/MW/TANB +RMB/MW*TANB)/2.D0
0069       A = (-MT/MW/TANB +RMB/MW*TANB)/2.D0
0070 C V AND A ARE (NORMALIZED) VECTOR AND AXIAL TBH^+ COUPLINGS
0071 C
0072 C REDEFINING P2 INGOING FROM OVERALL MOMENTUM CONSERVATION
0073 C (BECAUSE P2 INGOING WAS USED IN OUR GRAPH CALCULATION CONVENTIONS)
0074       DO 100 KK=1,4
0075       P2(KK)=P3(KK)-Q1(KK)-Q2(KK)+P1(KK)
0076   100 CONTINUE
0077 C DEFINING VARIOUS RELEVANT 4-SCALAR PRODUCTS:
0078       S = 2*PYTBHS(Q1,Q2)
0079       P1Q1=PYTBHS(Q1,P1)
0080       P1Q2=PYTBHS(P1,Q2)
0081       P2Q1=PYTBHS(P2,Q1)
0082       P2Q2=PYTBHS(P2,Q2)
0083       P1P2=PYTBHS(P1,P2)
0084 C
0085 C   TOP WIDTH CALCULATION
0086       CALL PYTBHB(MT,MB,MHP,BR,GAMT)
0087 C   GAMT IS THE TOP WIDTH: T->BH^+ AND/OR T->B W^+
0088 C THEN DEFINE TOP (RESONANT) PROPAGATOR:
0089       A1INV= S -2*P1Q1 -2*P1Q2
0090       A1 =A1INV/(A1INV**2+ (GAMT*MT)**2)
0091 C (I.E. INTRODUCE THE TOP WIDTH IN A1 TO REGULARISE THE POLE)
0092 C  NB:    A12 = A1*A1 BUT CORRECT EXPRESSION BELOW BECAUSE OF
0093 C  THE TOP WIDTH
0094       A12 = 1.D0/(A1INV**2+ (GAMT*MT)**2)
0095       A2 =1.D0/(S +2*P2Q1 +2*P2Q2)
0096 C NOTE A2 IS B PROPAGATOR, DOES NOT NEED A WIDTH
0097 C  NOW COMES THE AMP**2:
0098 C NB COLOR FACTOR (COMING FROM GRAPHS) ALREADY INCLUDED IN
0099 C THE EXPRESSIONS BELOW
0100       V18=0.D0
0101       A18=0.D0
0102       V18= 640*A1/3+640*A2/3+32*A1*A2*MB**2-368*A12*MB*MT-
0103      &512*A1*A2*MB*MT/3-
0104      &368*A2**2*MB*MT+32*A1*A2*MT**2+496*A12*P1P2/3+
0105      &320*A1*A2*P1P2+496*A2**2*P1P2/3+128*A1*MB*MT**3/(3*P1Q1**2)+
0106      &128*A1*MT**4/(3*P1Q1**2)-256*A12*MB*MT**5/(3*P1Q1**2)+
0107      &256*A1*MT**2*P1P2/(3*P1Q1**2)-256*A12*MT**4*P1P2/(3*P1Q1**2)+
0108      &8/(3*P1Q1)-32*A1*MB*MT/P1Q1-56*A2*MB*MT/(3*P1Q1)+
0109      &88*A1*MT**2/(3*P1Q1)+72*A2*MT**2/P1Q1+
0110      &704*A12*MB*MT**3/(3*P1Q1)-224*A1*A2*MB*MT**3/(3*P1Q1)+
0111      &104*A1*P1P2/(3*P1Q1)+48*A2*P1P2/P1Q1+
0112      &128*A1*A2*MB*MT*P1P2/(3*P1Q1)+512*A12*MT**2*P1P2/(3*P1Q1)-
0113      &448*A1*A2*MT**2*P1P2/(3*P1Q1)-32*A1*A2*P1P2**2/P1Q1-
0114      &656*A1*A2*P1Q1/3-224*A2**2*P1Q1+128*A1*MB*MT**3/(3*P1Q2**2)+
0115      &128*A1*MT**4/(3*P1Q2**2)-256*A12*MB*MT**5/(3*P1Q2**2)+
0116      &256*A1*MT**2*P1P2/(3*P1Q2**2)-256*A12*MT**4*P1P2/(3*P1Q2**2)+
0117      &256*A1*MT**2*P1Q1/(3*P1Q2**2)+256*A12*MB*MT**3*P1Q1/(3*P1Q2**2)+
0118      &8/(3*P1Q2)-32*A1*MB*MT/P1Q2-56*A2*MB*MT/(3*P1Q2)
0119       V18=V18+88*A1*MT**2/(3*P1Q2)+72*A2*MT**2/P1Q2+
0120      &704*A12*MB*MT**3/(3*P1Q2)-224*A1*A2*MB*MT**3/(3*P1Q2)+
0121      &104*A1*P1P2/(3*P1Q2)+48*A2*P1P2/P1Q2+
0122      &128*A1*A2*MB*MT*P1P2/(3*P1Q2)+512*A12*MT**2*P1P2/(3*P1Q2)-
0123      &448*A1*A2*MT**2*P1P2/(3*P1Q2)-32*A1*A2*P1P2**2/P1Q2-
0124      &32*A1*MB*MT**3/(3*P1Q1*P1Q2)-32*A1*MT**4/(3*P1Q1*P1Q2)+
0125      &64*A12*MB*MT**5/(3*P1Q1*P1Q2)+16*P1P2/(3*P1Q1*P1Q2)-
0126      &64*A1*MT**2*P1P2/(3*P1Q1*P1Q2)+64*A12*MT**4*P1P2/(3*P1Q1*P1Q2)+
0127      &112*A1*P1Q1/P1Q2+272*A2*P1Q1/(3*P1Q2)-
0128      &272*A1*A2*MB**2*P1Q1/(3*P1Q2)+208*A12*MB*MT*P1Q1/(3*P1Q2)-
0129      &400*A1*A2*MB*MT*P1Q1/(3*P1Q2)-80*A1*A2*MT**2*P1Q1/P1Q2+
0130      &96*A12*P1P2*P1Q1/P1Q2-320*A1*A2*P1P2*P1Q1/P1Q2-
0131      &544*A1*A2*P1Q1**2/(3*P1Q2)-656*A1*A2*P1Q2/3-224*A2**2*P1Q2+
0132      &256*A1*MT**2*P1Q2/(3*P1Q1**2)+256*A12*MB*MT**3*P1Q2/(3*P1Q1**2)+
0133      &112*A1*P1Q2/P1Q1+272*A2*P1Q2/(3*P1Q1)-
0134      &272*A1*A2*MB**2*P1Q2/(3*P1Q1)+208*A12*MB*MT*P1Q2/(3*P1Q1)-
0135      &400*A1*A2*MB*MT*P1Q2/(3*P1Q1)-80*A1*A2*MT**2*P1Q2/P1Q1
0136       V18=V18+96*A12*P1P2*P1Q2/P1Q1-320*A1*A2*P1P2*P1Q2/P1Q1-
0137      &544*A1*A2*P1Q2**2/(3*P1Q1)+128*A2*MB**4/(3*P2Q1**2)+
0138      &128*A2*MB**3*MT/(3*P2Q1**2)-256*A2**2*MB**5*MT/(3*P2Q1**2)+
0139      &256*A2*MB**2*P1P2/(3*P2Q1**2)-256*A2**2*MB**4*P1P2/(3*P2Q1**2)+
0140      &256*A2*MB**2*P1Q1/(3*P2Q1**2)-256*A2**2*MB**4*P1Q1/(3*P2Q1**2)-
0141      &64*MB**3*MT**3/(3*P1Q2**2*P2Q1**2)-
0142      &64*MB**2*MT**2*P1P2/(3*P1Q2**2*P2Q1**2)-
0143      &64*MB**2*MT**2*P1Q1/(3*P1Q2**2*P2Q1**2)+
0144      &64*MB**3*MT/(3*P1Q2*P2Q1**2)+
0145      &256*A2*MB**3*MT*P1P2/(3*P1Q2*P2Q1**2)+
0146      &256*A2*MB**2*P1P2**2/(3*P1Q2*P2Q1**2)+
0147      &256*A2*MB**3*MT*P1Q1/(3*P1Q2*P2Q1**2)+
0148      &512*A2*MB**2*P1P2*P1Q1/(3*P1Q2*P2Q1**2)+
0149      &256*A2*MB**2*P1Q1**2/(3*P1Q2*P2Q1**2)-
0150      &256*A2**2*MB**4*P1Q2/(3*P2Q1**2)-8/(3*P2Q1)-72*A1*MB**2/P2Q1-
0151      &88*A2*MB**2/(3*P2Q1)+56*A1*MB*MT/(3*P2Q1)+32*A2*MB*MT/P2Q1+
0152      &224*A1*A2*MB**3*MT/(3*P2Q1)-704*A2**2*MB**3*MT/(3*P2Q1)
0153       V18=V18-48*A1*P1P2/P2Q1-104*A2*P1P2/(3*P2Q1)+
0154      &448*A1*A2*MB**2*P1P2/(3*P2Q1)-512*A2**2*MB**2*P1P2/(3*P2Q1)-
0155      &128*A1*A2*MB*MT*P1P2/(3*P2Q1)+32*A1*A2*P1P2**2/P2Q1-
0156      &16*P1P2/(3*P1Q1*P2Q1)-32*A1*MB*MT*P1P2/(3*P1Q1*P2Q1)-
0157      &32*A2*MB*MT*P1P2/(3*P1Q1*P2Q1)-
0158      &64*A1*A2*MB*MT*P1P2**2/(3*P1Q1*P2Q1)-
0159      &64*A1*A2*P1P2**3/(3*P1Q1*P2Q1)-256*A2*P1Q1/(3*P2Q1)+
0160      &448*A1*A2*MB**2*P1Q1/(3*P2Q1)-368*A2**2*MB**2*P1Q1/(3*P2Q1)+
0161      &224*A1*A2*MB*MT*P1Q1/(3*P2Q1)+304*A1*A2*P1P2*P1Q1/(3*P2Q1)-
0162      &64*MB*MT**3/(3*P1Q2**2*P2Q1)-
0163      &256*A1*MB*MT**3*P1P2/(3*P1Q2**2*P2Q1)-
0164      &256*A1*MT**2*P1P2**2/(3*P1Q2**2*P2Q1)+
0165      &64*MT**2*P1Q1/(3*P1Q2**2*P2Q1)-
0166      &128*A1*MB**2*MT**2*P1Q1/(3*P1Q2**2*P2Q1)-
0167      &128*A1*MB*MT**3*P1Q1/(3*P1Q2**2*P2Q1)-
0168      &256*A1*MT**2*P1P2*P1Q1/(3*P1Q2**2*P2Q1)-4*MB**2/(3*P1Q2*P2Q1)+
0169      &64*MB*MT/(3*P1Q2*P2Q1)-128*A2*MB**3*MT/(3*P1Q2*P2Q1)
0170       V18=V18-4*MT**2/(3*P1Q2*P2Q1)-128*A1*MB**2*MT**2/(3*P1Q2*P2Q1)-
0171      &128*A2*MB**2*MT**2/(3*P1Q2*P2Q1)-128*A1*MB*MT**3/(3*P1Q2*P2Q1)-
0172      &112*A2*MB**2*P1P2/(3*P1Q2*P2Q1)-32*A1*MB*MT*P1P2/(3*P1Q2*P2Q1)-
0173      &32*A2*MB*MT*P1P2/(3*P1Q2*P2Q1)-112*A1*MT**2*P1P2/(3*P1Q2*P2Q1)-
0174      &48*A1*P1P2**2/(P1Q2*P2Q1)-48*A2*P1P2**2/(P1Q2*P2Q1)+
0175      &512*A1*A2*MB*MT*P1P2**2/(3*P1Q2*P2Q1)+
0176      &512*A1*A2*P1P2**3/(3*P1Q2*P2Q1)-8*MB*MT*P1P2/(3*P1Q1*P1Q2*P2Q1)-
0177      &8*MT**2*P1P2/(3*P1Q1*P1Q2*P2Q1)+
0178      &32*A1*MB*MT**3*P1P2/(3*P1Q1*P1Q2*P2Q1)-
0179      &16*P1P2**2/(3*P1Q1*P1Q2*P2Q1)+
0180      &32*A1*MT**2*P1P2**2/(3*P1Q1*P1Q2*P2Q1)+8*P1Q1/(3*P1Q2*P2Q1)-
0181      &160*A1*MB**2*P1Q1/(3*P1Q2*P2Q1)-272*A2*MB**2*P1Q1/(3*P1Q2*P2Q1)+
0182      &56*A1*MB*MT*P1Q1/(3*P1Q2*P2Q1)+200*A2*MB*MT*P1Q1/(3*P1Q2*P2Q1)-
0183      &48*A1*P1P2*P1Q1/(P1Q2*P2Q1)-256*A2*P1P2*P1Q1/(3*P1Q2*P2Q1)+
0184      &256*A1*A2*MB**2*P1P2*P1Q1/(3*P1Q2*P2Q1)+
0185      &256*A1*A2*MB*MT*P1P2*P1Q1/(P1Q2*P2Q1)+
0186      &1024*A1*A2*P1P2**2*P1Q1/(3*P1Q2*P2Q1)
0187       V18=V18-272*A2*P1Q1**2/(3*P1Q2*P2Q1)+
0188      &256*A1*A2*MB**2*P1Q1**2/(3*P1Q2*P2Q1)+
0189      &256*A1*A2*MB*MT*P1Q1**2/(3*P1Q2*P2Q1)+
0190      &512*A1*A2*P1P2*P1Q1**2/(3*P1Q2*P2Q1)+16*A2*P1Q2/(3*P2Q1)+
0191      &64*A1*A2*MB**2*P1Q2/P2Q1+32*A2**2*MB**2*P1Q2/(3*P2Q1)+
0192      &112*A1*A2*MB*MT*P1Q2/(3*P2Q1)+368*A1*A2*P1P2*P1Q2/(3*P2Q1)+
0193      &32*A2*P1P2*P1Q2/(3*P1Q1*P2Q1)-
0194      &32*A1*A2*MB**2*P1P2*P1Q2/(3*P1Q1*P2Q1)-
0195      &32*A1*A2*MB*MT*P1P2*P1Q2/(3*P1Q1*P2Q1)-
0196      &64*A1*A2*P1P2**2*P1Q2/(3*P1Q1*P2Q1)+224*A12*P2Q1+
0197      &656*A1*A2*P2Q1/3-256*A1*MT**2*P2Q1/(3*P1Q1**2)+
0198      &256*A12*MT**4*P2Q1/(3*P1Q1**2)-256*A1*P2Q1/(3*P1Q1)+
0199      &224*A1*A2*MB*MT*P2Q1/(3*P1Q1)-368*A12*MT**2*P2Q1/(3*P1Q1)+
0200      &448*A1*A2*MT**2*P2Q1/(3*P1Q1)+304*A1*A2*P1P2*P2Q1/(3*P1Q1)+
0201      &256*A12*MT**4*P2Q1/(3*P1Q2**2)+
0202      &256*A12*MT**2*P1Q1*P2Q1/(3*P1Q2**2)+16*A1*P2Q1/(3*P1Q2)+
0203      &112*A1*A2*MB*MT*P2Q1/(3*P1Q2)+32*A12*MT**2*P2Q1/(3*P1Q2)
0204       V18=V18+64*A1*A2*MT**2*P2Q1/P1Q2+368*A1*A2*P1P2*P2Q1/(3*P1Q2)+
0205      &16*A1*MT**2*P2Q1/(3*P1Q1*P1Q2)-64*A12*MT**4*P2Q1/(3*P1Q1*P1Q2)+
0206      &640*A12*P1Q1*P2Q1/(3*P1Q2)+544*A1*A2*P1Q1*P2Q1/(3*P1Q2)+
0207      &32*A12*P1Q2*P2Q1/P1Q1+944*A1*A2*P1Q2*P2Q1/(3*P1Q1)+
0208      &128*A2*MB**4/(3*P2Q2**2)+128*A2*MB**3*MT/(3*P2Q2**2)-
0209      &256*A2**2*MB**5*MT/(3*P2Q2**2)+256*A2*MB**2*P1P2/(3*P2Q2**2)-
0210      &256*A2**2*MB**4*P1P2/(3*P2Q2**2)-
0211      &64*MB**3*MT**3/(3*P1Q1**2*P2Q2**2)-
0212      &64*MB**2*MT**2*P1P2/(3*P1Q1**2*P2Q2**2)+
0213      &64*MB**3*MT/(3*P1Q1*P2Q2**2)+
0214      &256*A2*MB**3*MT*P1P2/(3*P1Q1*P2Q2**2)+
0215      &256*A2*MB**2*P1P2**2/(3*P1Q1*P2Q2**2)-
0216      &256*A2**2*MB**4*P1Q1/(3*P2Q2**2)+256*A2*MB**2*P1Q2/(3*P2Q2**2)-
0217      &256*A2**2*MB**4*P1Q2/(3*P2Q2**2)-
0218      &64*MB**2*MT**2*P1Q2/(3*P1Q1**2*P2Q2**2)+
0219      &256*A2*MB**3*MT*P1Q2/(3*P1Q1*P2Q2**2)+
0220      &512*A2*MB**2*P1P2*P1Q2/(3*P1Q1*P2Q2**2)
0221       V18=V18+256*A2*MB**2*P1Q2**2/(3*P1Q1*P2Q2**2)-
0222      &256*A2*MB**2*P2Q1/(3*P2Q2**2)-256*A2**2*MB**3*MT*P2Q1/(3*P2Q2**2)+
0223      &64*MB**2*MT**2*P2Q1/(3*P1Q1**2*P2Q2**2)+
0224      &64*MB**2*P2Q1/(3*P1Q1*P2Q2**2)-
0225      &128*A2*MB**3*MT*P2Q1/(3*P1Q1*P2Q2**2)-
0226      &128*A2*MB**2*MT**2*P2Q1/(3*P1Q1*P2Q2**2)-
0227      &256*A2*MB**2*P1P2*P2Q1/(3*P1Q1*P2Q2**2)+
0228      &256*A2**2*MB**2*P1Q1*P2Q1/(3*P2Q2**2)-
0229      &256*A2*MB**2*P1Q2*P2Q1/(3*P1Q1*P2Q2**2)-8/(3*P2Q2)-
0230      &72*A1*MB**2/P2Q2-88*A2*MB**2/(3*P2Q2)+56*A1*MB*MT/(3*P2Q2)+
0231      &32*A2*MB*MT/P2Q2+224*A1*A2*MB**3*MT/(3*P2Q2)-
0232      &704*A2**2*MB**3*MT/(3*P2Q2)-48*A1*P1P2/P2Q2-
0233      &104*A2*P1P2/(3*P2Q2)+448*A1*A2*MB**2*P1P2/(3*P2Q2)-
0234      &512*A2**2*MB**2*P1P2/(3*P2Q2)-128*A1*A2*MB*MT*P1P2/(3*P2Q2)+
0235      &32*A1*A2*P1P2**2/P2Q2-64*MB*MT**3/(3*P1Q1**2*P2Q2)-
0236      &256*A1*MB*MT**3*P1P2/(3*P1Q1**2*P2Q2)-
0237      &256*A1*MT**2*P1P2**2/(3*P1Q1**2*P2Q2)-4*MB**2/(3*P1Q1*P2Q2)
0238       V18=V18+64*MB*MT/(3*P1Q1*P2Q2)-128*A2*MB**3*MT/(3*P1Q1*P2Q2)-
0239      &4*MT**2/(3*P1Q1*P2Q2)-128*A1*MB**2*MT**2/(3*P1Q1*P2Q2)-
0240      &128*A2*MB**2*MT**2/(3*P1Q1*P2Q2)-128*A1*MB*MT**3/(3*P1Q1*P2Q2)-
0241      &112*A2*MB**2*P1P2/(3*P1Q1*P2Q2)-32*A1*MB*MT*P1P2/(3*P1Q1*P2Q2)-
0242      &32*A2*MB*MT*P1P2/(3*P1Q1*P2Q2)-112*A1*MT**2*P1P2/(3*P1Q1*P2Q2)-
0243      &48*A1*P1P2**2/(P1Q1*P2Q2)-48*A2*P1P2**2/(P1Q1*P2Q2)+
0244      &512*A1*A2*MB*MT*P1P2**2/(3*P1Q1*P2Q2)+
0245      &512*A1*A2*P1P2**3/(3*P1Q1*P2Q2)+16*A2*P1Q1/(3*P2Q2)+
0246      &64*A1*A2*MB**2*P1Q1/P2Q2+32*A2**2*MB**2*P1Q1/(3*P2Q2)+
0247      &112*A1*A2*MB*MT*P1Q1/(3*P2Q2)+368*A1*A2*P1P2*P1Q1/(3*P2Q2)-
0248      &16*P1P2/(3*P1Q2*P2Q2)-32*A1*MB*MT*P1P2/(3*P1Q2*P2Q2)-
0249      &32*A2*MB*MT*P1P2/(3*P1Q2*P2Q2)-
0250      &64*A1*A2*MB*MT*P1P2**2/(3*P1Q2*P2Q2)-
0251      &64*A1*A2*P1P2**3/(3*P1Q2*P2Q2)-8*MB*MT*P1P2/(3*P1Q1*P1Q2*P2Q2)-
0252      &8*MT**2*P1P2/(3*P1Q1*P1Q2*P2Q2)+
0253      &32*A1*MB*MT**3*P1P2/(3*P1Q1*P1Q2*P2Q2)-
0254      &16*P1P2**2/(3*P1Q1*P1Q2*P2Q2)
0255       V18=V18+32*A1*MT**2*P1P2**2/(3*P1Q1*P1Q2*P2Q2)+
0256      &32*A2*P1P2*P1Q1/(3*P1Q2*P2Q2)-
0257      &32*A1*A2*MB**2*P1P2*P1Q1/(3*P1Q2*P2Q2)-
0258      &32*A1*A2*MB*MT*P1P2*P1Q1/(3*P1Q2*P2Q2)-
0259      &64*A1*A2*P1P2**2*P1Q1/(3*P1Q2*P2Q2)-256*A2*P1Q2/(3*P2Q2)+
0260      &448*A1*A2*MB**2*P1Q2/(3*P2Q2)-368*A2**2*MB**2*P1Q2/(3*P2Q2)+
0261      &224*A1*A2*MB*MT*P1Q2/(3*P2Q2)+304*A1*A2*P1P2*P1Q2/(3*P2Q2)+
0262      &64*MT**2*P1Q2/(3*P1Q1**2*P2Q2)-
0263      &128*A1*MB**2*MT**2*P1Q2/(3*P1Q1**2*P2Q2)-
0264      &128*A1*MB*MT**3*P1Q2/(3*P1Q1**2*P2Q2)-
0265      &256*A1*MT**2*P1P2*P1Q2/(3*P1Q1**2*P2Q2)+8*P1Q2/(3*P1Q1*P2Q2)-
0266      &160*A1*MB**2*P1Q2/(3*P1Q1*P2Q2)-272*A2*MB**2*P1Q2/(3*P1Q1*P2Q2)+
0267      &56*A1*MB*MT*P1Q2/(3*P1Q1*P2Q2)+200*A2*MB*MT*P1Q2/(3*P1Q1*P2Q2)-
0268      &48*A1*P1P2*P1Q2/(P1Q1*P2Q2)-256*A2*P1P2*P1Q2/(3*P1Q1*P2Q2)+
0269      &256*A1*A2*MB**2*P1P2*P1Q2/(3*P1Q1*P2Q2)+
0270      &256*A1*A2*MB*MT*P1P2*P1Q2/(P1Q1*P2Q2)+
0271      &1024*A1*A2*P1P2**2*P1Q2/(3*P1Q1*P2Q2)
0272       V18=V18-272*A2*P1Q2**2/(3*P1Q1*P2Q2)+
0273      &256*A1*A2*MB**2*P1Q2**2/(3*P1Q1*P2Q2)+
0274      &256*A1*A2*MB*MT*P1Q2**2/(3*P1Q1*P2Q2)+
0275      &512*A1*A2*P1P2*P1Q2**2/(3*P1Q1*P2Q2)-32*A2*MB**4/(3*P2Q1*P2Q2)-
0276      &32*A2*MB**3*MT/(3*P2Q1*P2Q2)+64*A2**2*MB**5*MT/(3*P2Q1*P2Q2)+
0277      &16*P1P2/(3*P2Q1*P2Q2)-64*A2*MB**2*P1P2/(3*P2Q1*P2Q2)+
0278      &64*A2**2*MB**4*P1P2/(3*P2Q1*P2Q2)+8*MB**2*P1P2/(3*P1Q1*P2Q1*P2Q2)+
0279      &8*MB*MT*P1P2/(3*P1Q1*P2Q1*P2Q2)-
0280      &32*A2*MB**3*MT*P1P2/(3*P1Q1*P2Q1*P2Q2)+
0281      &16*P1P2**2/(3*P1Q1*P2Q1*P2Q2)-
0282      &32*A2*MB**2*P1P2**2/(3*P1Q1*P2Q1*P2Q2)-
0283      &16*A2*MB**2*P1Q1/(3*P2Q1*P2Q2)+64*A2**2*MB**4*P1Q1/(3*P2Q1*P2Q2)+
0284      &8*MB**2*P1P2/(3*P1Q2*P2Q1*P2Q2)+8*MB*MT*P1P2/(3*P1Q2*P2Q1*P2Q2)-
0285      &32*A2*MB**3*MT*P1P2/(3*P1Q2*P2Q1*P2Q2)+
0286      &16*P1P2**2/(3*P1Q2*P2Q1*P2Q2)-
0287      &32*A2*MB**2*P1P2**2/(3*P1Q2*P2Q1*P2Q2)+
0288      &16*MB*MT*P1P2**2/(3*P1Q1*P1Q2*P2Q1*P2Q2)
0289       V18=V18+16*P1P2**3/(3*P1Q1*P1Q2*P2Q1*P2Q2)-
0290      &32*A2*MB**2*P1P2*P1Q1/(3*P1Q2*P2Q1*P2Q2)-
0291      &16*A2*MB**2*P1Q2/(3*P2Q1*P2Q2)+64*A2**2*MB**4*P1Q2/(3*P2Q1*P2Q2)-
0292      &32*A2*MB**2*P1P2*P1Q2/(3*P1Q1*P2Q1*P2Q2)+272*A1*P2Q1/(3*P2Q2)+
0293      &112*A2*P2Q1/P2Q2-80*A1*A2*MB**2*P2Q1/P2Q2-
0294      &400*A1*A2*MB*MT*P2Q1/(3*P2Q2)+208*A2**2*MB*MT*P2Q1/(3*P2Q2)-
0295      &272*A1*A2*MT**2*P2Q1/(3*P2Q2)-320*A1*A2*P1P2*P2Q1/P2Q2+
0296      &96*A2**2*P1P2*P2Q1/P2Q2+256*A1*MB*MT**3*P2Q1/(3*P1Q1**2*P2Q2)+
0297      &512*A1*MT**2*P1P2*P2Q1/(3*P1Q1**2*P2Q2)-8*P2Q1/(3*P1Q1*P2Q2)-
0298      &200*A1*MB*MT*P2Q1/(3*P1Q1*P2Q2)-56*A2*MB*MT*P2Q1/(3*P1Q1*P2Q2)+
0299      &272*A1*MT**2*P2Q1/(3*P1Q1*P2Q2)+160*A2*MT**2*P2Q1/(3*P1Q1*P2Q2)+
0300      &256*A1*P1P2*P2Q1/(3*P1Q1*P2Q2)+48*A2*P1P2*P2Q1/(P1Q1*P2Q2)-
0301      &256*A1*A2*MB*MT*P1P2*P2Q1/(P1Q1*P2Q2)-
0302      &256*A1*A2*MT**2*P1P2*P2Q1/(3*P1Q1*P2Q2)-
0303      &1024*A1*A2*P1P2**2*P2Q1/(3*P1Q1*P2Q2)-
0304      &544*A1*A2*P1Q1*P2Q1/(3*P2Q2)-640*A2**2*P1Q1*P2Q1/(3*P2Q2)-
0305      &32*A1*P1P2*P2Q1/(3*P1Q2*P2Q2)
0306       V18=V18+32*A1*A2*MB*MT*P1P2*P2Q1/(3*P1Q2*P2Q2)+
0307      &32*A1*A2*MT**2*P1P2*P2Q1/(3*P1Q2*P2Q2)+
0308      &64*A1*A2*P1P2**2*P2Q1/(3*P1Q2*P2Q2)-
0309      &32*A1*MT**2*P1P2*P2Q1/(3*P1Q1*P1Q2*P2Q2)+
0310      &64*A1*A2*P1P2*P1Q1*P2Q1/(3*P1Q2*P2Q2)-
0311      &944*A1*A2*P1Q2*P2Q1/(3*P2Q2)-32*A2**2*P1Q2*P2Q1/P2Q2+
0312      &256*A1*MT**2*P1Q2*P2Q1/(3*P1Q1**2*P2Q2)+
0313      &96*A1*P1Q2*P2Q1/(P1Q1*P2Q2)+96*A2*P1Q2*P2Q1/(P1Q1*P2Q2)-
0314      &128*A1*A2*MB**2*P1Q2*P2Q1/(3*P1Q1*P2Q2)-
0315      &256*A1*A2*MB*MT*P1Q2*P2Q1/(P1Q1*P2Q2)-
0316      &128*A1*A2*MT**2*P1Q2*P2Q1/(3*P1Q1*P2Q2)-
0317      &512*A1*A2*P1P2*P1Q2*P2Q1/(P1Q1*P2Q2)-
0318      &512*A1*A2*P1Q2**2*P2Q1/(3*P1Q1*P2Q2)+544*A1*A2*P2Q1**2/(3*P2Q2)-
0319      &256*A1*MT**2*P2Q1**2/(3*P1Q1**2*P2Q2)-
0320      &272*A1*P2Q1**2/(3*P1Q1*P2Q2)+
0321      &256*A1*A2*MB*MT*P2Q1**2/(3*P1Q1*P2Q2)+
0322      &256*A1*A2*MT**2*P2Q1**2/(3*P1Q1*P2Q2)
0323       V18=V18+512*A1*A2*P1P2*P2Q1**2/(3*P1Q1*P2Q2)+
0324      &512*A1*A2*P1Q2*P2Q1**2/(3*P1Q1*P2Q2)+224*A12*P2Q2+
0325      &656*A1*A2*P2Q2/3+256*A12*MT**4*P2Q2/(3*P1Q1**2)+
0326      &16*A1*P2Q2/(3*P1Q1)+112*A1*A2*MB*MT*P2Q2/(3*P1Q1)+
0327      &32*A12*MT**2*P2Q2/(3*P1Q1)+64*A1*A2*MT**2*P2Q2/P1Q1+
0328      &368*A1*A2*P1P2*P2Q2/(3*P1Q1)-256*A1*MT**2*P2Q2/(3*P1Q2**2)+
0329      &256*A12*MT**4*P2Q2/(3*P1Q2**2)-256*A1*P2Q2/(3*P1Q2)+
0330      &224*A1*A2*MB*MT*P2Q2/(3*P1Q2)-368*A12*MT**2*P2Q2/(3*P1Q2)+
0331      &448*A1*A2*MT**2*P2Q2/(3*P1Q2)+304*A1*A2*P1P2*P2Q2/(3*P1Q2)+
0332      &16*A1*MT**2*P2Q2/(3*P1Q1*P1Q2)-64*A12*MT**4*P2Q2/(3*P1Q1*P1Q2)+
0333      &32*A12*P1Q1*P2Q2/P1Q2+944*A1*A2*P1Q1*P2Q2/(3*P1Q2)+
0334      &256*A12*MT**2*P1Q2*P2Q2/(3*P1Q1**2)+
0335      &640*A12*P1Q2*P2Q2/(3*P1Q1)+544*A1*A2*P1Q2*P2Q2/(3*P1Q1)-
0336      &256*A2*MB**2*P2Q2/(3*P2Q1**2)-256*A2**2*MB**3*MT*P2Q2/(3*P2Q1**2)+
0337      &64*MB**2*MT**2*P2Q2/(3*P1Q2**2*P2Q1**2)+
0338      &64*MB**2*P2Q2/(3*P1Q2*P2Q1**2)-
0339      &128*A2*MB**3*MT*P2Q2/(3*P1Q2*P2Q1**2)
0340       V18=V18-128*A2*MB**2*MT**2*P2Q2/(3*P1Q2*P2Q1**2)-
0341      &256*A2*MB**2*P1P2*P2Q2/(3*P1Q2*P2Q1**2)-
0342      &256*A2*MB**2*P1Q1*P2Q2/(3*P1Q2*P2Q1**2)+
0343      &256*A2**2*MB**2*P1Q2*P2Q2/(3*P2Q1**2)+272*A1*P2Q2/(3*P2Q1)+
0344      &112*A2*P2Q2/P2Q1-80*A1*A2*MB**2*P2Q2/P2Q1-
0345      &400*A1*A2*MB*MT*P2Q2/(3*P2Q1)+208*A2**2*MB*MT*P2Q2/(3*P2Q1)-
0346      &272*A1*A2*MT**2*P2Q2/(3*P2Q1)-320*A1*A2*P1P2*P2Q2/P2Q1+
0347      &96*A2**2*P1P2*P2Q2/P2Q1-32*A1*P1P2*P2Q2/(3*P1Q1*P2Q1)+
0348      &32*A1*A2*MB*MT*P1P2*P2Q2/(3*P1Q1*P2Q1)+
0349      &32*A1*A2*MT**2*P1P2*P2Q2/(3*P1Q1*P2Q1)+
0350      &64*A1*A2*P1P2**2*P2Q2/(3*P1Q1*P2Q1)-944*A1*A2*P1Q1*P2Q2/(3*P2Q1)-
0351      &32*A2**2*P1Q1*P2Q2/P2Q1+256*A1*MB*MT**3*P2Q2/(3*P1Q2**2*P2Q1)+
0352      &512*A1*MT**2*P1P2*P2Q2/(3*P1Q2**2*P2Q1)+
0353      &256*A1*MT**2*P1Q1*P2Q2/(3*P1Q2**2*P2Q1)-8*P2Q2/(3*P1Q2*P2Q1)-
0354      &200*A1*MB*MT*P2Q2/(3*P1Q2*P2Q1)-56*A2*MB*MT*P2Q2/(3*P1Q2*P2Q1)+
0355      &272*A1*MT**2*P2Q2/(3*P1Q2*P2Q1)+160*A2*MT**2*P2Q2/(3*P1Q2*P2Q1)+
0356      &256*A1*P1P2*P2Q2/(3*P1Q2*P2Q1)+48*A2*P1P2*P2Q2/(P1Q2*P2Q1)
0357       V18=V18-256*A1*A2*MB*MT*P1P2*P2Q2/(P1Q2*P2Q1)-
0358      &256*A1*A2*MT**2*P1P2*P2Q2/(3*P1Q2*P2Q1)-
0359      &1024*A1*A2*P1P2**2*P2Q2/(3*P1Q2*P2Q1)-
0360      &32*A1*MT**2*P1P2*P2Q2/(3*P1Q1*P1Q2*P2Q1)+
0361      &96*A1*P1Q1*P2Q2/(P1Q2*P2Q1)+96*A2*P1Q1*P2Q2/(P1Q2*P2Q1)-
0362      &128*A1*A2*MB**2*P1Q1*P2Q2/(3*P1Q2*P2Q1)-
0363      &256*A1*A2*MB*MT*P1Q1*P2Q2/(P1Q2*P2Q1)-
0364      &128*A1*A2*MT**2*P1Q1*P2Q2/(3*P1Q2*P2Q1)-
0365      &512*A1*A2*P1P2*P1Q1*P2Q2/(P1Q2*P2Q1)-
0366      &512*A1*A2*P1Q1**2*P2Q2/(3*P1Q2*P2Q1)-544*A1*A2*P1Q2*P2Q2/(3*P2Q1)-
0367      &640*A2**2*P1Q2*P2Q2/(3*P2Q1)+
0368      &64*A1*A2*P1P2*P1Q2*P2Q2/(3*P1Q1*P2Q1)+544*A1*A2*P2Q2**2/(3*P2Q1)-
0369      &256*A1*MT**2*P2Q2**2/(3*P1Q2**2*P2Q1)-
0370      &272*A1*P2Q2**2/(3*P1Q2*P2Q1)+
0371      &256*A1*A2*MB*MT*P2Q2**2/(3*P1Q2*P2Q1)+
0372      &256*A1*A2*MT**2*P2Q2**2/(3*P1Q2*P2Q1)+
0373      &512*A1*A2*P1P2*P2Q2**2/(3*P1Q2*P2Q1)
0374       V18=V18+512*A1*A2*P1Q1*P2Q2**2/(3*P1Q2*P2Q1)+
0375      &384*A12*MB*MT*P1Q1**2/S**2+
0376      &384*A12*P1P2*P1Q1**2/S**2+2688*A12*MB*MT*P1Q1*P1Q2/S**2+
0377      &2688*A12*P1P2*P1Q1*P1Q2/S**2+384*A12*MB*MT*P1Q2**2/S**2+
0378      &384*A12*P1P2*P1Q2**2/S**2+768*A1*A2*MB*MT*P1Q1*P2Q1/S**2+
0379      &768*A1*A2*P1P2*P1Q1*P2Q1/S**2+2688*A1*A2*MB*MT*P1Q2*P2Q1/S**2+
0380      &2688*A1*A2*P1P2*P1Q2*P2Q1/S**2-960*A12*P1Q1*P1Q2*P2Q1/S**2-
0381      &960*A1*A2*P1Q1*P1Q2*P2Q1/S**2+960*A12*P1Q2**2*P2Q1/S**2+
0382      &960*A1*A2*P1Q2**2*P2Q1/S**2+384*A2**2*MB*MT*P2Q1**2/S**2+
0383      &384*A2**2*P1P2*P2Q1**2/S**2-960*A1*A2*P1Q2*P2Q1**2/S**2-
0384      &960*A2**2*P1Q2*P2Q1**2/S**2+2688*A1*A2*MB*MT*P1Q1*P2Q2/S**2+
0385      &2688*A1*A2*P1P2*P1Q1*P2Q2/S**2+960*A12*P1Q1**2*P2Q2/S**2+
0386      &960*A1*A2*P1Q1**2*P2Q2/S**2+768*A1*A2*MB*MT*P1Q2*P2Q2/S**2+
0387      &768*A1*A2*P1P2*P1Q2*P2Q2/S**2-960*A12*P1Q1*P1Q2*P2Q2/S**2-
0388      &960*A1*A2*P1Q1*P1Q2*P2Q2/S**2+2688*A2**2*MB*MT*P2Q1*P2Q2/S**2+
0389      &2688*A2**2*P1P2*P2Q1*P2Q2/S**2+960*A1*A2*P1Q1*P2Q1*P2Q2/S**2+
0390      &960*A2**2*P1Q1*P2Q1*P2Q2/S**2+960*A1*A2*P1Q2*P2Q1*P2Q2/S**2+
0391      &960*A2**2*P1Q2*P2Q1*P2Q2/S**2+384*A2**2*MB*MT*P2Q2**2/S**2
0392       V18=V18+384*A2**2*P1P2*P2Q2**2/S**2-960*A1*A2*P1Q1*P2Q2**2/S**2-
0393      &960*A2**2*P1Q1*P2Q2**2/S**2+96*A1*MB*MT/S+96*A2*MB*MT/S-
0394      &768*A2**2*MB**3*MT/S-768*A12*MB*MT**3/S-192*A1*P1P2/S-
0395      &192*A2*P1P2/S-768*A2**2*MB**2*P1P2/S-2304*A1*A2*MB*MT*P1P2/S-
0396      &768*A12*MT**2*P1P2/S-2304*A1*A2*P1P2**2/S-
0397      &96*A1*MB*MT**3/(P1Q1*S)-192*A2*MB*MT*P1P2/(P1Q1*S)-
0398      &96*A1*MT**2*P1P2/(P1Q1*S)-192*A2*P1P2**2/(P1Q1*S)-192*A1*P1Q1/S-
0399      &144*A2*P1Q1/S-384*A1*A2*MB**2*P1Q1/S-480*A2**2*MB**2*P1Q1/S-
0400      &480*A12*MB*MT*P1Q1/S+96*A1*A2*MB*MT*P1Q1/S-
0401      &864*A12*P1P2*P1Q1/S-672*A1*A2*P1P2*P1Q1/S-96*A1*A2*P1Q1**2/S-
0402      &96*A1*MB*MT**3/(P1Q2*S)-192*A2*MB*MT*P1P2/(P1Q2*S)-
0403      &96*A1*MT**2*P1P2/(P1Q2*S)-192*A2*P1P2**2/(P1Q2*S)-
0404      &48*A1*MB*MT*P1Q1/(P1Q2*S)+96*A2*MB*MT*P1Q1/(P1Q2*S)-
0405      &48*A1*MT**2*P1Q1/(P1Q2*S)-192*A1*P1P2*P1Q1/(P1Q2*S)-
0406      &192*A2*P1P2*P1Q1/(P1Q2*S)+192*A1*A2*MB*MT*P1P2*P1Q1/(P1Q2*S)+
0407      &192*A1*A2*P1P2**2*P1Q1/(P1Q2*S)-192*A1*P1Q1**2/(P1Q2*S)-
0408      &192*A2*P1Q1**2/(P1Q2*S)+192*A1*A2*MB**2*P1Q1**2/(P1Q2*S)
0409       V18=V18-192*A12*MB*MT*P1Q1**2/(P1Q2*S)+
0410      &96*A1*A2*MB*MT*P1Q1**2/(P1Q2*S)+
0411      &192*A1*A2*P1P2*P1Q1**2/(P1Q2*S)-192*A1*P1Q2/S-144*A2*P1Q2/S-
0412      &384*A1*A2*MB**2*P1Q2/S-480*A2**2*MB**2*P1Q2/S-
0413      &480*A12*MB*MT*P1Q2/S+96*A1*A2*MB*MT*P1Q2/S-
0414      &864*A12*P1P2*P1Q2/S-672*A1*A2*P1P2*P1Q2/S-
0415      &48*A1*MB*MT*P1Q2/(P1Q1*S)+96*A2*MB*MT*P1Q2/(P1Q1*S)-
0416      &48*A1*MT**2*P1Q2/(P1Q1*S)-192*A1*P1P2*P1Q2/(P1Q1*S)-
0417      &192*A2*P1P2*P1Q2/(P1Q1*S)+192*A1*A2*MB*MT*P1P2*P1Q2/(P1Q1*S)+
0418      &192*A1*A2*P1P2**2*P1Q2/(P1Q1*S)-576*A1*A2*P1Q1*P1Q2/S-
0419      &96*A1*A2*P1Q2**2/S-192*A1*P1Q2**2/(P1Q1*S)-
0420      &192*A2*P1Q2**2/(P1Q1*S)+192*A1*A2*MB**2*P1Q2**2/(P1Q1*S)-
0421      &192*A12*MB*MT*P1Q2**2/(P1Q1*S)+96*A1*A2*MB*MT*P1Q2**2/(P1Q1*S)+
0422      &192*A1*A2*P1P2*P1Q2**2/(P1Q1*S)+96*A2*MB**3*MT/(P2Q1*S)+
0423      &96*A2*MB**2*P1P2/(P2Q1*S)+192*A1*MB*MT*P1P2/(P2Q1*S)+
0424      &192*A1*P1P2**2/(P2Q1*S)+96*A1*MB**2*P1Q1/(P2Q1*S)+
0425      &192*A2*MB**2*P1Q1/(P2Q1*S)+96*A1*MB*MT*P1Q1/(P2Q1*S)+
0426      &192*A1*A2*MB**3*MT*P1Q1/(P2Q1*S)+192*A1*P1P2*P1Q1/(P2Q1*S)
0427       V18=V18+192*A1*A2*MB**2*P1P2*P1Q1/(P2Q1*S)+
0428      &96*A1*A2*MB**2*P1Q1**2/(P2Q1*S)+
0429      &192*A2*MB**3*MT*P1Q1/(P1Q2*P2Q1*S)+
0430      &192*A2*MB**2*P1P2*P1Q1/(P1Q2*P2Q1*S)+
0431      &96*A1*MB*MT*P1P2*P1Q1/(P1Q2*P2Q1*S)+
0432      &96*A1*P1P2**2*P1Q1/(P1Q2*P2Q1*S)+
0433      &96*A1*MB**2*P1Q1**2/(P1Q2*P2Q1*S)+
0434      &192*A2*MB**2*P1Q1**2/(P1Q2*P2Q1*S)+
0435      &48*A1*MB*MT*P1Q1**2/(P1Q2*P2Q1*S)+
0436      &96*A1*P1P2*P1Q1**2/(P1Q2*P2Q1*S)+96*A1*MB**2*P1Q2/(P2Q1*S)+
0437      &48*A2*MB**2*P1Q2/(P2Q1*S)-192*A1*A2*MB**3*MT*P1Q2/(P2Q1*S)-
0438      &192*A1*A2*MB**2*P1P2*P1Q2/(P2Q1*S)-
0439      &96*A1*A2*MB**2*P1Q2**2/(P2Q1*S)+144*A1*P2Q1/S+192*A2*P2Q1/S-
0440      &96*A1*A2*MB*MT*P2Q1/S+480*A2**2*MB*MT*P2Q1/S+
0441      &480*A12*MT**2*P2Q1/S+384*A1*A2*MT**2*P2Q1/S+
0442      &672*A1*A2*P1P2*P2Q1/S+864*A2**2*P1P2*P2Q1/S+
0443      &96*A2*MB*MT*P2Q1/(P1Q1*S)+192*A1*MT**2*P2Q1/(P1Q1*S)
0444       V18=V18+96*A2*MT**2*P2Q1/(P1Q1*S)+
0445      &192*A1*A2*MB*MT**3*P2Q1/(P1Q1*S)+
0446      &192*A2*P1P2*P2Q1/(P1Q1*S)+192*A1*A2*MT**2*P1P2*P2Q1/(P1Q1*S)-
0447      &192*A12*P1Q1*P2Q1/S-192*A2**2*P1Q1*P2Q1/S+
0448      &48*A1*MT**2*P2Q1/(P1Q2*S)+96*A2*MT**2*P2Q1/(P1Q2*S)-
0449      &192*A1*A2*MB*MT**3*P2Q1/(P1Q2*S)-
0450      &192*A1*A2*MT**2*P1P2*P2Q1/(P1Q2*S)-
0451      &96*A1*A2*MB*MT*P1Q1*P2Q1/(P1Q2*S)-
0452      &192*A12*MT**2*P1Q1*P2Q1/(P1Q2*S)-
0453      &96*A1*A2*MT**2*P1Q1*P2Q1/(P1Q2*S)-
0454      &384*A1*A2*P1P2*P1Q1*P2Q1/(P1Q2*S)-384*A12*P1Q1**2*P2Q1/(P1Q2*S)-
0455      &384*A1*A2*P1Q1**2*P2Q1/(P1Q2*S)-480*A12*P1Q2*P2Q1/S-
0456      &960*A1*A2*P1Q2*P2Q1/S-480*A2**2*P1Q2*P2Q1/S+
0457      &144*A1*P1Q2*P2Q1/(P1Q1*S)+96*A2*P1Q2*P2Q1/(P1Q1*S)-
0458      &384*A1*A2*MB*MT*P1Q2*P2Q1/(P1Q1*S)-
0459      &96*A12*MT**2*P1Q2*P2Q1/(P1Q1*S)+
0460      &96*A1*A2*MT**2*P1Q2*P2Q1/(P1Q1*S)-
0461      &576*A1*A2*P1P2*P1Q2*P2Q1/(P1Q1*S)-192*A12*P1Q2**2*P2Q1/(P1Q1*S)
0462       V18=V18-384*A1*A2*P1Q2**2*P2Q1/(P1Q1*S)-96*A1*A2*P2Q1**2/S-
0463      &96*A1*A2*MT**2*P2Q1**2/(P1Q1*S)+96*A1*A2*MT**2*P2Q1**2/(P1Q2*S)+
0464      &288*A1*A2*P1Q2*P2Q1**2/(P1Q1*S)+96*A2*MB**3*MT/(P2Q2*S)+
0465      &96*A2*MB**2*P1P2/(P2Q2*S)+192*A1*MB*MT*P1P2/(P2Q2*S)+
0466      &192*A1*P1P2**2/(P2Q2*S)+96*A1*MB**2*P1Q1/(P2Q2*S)+
0467      &48*A2*MB**2*P1Q1/(P2Q2*S)-192*A1*A2*MB**3*MT*P1Q1/(P2Q2*S)-
0468      &192*A1*A2*MB**2*P1P2*P1Q1/(P2Q2*S)-
0469      &96*A1*A2*MB**2*P1Q1**2/(P2Q2*S)+96*A1*MB**2*P1Q2/(P2Q2*S)+
0470      &192*A2*MB**2*P1Q2/(P2Q2*S)+96*A1*MB*MT*P1Q2/(P2Q2*S)+
0471      &192*A1*A2*MB**3*MT*P1Q2/(P2Q2*S)+192*A1*P1P2*P1Q2/(P2Q2*S)+
0472      &192*A1*A2*MB**2*P1P2*P1Q2/(P2Q2*S)+
0473      &192*A2*MB**3*MT*P1Q2/(P1Q1*P2Q2*S)+
0474      &192*A2*MB**2*P1P2*P1Q2/(P1Q1*P2Q2*S)+
0475      &96*A1*MB*MT*P1P2*P1Q2/(P1Q1*P2Q2*S)+
0476      &96*A1*P1P2**2*P1Q2/(P1Q1*P2Q2*S)+96*A1*A2*MB**2*P1Q2**2/(P2Q2*S)+
0477      &96*A1*MB**2*P1Q2**2/(P1Q1*P2Q2*S)+
0478      &192*A2*MB**2*P1Q2**2/(P1Q1*P2Q2*S)
0479       V18=V18+48*A1*MB*MT*P1Q2**2/(P1Q1*P2Q2*S)+
0480      &96*A1*P1P2*P1Q2**2/(P1Q1*P2Q2*S)-48*A2*MB**2*P2Q1/(P2Q2*S)+
0481      &96*A1*MB*MT*P2Q1/(P2Q2*S)-48*A2*MB*MT*P2Q1/(P2Q2*S)-
0482      &192*A1*P1P2*P2Q1/(P2Q2*S)-192*A2*P1P2*P2Q1/(P2Q2*S)+
0483      &192*A1*A2*MB*MT*P1P2*P2Q1/(P2Q2*S)+
0484      &192*A1*A2*P1P2**2*P2Q1/(P2Q2*S)-
0485      &192*A1*MB*MT**3*P2Q1/(P1Q1*P2Q2*S)-
0486      &96*A2*MB*MT*P1P2*P2Q1/(P1Q1*P2Q2*S)-
0487      &192*A1*MT**2*P1P2*P2Q1/(P1Q1*P2Q2*S)-
0488      &96*A2*P1P2**2*P2Q1/(P1Q1*P2Q2*S)+
0489      &96*A1*A2*MB**2*P1Q1*P2Q1/(P2Q2*S)+
0490      &192*A2**2*MB**2*P1Q1*P2Q1/(P2Q2*S)+
0491      &96*A1*A2*MB*MT*P1Q1*P2Q1/(P2Q2*S)+
0492      &384*A1*A2*P1P2*P1Q1*P2Q1/(P2Q2*S)-96*A1*P1Q2*P2Q1/(P2Q2*S)-
0493      &144*A2*P1Q2*P2Q1/(P2Q2*S)-96*A1*A2*MB**2*P1Q2*P2Q1/(P2Q2*S)+
0494      &96*A2**2*MB**2*P1Q2*P2Q1/(P2Q2*S)+
0495      &384*A1*A2*MB*MT*P1Q2*P2Q1/(P2Q2*S)
0496       V18=V18+576*A1*A2*P1P2*P1Q2*P2Q1/(P2Q2*S)-
0497      &96*A2*MB**2*P1Q2*P2Q1/(P1Q1*P2Q2*S)+
0498      &48*A1*MB*MT*P1Q2*P2Q1/(P1Q1*P2Q2*S)+
0499      &48*A2*MB*MT*P1Q2*P2Q1/(P1Q1*P2Q2*S)-
0500      &96*A1*MT**2*P1Q2*P2Q1/(P1Q1*P2Q2*S)-
0501      &96*A1*P1P2*P1Q2*P2Q1/(P1Q1*P2Q2*S)-
0502      &96*A2*P1P2*P1Q2*P2Q1/(P1Q1*P2Q2*S)+
0503      &96*A1*A2*P1Q1*P1Q2*P2Q1/(P2Q2*S)+288*A1*A2*P1Q2**2*P2Q1/(P2Q2*S)-
0504      &96*A1*P1Q2**2*P2Q1/(P1Q1*P2Q2*S)-96*A2*P1Q2**2*P2Q1/(P1Q1*P2Q2*S)+
0505      &192*A1*P2Q1**2/(P2Q2*S)+192*A2*P2Q1**2/(P2Q2*S)-
0506      &96*A1*A2*MB*MT*P2Q1**2/(P2Q2*S)+192*A2**2*MB*MT*P2Q1**2/(P2Q2*S)-
0507      &192*A1*A2*MT**2*P2Q1**2/(P2Q2*S)-192*A1*A2*P1P2*P2Q1**2/(P2Q2*S)+
0508      &48*A2*MB*MT*P2Q1**2/(P1Q1*P2Q2*S)+
0509      &192*A1*MT**2*P2Q1**2/(P1Q1*P2Q2*S)+
0510      &96*A2*MT**2*P2Q1**2/(P1Q1*P2Q2*S)+
0511      &96*A2*P1P2*P2Q1**2/(P1Q1*P2Q2*S)-384*A1*A2*P1Q1*P2Q1**2/(P2Q2*S)-
0512      &384*A2**2*P1Q1*P2Q1**2/(P2Q2*S)-384*A1*A2*P1Q2*P2Q1**2/(P2Q2*S)
0513       V18=V18-192*A2**2*P1Q2*P2Q1**2/(P2Q2*S)+
0514      &96*A1*P1Q2*P2Q1**2/(P1Q1*P2Q2*S)+
0515      &96*A2*P1Q2*P2Q1**2/(P1Q1*P2Q2*S)+144*A1*P2Q2/S+192*A2*P2Q2/S-
0516      &96*A1*A2*MB*MT*P2Q2/S+480*A2**2*MB*MT*P2Q2/S+
0517      &480*A12*MT**2*P2Q2/S+384*A1*A2*MT**2*P2Q2/S+
0518      &672*A1*A2*P1P2*P2Q2/S+864*A2**2*P1P2*P2Q2/S+
0519      &48*A1*MT**2*P2Q2/(P1Q1*S)+96*A2*MT**2*P2Q2/(P1Q1*S)-
0520      &192*A1*A2*MB*MT**3*P2Q2/(P1Q1*S)-
0521      &192*A1*A2*MT**2*P1P2*P2Q2/(P1Q1*S)-480*A12*P1Q1*P2Q2/S-
0522      &960*A1*A2*P1Q1*P2Q2/S-480*A2**2*P1Q1*P2Q2/S+
0523      &96*A2*MB*MT*P2Q2/(P1Q2*S)+192*A1*MT**2*P2Q2/(P1Q2*S)+
0524      &96*A2*MT**2*P2Q2/(P1Q2*S)+192*A1*A2*MB*MT**3*P2Q2/(P1Q2*S)+
0525      &192*A2*P1P2*P2Q2/(P1Q2*S)+192*A1*A2*MT**2*P1P2*P2Q2/(P1Q2*S)+
0526      &144*A1*P1Q1*P2Q2/(P1Q2*S)+96*A2*P1Q1*P2Q2/(P1Q2*S)-
0527      &384*A1*A2*MB*MT*P1Q1*P2Q2/(P1Q2*S)-
0528      &96*A12*MT**2*P1Q1*P2Q2/(P1Q2*S)+
0529      &96*A1*A2*MT**2*P1Q1*P2Q2/(P1Q2*S)
0530       V18=V18-576*A1*A2*P1P2*P1Q1*P2Q2/(P1Q2*S)-
0531      &192*A12*P1Q1**2*P2Q2/(P1Q2*S)-
0532      &384*A1*A2*P1Q1**2*P2Q2/(P1Q2*S)-192*A12*P1Q2*P2Q2/S-
0533      &192*A2**2*P1Q2*P2Q2/S-96*A1*A2*MB*MT*P1Q2*P2Q2/(P1Q1*S)-
0534      &192*A12*MT**2*P1Q2*P2Q2/(P1Q1*S)-
0535      &96*A1*A2*MT**2*P1Q2*P2Q2/(P1Q1*S)-
0536      &384*A1*A2*P1P2*P1Q2*P2Q2/(P1Q1*S)-384*A12*P1Q2**2*P2Q2/(P1Q1*S)-
0537      &384*A1*A2*P1Q2**2*P2Q2/(P1Q1*S)-48*A2*MB**2*P2Q2/(P2Q1*S)+
0538      &96*A1*MB*MT*P2Q2/(P2Q1*S)-48*A2*MB*MT*P2Q2/(P2Q1*S)-
0539      &192*A1*P1P2*P2Q2/(P2Q1*S)-192*A2*P1P2*P2Q2/(P2Q1*S)+
0540      &192*A1*A2*MB*MT*P1P2*P2Q2/(P2Q1*S)+
0541      &192*A1*A2*P1P2**2*P2Q2/(P2Q1*S)-96*A1*P1Q1*P2Q2/(P2Q1*S)-
0542      &144*A2*P1Q1*P2Q2/(P2Q1*S)-96*A1*A2*MB**2*P1Q1*P2Q2/(P2Q1*S)+
0543      &96*A2**2*MB**2*P1Q1*P2Q2/(P2Q1*S)+
0544      &384*A1*A2*MB*MT*P1Q1*P2Q2/(P2Q1*S)+
0545      &576*A1*A2*P1P2*P1Q1*P2Q2/(P2Q1*S)+288*A1*A2*P1Q1**2*P2Q2/(P2Q1*S)-
0546      &192*A1*MB*MT**3*P2Q2/(P1Q2*P2Q1*S)
0547       V18=V18-96*A2*MB*MT*P1P2*P2Q2/(P1Q2*P2Q1*S)-
0548      &192*A1*MT**2*P1P2*P2Q2/(P1Q2*P2Q1*S)-
0549      &96*A2*P1P2**2*P2Q2/(P1Q2*P2Q1*S)-
0550      &96*A2*MB**2*P1Q1*P2Q2/(P1Q2*P2Q1*S)+
0551      &48*A1*MB*MT*P1Q1*P2Q2/(P1Q2*P2Q1*S)
0552  
0553       V18BIS=
0554      &48*A2*MB*MT*P1Q1*P2Q2/(P1Q2*P2Q1*S)-
0555      &96*A1*MT**2*P1Q1*P2Q2/(P1Q2*P2Q1*S)-
0556      &96*A1*P1P2*P1Q1*P2Q2/(P1Q2*P2Q1*S)-
0557      &96*A2*P1P2*P1Q1*P2Q2/(P1Q2*P2Q1*S)-
0558      &96*A1*P1Q1**2*P2Q2/(P1Q2*P2Q1*S)-96*A2*P1Q1**2*P2Q2/(P1Q2*P2Q1*S)+
0559      &96*A1*A2*MB**2*P1Q2*P2Q2/(P2Q1*S)+
0560      &192*A2**2*MB**2*P1Q2*P2Q2/(P2Q1*S)+
0561      &96*A1*A2*MB*MT*P1Q2*P2Q2/(P2Q1*S)+
0562      &384*A1*A2*P1P2*P1Q2*P2Q2/(P2Q1*S)+
0563      &96*A1*A2*P1Q1*P1Q2*P2Q2/(P2Q1*S)-576*A1*A2*P2Q1*P2Q2/S+
0564      &96*A1*A2*P1Q1*P2Q1*P2Q2/(P1Q2*S)+96*A1*A2*P1Q2*P2Q1*P2Q2/(P1Q1*S)-
0565      &96*A1*A2*P2Q2**2/S+96*A1*A2*MT**2*P2Q2**2/(P1Q1*S)-
0566      &96*A1*A2*MT**2*P2Q2**2/(P1Q2*S)+288*A1*A2*P1Q1*P2Q2**2/(P1Q2*S)+
0567      &192*A1*P2Q2**2/(P2Q1*S)+192*A2*P2Q2**2/(P2Q1*S)-
0568      &96*A1*A2*MB*MT*P2Q2**2/(P2Q1*S)+192*A2**2*MB*MT*P2Q2**2/(P2Q1*S)-
0569      &192*A1*A2*MT**2*P2Q2**2/(P2Q1*S)-192*A1*A2*P1P2*P2Q2**2/(P2Q1*S)
0570       V18BIS=V18BIS-384*A1*A2*P1Q1*P2Q2**2/(P2Q1*S)-
0571      &192*A2**2*P1Q1*P2Q2**2/(P2Q1*S)+
0572      &48*A2*MB*MT*P2Q2**2/(P1Q2*P2Q1*S)+
0573      &192*A1*MT**2*P2Q2**2/(P1Q2*P2Q1*S)+
0574      &96*A2*MT**2*P2Q2**2/(P1Q2*P2Q1*S)+
0575      &96*A2*P1P2*P2Q2**2/(P1Q2*P2Q1*S)+96*A1*P1Q1*P2Q2**2/(P1Q2*P2Q1*S)+
0576      &96*A2*P1Q1*P2Q2**2/(P1Q2*P2Q1*S)-384*A1*A2*P1Q2*P2Q2**2/(P2Q1*S)-
0577      &384*A2**2*P1Q2*P2Q2**2/(P2Q1*S)+512*A1*A2*S/3-
0578      &128*A1*MT**2*S/(3*P1Q1**2)-128*A12*MB*MT**3*S/(3*P1Q1**2)-
0579      &152*A1*S/(3*P1Q1)+152*A12*MB*MT*S/(3*P1Q1)+
0580      &128*A1*A2*MB*MT*S/(3*P1Q1)+112*A1*A2*MT**2*S/(3*P1Q1)-
0581      &16*A12*P1P2*S/P1Q1+152*A1*A2*P1P2*S/(3*P1Q1)-
0582      &128*A1*MT**2*S/(3*P1Q2**2)-128*A12*MB*MT**3*S/(3*P1Q2**2)-
0583      &152*A1*S/(3*P1Q2)+152*A12*MB*MT*S/(3*P1Q2)+
0584      &128*A1*A2*MB*MT*S/(3*P1Q2)+112*A1*A2*MT**2*S/(3*P1Q2)-
0585      &16*A12*P1P2*S/P1Q2+152*A1*A2*P1P2*S/(3*P1Q2)-
0586      &16*A1*MB*MT*S/(3*P1Q1*P1Q2)+32*A12*MB*MT**3*S/(3*P1Q1*P1Q2)
0587       V18BIS=V18BIS-16*A1*P1P2*S/(3*P1Q1*P1Q2)+
0588      &272*A1*A2*P1Q1*S/(3*P1Q2)+
0589      &272*A1*A2*P1Q2*S/(3*P1Q1)-128*A2*MB**2*S/(3*P2Q1**2)-
0590      &128*A2**2*MB**3*MT*S/(3*P2Q1**2)+
0591      &32*MB**2*MT**2*S/(3*P1Q2**2*P2Q1**2)+32*MB**2*S/(3*P1Q2*P2Q1**2)-
0592      &64*A2*MB**3*MT*S/(3*P1Q2*P2Q1**2)-
0593      &64*A2*MB**2*MT**2*S/(3*P1Q2*P2Q1**2)-
0594      &128*A2*MB**2*P1P2*S/(3*P1Q2*P2Q1**2)-
0595      &128*A2*MB**2*P1Q1*S/(3*P1Q2*P2Q1**2)+
0596      &128*A2**2*MB**2*P1Q2*S/(3*P2Q1**2)+152*A2*S/(3*P2Q1)-
0597      &112*A1*A2*MB**2*S/(3*P2Q1)-128*A1*A2*MB*MT*S/(3*P2Q1)-
0598      &152*A2**2*MB*MT*S/(3*P2Q1)-152*A1*A2*P1P2*S/(3*P2Q1)+
0599      &16*A2**2*P1P2*S/P2Q1+8*A1*A2*MB**3*MT*S/(3*P1Q1*P2Q1)+
0600      &16*A1*A2*MB**2*MT**2*S/(3*P1Q1*P2Q1)+
0601      &8*A1*A2*MB*MT**3*S/(3*P1Q1*P2Q1)-8*A1*P1P2*S/(3*P1Q1*P2Q1)-
0602      &8*A2*P1P2*S/(3*P1Q1*P2Q1)+8*A1*A2*MB**2*P1P2*S/(3*P1Q1*P2Q1)+
0603      &16*A1*A2*MB*MT*P1P2*S/(3*P1Q1*P2Q1)
0604       V18BIS=V18BIS+8*A1*A2*MT**2*P1P2*S/(3*P1Q1*P2Q1)+
0605      &32*A1*A2*P1P2**2*S/(3*P1Q1*P2Q1)-32*A2**2*P1Q1*S/(3*P2Q1)-
0606      &32*MT**2*S/(3*P1Q2**2*P2Q1)+64*A1*MB**2*MT**2*S/(3*P1Q2**2*P2Q1)+
0607      &64*A1*MB*MT**3*S/(3*P1Q2**2*P2Q1)+
0608      &128*A1*MT**2*P1P2*S/(3*P1Q2**2*P2Q1)-12*S/(P1Q2*P2Q1)+
0609      &24*A1*MB**2*S/(P1Q2*P2Q1)-64*A1*A2*MB**3*MT*S/(3*P1Q2*P2Q1)+
0610      &24*A2*MT**2*S/(P1Q2*P2Q1)-128*A1*A2*MB**2*MT**2*S/(3*P1Q2*P2Q1)-
0611      &64*A1*A2*MB*MT**3*S/(3*P1Q2*P2Q1)+56*A1*P1P2*S/(3*P1Q2*P2Q1)+
0612      &56*A2*P1P2*S/(3*P1Q2*P2Q1)-64*A1*A2*MB**2*P1P2*S/(3*P1Q2*P2Q1)-
0613      &128*A1*A2*MB*MT*P1P2*S/(3*P1Q2*P2Q1)-
0614      &64*A1*A2*MT**2*P1P2*S/(3*P1Q2*P2Q1)-
0615      &256*A1*A2*P1P2**2*S/(3*P1Q2*P2Q1)+4*P1P2*S/(3*P1Q1*P1Q2*P2Q1)+
0616      &8*A1*MB*MT*P1P2*S/(3*P1Q1*P1Q2*P2Q1)-
0617      &8*A1*MT**2*P1P2*S/(3*P1Q1*P1Q2*P2Q1)+136*A2*P1Q1*S/(3*P1Q2*P2Q1)-
0618      &128*A1*A2*MB**2*P1Q1*S/(3*P1Q2*P2Q1)-
0619      &128*A1*A2*MB*MT*P1Q1*S/(3*P1Q2*P2Q1)-
0620      &256*A1*A2*P1P2*P1Q1*S/(3*P1Q2*P2Q1)-160*A2**2*P1Q2*S/(3*P2Q1)
0621       V18BIS=V18BIS+16*A1*A2*P1P2*P1Q2*S/(3*P1Q1*P2Q1)-
0622      &32*A12*P2Q1*S/(3*P1Q1)-
0623      &128*A12*MT**2*P2Q1*S/(3*P1Q2**2)-160*A12*P2Q1*S/(3*P1Q2)-
0624      &128*A2*MB**2*S/(3*P2Q2**2)-128*A2**2*MB**3*MT*S/(3*P2Q2**2)+
0625      &32*MB**2*MT**2*S/(3*P1Q1**2*P2Q2**2)+32*MB**2*S/(3*P1Q1*P2Q2**2)-
0626      &64*A2*MB**3*MT*S/(3*P1Q1*P2Q2**2)-
0627      &64*A2*MB**2*MT**2*S/(3*P1Q1*P2Q2**2)-
0628      &128*A2*MB**2*P1P2*S/(3*P1Q1*P2Q2**2)+
0629      &128*A2**2*MB**2*P1Q1*S/(3*P2Q2**2)-
0630      &128*A2*MB**2*P1Q2*S/(3*P1Q1*P2Q2**2)+152*A2*S/(3*P2Q2)-
0631      &112*A1*A2*MB**2*S/(3*P2Q2)-128*A1*A2*MB*MT*S/(3*P2Q2)-
0632      &152*A2**2*MB*MT*S/(3*P2Q2)-152*A1*A2*P1P2*S/(3*P2Q2)+
0633      &16*A2**2*P1P2*S/P2Q2-32*MT**2*S/(3*P1Q1**2*P2Q2)+
0634      &64*A1*MB**2*MT**2*S/(3*P1Q1**2*P2Q2)+
0635      &64*A1*MB*MT**3*S/(3*P1Q1**2*P2Q2)+
0636      &128*A1*MT**2*P1P2*S/(3*P1Q1**2*P2Q2)-12*S/(P1Q1*P2Q2)+
0637      &24*A1*MB**2*S/(P1Q1*P2Q2)-64*A1*A2*MB**3*MT*S/(3*P1Q1*P2Q2)
0638       V18BIS=V18BIS+24*A2*MT**2*S/(P1Q1*P2Q2)-
0639      &128*A1*A2*MB**2*MT**2*S/(3*P1Q1*P2Q2)-
0640      &64*A1*A2*MB*MT**3*S/(3*P1Q1*P2Q2)+56*A1*P1P2*S/(3*P1Q1*P2Q2)+
0641      &56*A2*P1P2*S/(3*P1Q1*P2Q2)-64*A1*A2*MB**2*P1P2*S/(3*P1Q1*P2Q2)-
0642      &128*A1*A2*MB*MT*P1P2*S/(3*P1Q1*P2Q2)-
0643      &64*A1*A2*MT**2*P1P2*S/(3*P1Q1*P2Q2)-
0644      &256*A1*A2*P1P2**2*S/(3*P1Q1*P2Q2)-160*A2**2*P1Q1*S/(3*P2Q2)+
0645      &8*A1*A2*MB**3*MT*S/(3*P1Q2*P2Q2)+
0646      &16*A1*A2*MB**2*MT**2*S/(3*P1Q2*P2Q2)+
0647      &8*A1*A2*MB*MT**3*S/(3*P1Q2*P2Q2)-8*A1*P1P2*S/(3*P1Q2*P2Q2)-
0648      &8*A2*P1P2*S/(3*P1Q2*P2Q2)+8*A1*A2*MB**2*P1P2*S/(3*P1Q2*P2Q2)+
0649      &16*A1*A2*MB*MT*P1P2*S/(3*P1Q2*P2Q2)+
0650      &8*A1*A2*MT**2*P1P2*S/(3*P1Q2*P2Q2)+
0651      &32*A1*A2*P1P2**2*S/(3*P1Q2*P2Q2)+4*P1P2*S/(3*P1Q1*P1Q2*P2Q2)+
0652      &8*A1*MB*MT*P1P2*S/(3*P1Q1*P1Q2*P2Q2)-
0653      &8*A1*MT**2*P1P2*S/(3*P1Q1*P1Q2*P2Q2)+
0654      &16*A1*A2*P1P2*P1Q1*S/(3*P1Q2*P2Q2)-32*A2**2*P1Q2*S/(3*P2Q2)
0655       V18BIS=V18BIS+136*A2*P1Q2*S/(3*P1Q1*P2Q2)-
0656      &128*A1*A2*MB**2*P1Q2*S/(3*P1Q1*P2Q2)-
0657      &128*A1*A2*MB*MT*P1Q2*S/(3*P1Q1*P2Q2)-
0658      &256*A1*A2*P1P2*P1Q2*S/(3*P1Q1*P2Q2)-16*A2*MB*MT*S/(3*P2Q1*P2Q2)+
0659      &32*A2**2*MB**3*MT*S/(3*P2Q1*P2Q2)-16*A2*P1P2*S/(3*P2Q1*P2Q2)-
0660      &4*P1P2*S/(3*P1Q1*P2Q1*P2Q2)+8*A2*MB**2*P1P2*S/(3*P1Q1*P2Q1*P2Q2)-
0661      &8*A2*MB*MT*P1P2*S/(3*P1Q1*P2Q1*P2Q2)-4*P1P2*S/(3*P1Q2*P2Q1*P2Q2)+
0662      &8*A2*MB**2*P1P2*S/(3*P1Q2*P2Q1*P2Q2)-
0663      &8*A2*MB*MT*P1P2*S/(3*P1Q2*P2Q1*P2Q2)+
0664      &2*MB**3*MT*S/(3*P1Q1*P1Q2*P2Q1*P2Q2)+
0665      &4*MB**2*MT**2*S/(3*P1Q1*P1Q2*P2Q1*P2Q2)+
0666      &2*MB*MT**3*S/(3*P1Q1*P1Q2*P2Q1*P2Q2)-
0667      &2*MB**2*P1P2*S/(3*P1Q1*P1Q2*P2Q1*P2Q2)-
0668      &4*MB*MT*P1P2*S/(3*P1Q1*P1Q2*P2Q1*P2Q2)-
0669      &2*MT**2*P1P2*S/(3*P1Q1*P1Q2*P2Q1*P2Q2)-
0670      &8*P1P2**2*S/(3*P1Q1*P1Q2*P2Q1*P2Q2)+
0671      &8*A2*P1P2*P1Q1*S/(3*P1Q2*P2Q1*P2Q2)
0672       V18BIS=V18BIS+8*A2*P1P2*P1Q2*S/(3*P1Q1*P2Q1*P2Q2)+
0673      &272*A1*A2*P2Q1*S/(3*P2Q2)-
0674      &128*A1*MT**2*P2Q1*S/(3*P1Q1**2*P2Q2)-136*A1*P2Q1*S/(3*P1Q1*P2Q2)+
0675      &128*A1*A2*MB*MT*P2Q1*S/(3*P1Q1*P2Q2)+
0676      &128*A1*A2*MT**2*P2Q1*S/(3*P1Q1*P2Q2)+
0677      &256*A1*A2*P1P2*P2Q1*S/(3*P1Q1*P2Q2)-
0678      &16*A1*A2*P1P2*P2Q1*S/(3*P1Q2*P2Q2)+
0679      &8*A1*P1P2*P2Q1*S/(3*P1Q1*P1Q2*P2Q2)+
0680      &256*A1*A2*P1Q2*P2Q1*S/(3*P1Q1*P2Q2)-
0681      &128*A12*MT**2*P2Q2*S/(3*P1Q1**2)-160*A12*P2Q2*S/(3*P1Q1)-
0682      &32*A12*P2Q2*S/(3*P1Q2)+272*A1*A2*P2Q2*S/(3*P2Q1)-
0683      &16*A1*A2*P1P2*P2Q2*S/(3*P1Q1*P2Q1)-
0684      &128*A1*MT**2*P2Q2*S/(3*P1Q2**2*P2Q1)-136*A1*P2Q2*S/(3*P1Q2*P2Q1)+
0685      &128*A1*A2*MB*MT*P2Q2*S/(3*P1Q2*P2Q1)+
0686      &128*A1*A2*MT**2*P2Q2*S/(3*P1Q2*P2Q1)+
0687      &256*A1*A2*P1P2*P2Q2*S/(3*P1Q2*P2Q1)+
0688      &8*A1*P1P2*P2Q2*S/(3*P1Q1*P1Q2*P2Q1)
0689       V18BIS=V18BIS+256*A1*A2*P1Q1*P2Q2*S/(3*P1Q2*P2Q1)+
0690      &8*A12*MB*MT*S**2/(3*P1Q1*P1Q2)+16*A12*P1P2*S**2/(3*P1Q1*P1Q2)-
0691      &8*A1*A2*P1P2*S**2/(3*P1Q1*P2Q1)+4*A1*P1P2*S**2/(3*P1Q1*P1Q2*P2Q1)-
0692      &8*A1*A2*P1P2*S**2/(3*P1Q2*P2Q2)+4*A1*P1P2*S**2/(3*P1Q1*P1Q2*P2Q2)+
0693      &8*A2**2*MB*MT*S**2/(3*P2Q1*P2Q2)+16*A2**2*P1P2*S**2/(3*P2Q1*P2Q2)-
0694      &4*A2*P1P2*S**2/(3*P1Q1*P2Q1*P2Q2)-
0695      &4*A2*P1P2*S**2/(3*P1Q2*P2Q1*P2Q2)+
0696      &2*P1P2*S**2/(3*P1Q1*P1Q2*P2Q1*P2Q2)
0697 C
0698  
0699       A18 = 640*A1/3+640*A2/3+32*A1*A2*MB**2+368*A12*MB*MT+
0700      &512*A1*A2*MB*MT/3+
0701      &368*A2**2*MB*MT+32*A1*A2*MT**2+496*A12*P1P2/3+
0702      &320*A1*A2*P1P2+496*A2**2*P1P2/3-128*A1*MB*MT**3/(3*P1Q1**2)+
0703      &128*A1*MT**4/(3*P1Q1**2)+256*A12*MB*MT**5/(3*P1Q1**2)+
0704      &256*A1*MT**2*P1P2/(3*P1Q1**2)-256*A12*MT**4*P1P2/(3*P1Q1**2)+
0705      &8/(3*P1Q1)+32*A1*MB*MT/P1Q1+56*A2*MB*MT/(3*P1Q1)+
0706      &88*A1*MT**2/(3*P1Q1)+72*A2*MT**2/P1Q1-
0707      &704*A12*MB*MT**3/(3*P1Q1)+224*A1*A2*MB*MT**3/(3*P1Q1)+
0708      &104*A1*P1P2/(3*P1Q1)+48*A2*P1P2/P1Q1-
0709      &128*A1*A2*MB*MT*P1P2/(3*P1Q1)+512*A12*MT**2*P1P2/(3*P1Q1)-
0710      &448*A1*A2*MT**2*P1P2/(3*P1Q1)-32*A1*A2*P1P2**2/P1Q1-
0711      &656*A1*A2*P1Q1/3-224*A2**2*P1Q1-128*A1*MB*MT**3/(3*P1Q2**2)+
0712      &128*A1*MT**4/(3*P1Q2**2)+256*A12*MB*MT**5/(3*P1Q2**2)+
0713      &256*A1*MT**2*P1P2/(3*P1Q2**2)-256*A12*MT**4*P1P2/(3*P1Q2**2)+
0714      &256*A1*MT**2*P1Q1/(3*P1Q2**2)-256*A12*MB*MT**3*P1Q1/(3*P1Q2**2)+
0715      &8/(3*P1Q2)+32*A1*MB*MT/P1Q2+56*A2*MB*MT/(3*P1Q2)
0716       A18=A18+88*A1*MT**2/(3*P1Q2)+72*A2*MT**2/P1Q2-
0717      &704*A12*MB*MT**3/(3*P1Q2)+224*A1*A2*MB*MT**3/(3*P1Q2)+
0718      &104*A1*P1P2/(3*P1Q2)+48*A2*P1P2/P1Q2-
0719      &128*A1*A2*MB*MT*P1P2/(3*P1Q2)+512*A12*MT**2*P1P2/(3*P1Q2)-
0720      &448*A1*A2*MT**2*P1P2/(3*P1Q2)-32*A1*A2*P1P2**2/P1Q2+
0721      &32*A1*MB*MT**3/(3*P1Q1*P1Q2)-32*A1*MT**4/(3*P1Q1*P1Q2)-
0722      &64*A12*MB*MT**5/(3*P1Q1*P1Q2)+16*P1P2/(3*P1Q1*P1Q2)-
0723      &64*A1*MT**2*P1P2/(3*P1Q1*P1Q2)+64*A12*MT**4*P1P2/(3*P1Q1*P1Q2)+
0724      &112*A1*P1Q1/P1Q2+272*A2*P1Q1/(3*P1Q2)-
0725      &272*A1*A2*MB**2*P1Q1/(3*P1Q2)-208*A12*MB*MT*P1Q1/(3*P1Q2)+
0726      &400*A1*A2*MB*MT*P1Q1/(3*P1Q2)-80*A1*A2*MT**2*P1Q1/P1Q2+
0727      &96*A12*P1P2*P1Q1/P1Q2-320*A1*A2*P1P2*P1Q1/P1Q2-
0728      &544*A1*A2*P1Q1**2/(3*P1Q2)-656*A1*A2*P1Q2/3-224*A2**2*P1Q2+
0729      &256*A1*MT**2*P1Q2/(3*P1Q1**2)-256*A12*MB*MT**3*P1Q2/(3*P1Q1**2)+
0730      &112*A1*P1Q2/P1Q1+272*A2*P1Q2/(3*P1Q1)-
0731      &272*A1*A2*MB**2*P1Q2/(3*P1Q1)-208*A12*MB*MT*P1Q2/(3*P1Q1)+
0732      &400*A1*A2*MB*MT*P1Q2/(3*P1Q1)-80*A1*A2*MT**2*P1Q2/P1Q1
0733       A18=A18+96*A12*P1P2*P1Q2/P1Q1-320*A1*A2*P1P2*P1Q2/P1Q1-
0734      &544*A1*A2*P1Q2**2/(3*P1Q1)+128*A2*MB**4/(3*P2Q1**2)-
0735      &128*A2*MB**3*MT/(3*P2Q1**2)+256*A2**2*MB**5*MT/(3*P2Q1**2)+
0736      &256*A2*MB**2*P1P2/(3*P2Q1**2)-256*A2**2*MB**4*P1P2/(3*P2Q1**2)+
0737      &256*A2*MB**2*P1Q1/(3*P2Q1**2)-256*A2**2*MB**4*P1Q1/(3*P2Q1**2)+
0738      &64*MB**3*MT**3/(3*P1Q2**2*P2Q1**2)-
0739      &64*MB**2*MT**2*P1P2/(3*P1Q2**2*P2Q1**2)-
0740      &64*MB**2*MT**2*P1Q1/(3*P1Q2**2*P2Q1**2)-
0741      &64*MB**3*MT/(3*P1Q2*P2Q1**2)-
0742      &256*A2*MB**3*MT*P1P2/(3*P1Q2*P2Q1**2)+
0743      &256*A2*MB**2*P1P2**2/(3*P1Q2*P2Q1**2)-
0744      &256*A2*MB**3*MT*P1Q1/(3*P1Q2*P2Q1**2)+
0745      &512*A2*MB**2*P1P2*P1Q1/(3*P1Q2*P2Q1**2)+
0746      &256*A2*MB**2*P1Q1**2/(3*P1Q2*P2Q1**2)-
0747      &256*A2**2*MB**4*P1Q2/(3*P2Q1**2)-8/(3*P2Q1)-72*A1*MB**2/P2Q1-
0748      &88*A2*MB**2/(3*P2Q1)-56*A1*MB*MT/(3*P2Q1)-32*A2*MB*MT/P2Q1-
0749      &224*A1*A2*MB**3*MT/(3*P2Q1)+704*A2**2*MB**3*MT/(3*P2Q1)
0750       A18=A18-48*A1*P1P2/P2Q1-104*A2*P1P2/(3*P2Q1)+
0751      &448*A1*A2*MB**2*P1P2/(3*P2Q1)-512*A2**2*MB**2*P1P2/(3*P2Q1)+
0752      &128*A1*A2*MB*MT*P1P2/(3*P2Q1)+32*A1*A2*P1P2**2/P2Q1-
0753      &16*P1P2/(3*P1Q1*P2Q1)+32*A1*MB*MT*P1P2/(3*P1Q1*P2Q1)+
0754      &32*A2*MB*MT*P1P2/(3*P1Q1*P2Q1)+
0755      &64*A1*A2*MB*MT*P1P2**2/(3*P1Q1*P2Q1)-
0756      &64*A1*A2*P1P2**3/(3*P1Q1*P2Q1)-256*A2*P1Q1/(3*P2Q1)+
0757      &448*A1*A2*MB**2*P1Q1/(3*P2Q1)-368*A2**2*MB**2*P1Q1/(3*P2Q1)-
0758      &224*A1*A2*MB*MT*P1Q1/(3*P2Q1)+304*A1*A2*P1P2*P1Q1/(3*P2Q1)+
0759      &64*MB*MT**3/(3*P1Q2**2*P2Q1)+
0760      &256*A1*MB*MT**3*P1P2/(3*P1Q2**2*P2Q1)-
0761      &256*A1*MT**2*P1P2**2/(3*P1Q2**2*P2Q1)+
0762      &64*MT**2*P1Q1/(3*P1Q2**2*P2Q1)-
0763      &128*A1*MB**2*MT**2*P1Q1/(3*P1Q2**2*P2Q1)+
0764      &128*A1*MB*MT**3*P1Q1/(3*P1Q2**2*P2Q1)-
0765      &256*A1*MT**2*P1P2*P1Q1/(3*P1Q2**2*P2Q1)-4*MB**2/(3*P1Q2*P2Q1)-
0766      &64*MB*MT/(3*P1Q2*P2Q1)+128*A2*MB**3*MT/(3*P1Q2*P2Q1)
0767       A18=A18-4*MT**2/(3*P1Q2*P2Q1)-128*A1*MB**2*MT**2/(3*P1Q2*P2Q1)-
0768      &128*A2*MB**2*MT**2/(3*P1Q2*P2Q1)+128*A1*MB*MT**3/(3*P1Q2*P2Q1)-
0769      &112*A2*MB**2*P1P2/(3*P1Q2*P2Q1)+32*A1*MB*MT*P1P2/(3*P1Q2*P2Q1)+
0770      &32*A2*MB*MT*P1P2/(3*P1Q2*P2Q1)-112*A1*MT**2*P1P2/(3*P1Q2*P2Q1)-
0771      &48*A1*P1P2**2/(P1Q2*P2Q1)-48*A2*P1P2**2/(P1Q2*P2Q1)-
0772      &512*A1*A2*MB*MT*P1P2**2/(3*P1Q2*P2Q1)+
0773      &512*A1*A2*P1P2**3/(3*P1Q2*P2Q1)+8*MB*MT*P1P2/(3*P1Q1*P1Q2*P2Q1)-
0774      &8*MT**2*P1P2/(3*P1Q1*P1Q2*P2Q1)-
0775      &32*A1*MB*MT**3*P1P2/(3*P1Q1*P1Q2*P2Q1)-
0776      &16*P1P2**2/(3*P1Q1*P1Q2*P2Q1)+
0777      &32*A1*MT**2*P1P2**2/(3*P1Q1*P1Q2*P2Q1)+8*P1Q1/(3*P1Q2*P2Q1)-
0778      &160*A1*MB**2*P1Q1/(3*P1Q2*P2Q1)-272*A2*MB**2*P1Q1/(3*P1Q2*P2Q1)-
0779      &56*A1*MB*MT*P1Q1/(3*P1Q2*P2Q1)-200*A2*MB*MT*P1Q1/(3*P1Q2*P2Q1)-
0780      &48*A1*P1P2*P1Q1/(P1Q2*P2Q1)-256*A2*P1P2*P1Q1/(3*P1Q2*P2Q1)+
0781      &256*A1*A2*MB**2*P1P2*P1Q1/(3*P1Q2*P2Q1)-
0782      &256*A1*A2*MB*MT*P1P2*P1Q1/(P1Q2*P2Q1)+
0783      &1024*A1*A2*P1P2**2*P1Q1/(3*P1Q2*P2Q1)
0784       A18=A18-272*A2*P1Q1**2/(3*P1Q2*P2Q1)+
0785      &256*A1*A2*MB**2*P1Q1**2/(3*P1Q2*P2Q1)-
0786      &256*A1*A2*MB*MT*P1Q1**2/(3*P1Q2*P2Q1)+
0787      &512*A1*A2*P1P2*P1Q1**2/(3*P1Q2*P2Q1)+16*A2*P1Q2/(3*P2Q1)+
0788      &64*A1*A2*MB**2*P1Q2/P2Q1+32*A2**2*MB**2*P1Q2/(3*P2Q1)-
0789      &112*A1*A2*MB*MT*P1Q2/(3*P2Q1)+368*A1*A2*P1P2*P1Q2/(3*P2Q1)+
0790      &32*A2*P1P2*P1Q2/(3*P1Q1*P2Q1)-
0791      &32*A1*A2*MB**2*P1P2*P1Q2/(3*P1Q1*P2Q1)+
0792      &32*A1*A2*MB*MT*P1P2*P1Q2/(3*P1Q1*P2Q1)-
0793      &64*A1*A2*P1P2**2*P1Q2/(3*P1Q1*P2Q1)+224*A12*P2Q1+
0794      &656*A1*A2*P2Q1/3-256*A1*MT**2*P2Q1/(3*P1Q1**2)+
0795      &256*A12*MT**4*P2Q1/(3*P1Q1**2)-256*A1*P2Q1/(3*P1Q1)-
0796      &224*A1*A2*MB*MT*P2Q1/(3*P1Q1)-368*A12*MT**2*P2Q1/(3*P1Q1)+
0797      &448*A1*A2*MT**2*P2Q1/(3*P1Q1)+304*A1*A2*P1P2*P2Q1/(3*P1Q1)+
0798      &256*A12*MT**4*P2Q1/(3*P1Q2**2)+
0799      &256*A12*MT**2*P1Q1*P2Q1/(3*P1Q2**2)+16*A1*P2Q1/(3*P1Q2)-
0800      &112*A1*A2*MB*MT*P2Q1/(3*P1Q2)+32*A12*MT**2*P2Q1/(3*P1Q2)
0801       A18=A18+64*A1*A2*MT**2*P2Q1/P1Q2+368*A1*A2*P1P2*P2Q1/(3*P1Q2)+
0802      &16*A1*MT**2*P2Q1/(3*P1Q1*P1Q2)-64*A12*MT**4*P2Q1/(3*P1Q1*P1Q2)+
0803      &640*A12*P1Q1*P2Q1/(3*P1Q2)+544*A1*A2*P1Q1*P2Q1/(3*P1Q2)+
0804      &32*A12*P1Q2*P2Q1/P1Q1+944*A1*A2*P1Q2*P2Q1/(3*P1Q1)+
0805      &128*A2*MB**4/(3*P2Q2**2)-128*A2*MB**3*MT/(3*P2Q2**2)+
0806      &256*A2**2*MB**5*MT/(3*P2Q2**2)+256*A2*MB**2*P1P2/(3*P2Q2**2)-
0807      &256*A2**2*MB**4*P1P2/(3*P2Q2**2)+
0808      &64*MB**3*MT**3/(3*P1Q1**2*P2Q2**2)-
0809      &64*MB**2*MT**2*P1P2/(3*P1Q1**2*P2Q2**2)-
0810      &64*MB**3*MT/(3*P1Q1*P2Q2**2)-
0811      &256*A2*MB**3*MT*P1P2/(3*P1Q1*P2Q2**2)+
0812      &256*A2*MB**2*P1P2**2/(3*P1Q1*P2Q2**2)-
0813      &256*A2**2*MB**4*P1Q1/(3*P2Q2**2)+256*A2*MB**2*P1Q2/(3*P2Q2**2)-
0814      &256*A2**2*MB**4*P1Q2/(3*P2Q2**2)-
0815      &64*MB**2*MT**2*P1Q2/(3*P1Q1**2*P2Q2**2)-
0816      &256*A2*MB**3*MT*P1Q2/(3*P1Q1*P2Q2**2)+
0817      &512*A2*MB**2*P1P2*P1Q2/(3*P1Q1*P2Q2**2)
0818       A18=A18+256*A2*MB**2*P1Q2**2/(3*P1Q1*P2Q2**2)-
0819      &256*A2*MB**2*P2Q1/(3*P2Q2**2)+256*A2**2*MB**3*MT*P2Q1/(3*P2Q2**2)+
0820      &64*MB**2*MT**2*P2Q1/(3*P1Q1**2*P2Q2**2)+
0821      &64*MB**2*P2Q1/(3*P1Q1*P2Q2**2)+
0822      &128*A2*MB**3*MT*P2Q1/(3*P1Q1*P2Q2**2)-
0823      &128*A2*MB**2*MT**2*P2Q1/(3*P1Q1*P2Q2**2)-
0824      &256*A2*MB**2*P1P2*P2Q1/(3*P1Q1*P2Q2**2)+
0825      &256*A2**2*MB**2*P1Q1*P2Q1/(3*P2Q2**2)-
0826      &256*A2*MB**2*P1Q2*P2Q1/(3*P1Q1*P2Q2**2)-8/(3*P2Q2)-
0827      &72*A1*MB**2/P2Q2-88*A2*MB**2/(3*P2Q2)-56*A1*MB*MT/(3*P2Q2)-
0828      &32*A2*MB*MT/P2Q2-224*A1*A2*MB**3*MT/(3*P2Q2)+
0829      &704*A2**2*MB**3*MT/(3*P2Q2)-48*A1*P1P2/P2Q2-
0830      &104*A2*P1P2/(3*P2Q2)+448*A1*A2*MB**2*P1P2/(3*P2Q2)-
0831      &512*A2**2*MB**2*P1P2/(3*P2Q2)+128*A1*A2*MB*MT*P1P2/(3*P2Q2)+
0832      &32*A1*A2*P1P2**2/P2Q2+64*MB*MT**3/(3*P1Q1**2*P2Q2)+
0833      &256*A1*MB*MT**3*P1P2/(3*P1Q1**2*P2Q2)-
0834      &256*A1*MT**2*P1P2**2/(3*P1Q1**2*P2Q2)-4*MB**2/(3*P1Q1*P2Q2)
0835       A18=A18-64*MB*MT/(3*P1Q1*P2Q2)+128*A2*MB**3*MT/(3*P1Q1*P2Q2)-
0836      &4*MT**2/(3*P1Q1*P2Q2)-128*A1*MB**2*MT**2/(3*P1Q1*P2Q2)-
0837      &128*A2*MB**2*MT**2/(3*P1Q1*P2Q2)+128*A1*MB*MT**3/(3*P1Q1*P2Q2)-
0838      &112*A2*MB**2*P1P2/(3*P1Q1*P2Q2)+32*A1*MB*MT*P1P2/(3*P1Q1*P2Q2)+
0839      &32*A2*MB*MT*P1P2/(3*P1Q1*P2Q2)-112*A1*MT**2*P1P2/(3*P1Q1*P2Q2)-
0840      &48*A1*P1P2**2/(P1Q1*P2Q2)-48*A2*P1P2**2/(P1Q1*P2Q2)-
0841      &512*A1*A2*MB*MT*P1P2**2/(3*P1Q1*P2Q2)+
0842      &512*A1*A2*P1P2**3/(3*P1Q1*P2Q2)+16*A2*P1Q1/(3*P2Q2)+
0843      &64*A1*A2*MB**2*P1Q1/P2Q2+32*A2**2*MB**2*P1Q1/(3*P2Q2)-
0844      &112*A1*A2*MB*MT*P1Q1/(3*P2Q2)+368*A1*A2*P1P2*P1Q1/(3*P2Q2)-
0845      &16*P1P2/(3*P1Q2*P2Q2)+32*A1*MB*MT*P1P2/(3*P1Q2*P2Q2)+
0846      &32*A2*MB*MT*P1P2/(3*P1Q2*P2Q2)+
0847      &64*A1*A2*MB*MT*P1P2**2/(3*P1Q2*P2Q2)-
0848      &64*A1*A2*P1P2**3/(3*P1Q2*P2Q2)+8*MB*MT*P1P2/(3*P1Q1*P1Q2*P2Q2)-
0849      &8*MT**2*P1P2/(3*P1Q1*P1Q2*P2Q2)-
0850      &32*A1*MB*MT**3*P1P2/(3*P1Q1*P1Q2*P2Q2)-
0851      &16*P1P2**2/(3*P1Q1*P1Q2*P2Q2)
0852       A18=A18+32*A1*MT**2*P1P2**2/(3*P1Q1*P1Q2*P2Q2)+
0853      &32*A2*P1P2*P1Q1/(3*P1Q2*P2Q2)-
0854      &32*A1*A2*MB**2*P1P2*P1Q1/(3*P1Q2*P2Q2)+
0855      &32*A1*A2*MB*MT*P1P2*P1Q1/(3*P1Q2*P2Q2)-
0856      &64*A1*A2*P1P2**2*P1Q1/(3*P1Q2*P2Q2)-256*A2*P1Q2/(3*P2Q2)+
0857      &448*A1*A2*MB**2*P1Q2/(3*P2Q2)-368*A2**2*MB**2*P1Q2/(3*P2Q2)-
0858      &224*A1*A2*MB*MT*P1Q2/(3*P2Q2)+304*A1*A2*P1P2*P1Q2/(3*P2Q2)+
0859      &64*MT**2*P1Q2/(3*P1Q1**2*P2Q2)-
0860      &128*A1*MB**2*MT**2*P1Q2/(3*P1Q1**2*P2Q2)+
0861      &128*A1*MB*MT**3*P1Q2/(3*P1Q1**2*P2Q2)-
0862      &256*A1*MT**2*P1P2*P1Q2/(3*P1Q1**2*P2Q2)+8*P1Q2/(3*P1Q1*P2Q2)-
0863      &160*A1*MB**2*P1Q2/(3*P1Q1*P2Q2)-272*A2*MB**2*P1Q2/(3*P1Q1*P2Q2)-
0864      &56*A1*MB*MT*P1Q2/(3*P1Q1*P2Q2)-200*A2*MB*MT*P1Q2/(3*P1Q1*P2Q2)-
0865      &48*A1*P1P2*P1Q2/(P1Q1*P2Q2)-256*A2*P1P2*P1Q2/(3*P1Q1*P2Q2)+
0866      &256*A1*A2*MB**2*P1P2*P1Q2/(3*P1Q1*P2Q2)-
0867      &256*A1*A2*MB*MT*P1P2*P1Q2/(P1Q1*P2Q2)+
0868      &1024*A1*A2*P1P2**2*P1Q2/(3*P1Q1*P2Q2)
0869       A18=A18-272*A2*P1Q2**2/(3*P1Q1*P2Q2)+
0870      &256*A1*A2*MB**2*P1Q2**2/(3*P1Q1*P2Q2)-
0871      &256*A1*A2*MB*MT*P1Q2**2/(3*P1Q1*P2Q2)+
0872      &512*A1*A2*P1P2*P1Q2**2/(3*P1Q1*P2Q2)-32*A2*MB**4/(3*P2Q1*P2Q2)+
0873      &32*A2*MB**3*MT/(3*P2Q1*P2Q2)-64*A2**2*MB**5*MT/(3*P2Q1*P2Q2)+
0874      &16*P1P2/(3*P2Q1*P2Q2)-64*A2*MB**2*P1P2/(3*P2Q1*P2Q2)+
0875      &64*A2**2*MB**4*P1P2/(3*P2Q1*P2Q2)+8*MB**2*P1P2/(3*P1Q1*P2Q1*P2Q2)-
0876      &8*MB*MT*P1P2/(3*P1Q1*P2Q1*P2Q2)+
0877      &32*A2*MB**3*MT*P1P2/(3*P1Q1*P2Q1*P2Q2)+
0878      &16*P1P2**2/(3*P1Q1*P2Q1*P2Q2)-
0879      &32*A2*MB**2*P1P2**2/(3*P1Q1*P2Q1*P2Q2)-
0880      &16*A2*MB**2*P1Q1/(3*P2Q1*P2Q2)+64*A2**2*MB**4*P1Q1/(3*P2Q1*P2Q2)+
0881      &8*MB**2*P1P2/(3*P1Q2*P2Q1*P2Q2)-8*MB*MT*P1P2/(3*P1Q2*P2Q1*P2Q2)+
0882      &32*A2*MB**3*MT*P1P2/(3*P1Q2*P2Q1*P2Q2)+
0883      &16*P1P2**2/(3*P1Q2*P2Q1*P2Q2)-
0884      &32*A2*MB**2*P1P2**2/(3*P1Q2*P2Q1*P2Q2)-
0885      &16*MB*MT*P1P2**2/(3*P1Q1*P1Q2*P2Q1*P2Q2)
0886       A18=A18+16*P1P2**3/(3*P1Q1*P1Q2*P2Q1*P2Q2)-
0887      &32*A2*MB**2*P1P2*P1Q1/(3*P1Q2*P2Q1*P2Q2)-
0888      &16*A2*MB**2*P1Q2/(3*P2Q1*P2Q2)+64*A2**2*MB**4*P1Q2/(3*P2Q1*P2Q2)-
0889      &32*A2*MB**2*P1P2*P1Q2/(3*P1Q1*P2Q1*P2Q2)+272*A1*P2Q1/(3*P2Q2)+
0890      &112*A2*P2Q1/P2Q2-80*A1*A2*MB**2*P2Q1/P2Q2+
0891      &400*A1*A2*MB*MT*P2Q1/(3*P2Q2)-208*A2**2*MB*MT*P2Q1/(3*P2Q2)-
0892      &272*A1*A2*MT**2*P2Q1/(3*P2Q2)-320*A1*A2*P1P2*P2Q1/P2Q2+
0893      &96*A2**2*P1P2*P2Q1/P2Q2-256*A1*MB*MT**3*P2Q1/(3*P1Q1**2*P2Q2)+
0894      &512*A1*MT**2*P1P2*P2Q1/(3*P1Q1**2*P2Q2)-8*P2Q1/(3*P1Q1*P2Q2)+
0895      &200*A1*MB*MT*P2Q1/(3*P1Q1*P2Q2)+56*A2*MB*MT*P2Q1/(3*P1Q1*P2Q2)+
0896      &272*A1*MT**2*P2Q1/(3*P1Q1*P2Q2)+160*A2*MT**2*P2Q1/(3*P1Q1*P2Q2)+
0897      &256*A1*P1P2*P2Q1/(3*P1Q1*P2Q2)+48*A2*P1P2*P2Q1/(P1Q1*P2Q2)+
0898      &256*A1*A2*MB*MT*P1P2*P2Q1/(P1Q1*P2Q2)-
0899      &256*A1*A2*MT**2*P1P2*P2Q1/(3*P1Q1*P2Q2)-
0900      &1024*A1*A2*P1P2**2*P2Q1/(3*P1Q1*P2Q2)-
0901      &544*A1*A2*P1Q1*P2Q1/(3*P2Q2)-640*A2**2*P1Q1*P2Q1/(3*P2Q2)-
0902      &32*A1*P1P2*P2Q1/(3*P1Q2*P2Q2)
0903       A18=A18-32*A1*A2*MB*MT*P1P2*P2Q1/(3*P1Q2*P2Q2)+
0904      &32*A1*A2*MT**2*P1P2*P2Q1/(3*P1Q2*P2Q2)+
0905      &64*A1*A2*P1P2**2*P2Q1/(3*P1Q2*P2Q2)-
0906      &32*A1*MT**2*P1P2*P2Q1/(3*P1Q1*P1Q2*P2Q2)+
0907      &64*A1*A2*P1P2*P1Q1*P2Q1/(3*P1Q2*P2Q2)-
0908      &944*A1*A2*P1Q2*P2Q1/(3*P2Q2)-32*A2**2*P1Q2*P2Q1/P2Q2+
0909      &256*A1*MT**2*P1Q2*P2Q1/(3*P1Q1**2*P2Q2)+
0910      &96*A1*P1Q2*P2Q1/(P1Q1*P2Q2)+96*A2*P1Q2*P2Q1/(P1Q1*P2Q2)-
0911      &128*A1*A2*MB**2*P1Q2*P2Q1/(3*P1Q1*P2Q2)+
0912      &256*A1*A2*MB*MT*P1Q2*P2Q1/(P1Q1*P2Q2)-
0913      &128*A1*A2*MT**2*P1Q2*P2Q1/(3*P1Q1*P2Q2)-
0914      &512*A1*A2*P1P2*P1Q2*P2Q1/(P1Q1*P2Q2)-
0915      &512*A1*A2*P1Q2**2*P2Q1/(3*P1Q1*P2Q2)+544*A1*A2*P2Q1**2/(3*P2Q2)-
0916      &256*A1*MT**2*P2Q1**2/(3*P1Q1**2*P2Q2)-
0917      &272*A1*P2Q1**2/(3*P1Q1*P2Q2)-
0918      &256*A1*A2*MB*MT*P2Q1**2/(3*P1Q1*P2Q2)+
0919      &256*A1*A2*MT**2*P2Q1**2/(3*P1Q1*P2Q2)
0920       A18=A18+512*A1*A2*P1P2*P2Q1**2/(3*P1Q1*P2Q2)+
0921      &512*A1*A2*P1Q2*P2Q1**2/(3*P1Q1*P2Q2)+224*A12*P2Q2+
0922      &656*A1*A2*P2Q2/3+256*A12*MT**4*P2Q2/(3*P1Q1**2)+
0923      &16*A1*P2Q2/(3*P1Q1)-112*A1*A2*MB*MT*P2Q2/(3*P1Q1)+
0924      &32*A12*MT**2*P2Q2/(3*P1Q1)+64*A1*A2*MT**2*P2Q2/P1Q1+
0925      &368*A1*A2*P1P2*P2Q2/(3*P1Q1)-256*A1*MT**2*P2Q2/(3*P1Q2**2)+
0926      &256*A12*MT**4*P2Q2/(3*P1Q2**2)-256*A1*P2Q2/(3*P1Q2)-
0927      &224*A1*A2*MB*MT*P2Q2/(3*P1Q2)-368*A12*MT**2*P2Q2/(3*P1Q2)+
0928      &448*A1*A2*MT**2*P2Q2/(3*P1Q2)+304*A1*A2*P1P2*P2Q2/(3*P1Q2)+
0929      &16*A1*MT**2*P2Q2/(3*P1Q1*P1Q2)-64*A12*MT**4*P2Q2/(3*P1Q1*P1Q2)+
0930      &32*A12*P1Q1*P2Q2/P1Q2+944*A1*A2*P1Q1*P2Q2/(3*P1Q2)+
0931      &256*A12*MT**2*P1Q2*P2Q2/(3*P1Q1**2)+
0932      &640*A12*P1Q2*P2Q2/(3*P1Q1)+544*A1*A2*P1Q2*P2Q2/(3*P1Q1)-
0933      &256*A2*MB**2*P2Q2/(3*P2Q1**2)+256*A2**2*MB**3*MT*P2Q2/(3*P2Q1**2)+
0934      &64*MB**2*MT**2*P2Q2/(3*P1Q2**2*P2Q1**2)+
0935      &64*MB**2*P2Q2/(3*P1Q2*P2Q1**2)+
0936      &128*A2*MB**3*MT*P2Q2/(3*P1Q2*P2Q1**2)
0937       A18=A18-128*A2*MB**2*MT**2*P2Q2/(3*P1Q2*P2Q1**2)-
0938      &256*A2*MB**2*P1P2*P2Q2/(3*P1Q2*P2Q1**2)-
0939      &256*A2*MB**2*P1Q1*P2Q2/(3*P1Q2*P2Q1**2)+
0940      &256*A2**2*MB**2*P1Q2*P2Q2/(3*P2Q1**2)+272*A1*P2Q2/(3*P2Q1)+
0941      &112*A2*P2Q2/P2Q1-80*A1*A2*MB**2*P2Q2/P2Q1+
0942      &400*A1*A2*MB*MT*P2Q2/(3*P2Q1)-208*A2**2*MB*MT*P2Q2/(3*P2Q1)-
0943      &272*A1*A2*MT**2*P2Q2/(3*P2Q1)-320*A1*A2*P1P2*P2Q2/P2Q1+
0944      &96*A2**2*P1P2*P2Q2/P2Q1-32*A1*P1P2*P2Q2/(3*P1Q1*P2Q1)-
0945      &32*A1*A2*MB*MT*P1P2*P2Q2/(3*P1Q1*P2Q1)+
0946      &32*A1*A2*MT**2*P1P2*P2Q2/(3*P1Q1*P2Q1)+
0947      &64*A1*A2*P1P2**2*P2Q2/(3*P1Q1*P2Q1)-944*A1*A2*P1Q1*P2Q2/(3*P2Q1)-
0948      &32*A2**2*P1Q1*P2Q2/P2Q1-256*A1*MB*MT**3*P2Q2/(3*P1Q2**2*P2Q1)+
0949      &512*A1*MT**2*P1P2*P2Q2/(3*P1Q2**2*P2Q1)+
0950      &256*A1*MT**2*P1Q1*P2Q2/(3*P1Q2**2*P2Q1)-8*P2Q2/(3*P1Q2*P2Q1)+
0951      &200*A1*MB*MT*P2Q2/(3*P1Q2*P2Q1)+56*A2*MB*MT*P2Q2/(3*P1Q2*P2Q1)+
0952      &272*A1*MT**2*P2Q2/(3*P1Q2*P2Q1)+160*A2*MT**2*P2Q2/(3*P1Q2*P2Q1)+
0953      &256*A1*P1P2*P2Q2/(3*P1Q2*P2Q1)+48*A2*P1P2*P2Q2/(P1Q2*P2Q1)
0954       A18=A18+256*A1*A2*MB*MT*P1P2*P2Q2/(P1Q2*P2Q1)-
0955      &256*A1*A2*MT**2*P1P2*P2Q2/(3*P1Q2*P2Q1)-
0956      &1024*A1*A2*P1P2**2*P2Q2/(3*P1Q2*P2Q1)-
0957      &32*A1*MT**2*P1P2*P2Q2/(3*P1Q1*P1Q2*P2Q1)+
0958      &96*A1*P1Q1*P2Q2/(P1Q2*P2Q1)+96*A2*P1Q1*P2Q2/(P1Q2*P2Q1)-
0959      &128*A1*A2*MB**2*P1Q1*P2Q2/(3*P1Q2*P2Q1)+
0960      &256*A1*A2*MB*MT*P1Q1*P2Q2/(P1Q2*P2Q1)-
0961      &128*A1*A2*MT**2*P1Q1*P2Q2/(3*P1Q2*P2Q1)-
0962      &512*A1*A2*P1P2*P1Q1*P2Q2/(P1Q2*P2Q1)-
0963      &512*A1*A2*P1Q1**2*P2Q2/(3*P1Q2*P2Q1)-544*A1*A2*P1Q2*P2Q2/(3*P2Q1)-
0964      &640*A2**2*P1Q2*P2Q2/(3*P2Q1)+
0965      &64*A1*A2*P1P2*P1Q2*P2Q2/(3*P1Q1*P2Q1)+544*A1*A2*P2Q2**2/(3*P2Q1)-
0966      &256*A1*MT**2*P2Q2**2/(3*P1Q2**2*P2Q1)-
0967      &272*A1*P2Q2**2/(3*P1Q2*P2Q1)-
0968      &256*A1*A2*MB*MT*P2Q2**2/(3*P1Q2*P2Q1)+
0969      &256*A1*A2*MT**2*P2Q2**2/(3*P1Q2*P2Q1)+
0970      &512*A1*A2*P1P2*P2Q2**2/(3*P1Q2*P2Q1)
0971       A18=A18+512*A1*A2*P1Q1*P2Q2**2/(3*P1Q2*P2Q1)-
0972      &384*A12*MB*MT*P1Q1**2/S**2+
0973      &384*A12*P1P2*P1Q1**2/S**2-2688*A12*MB*MT*P1Q1*P1Q2/S**2+
0974      &2688*A12*P1P2*P1Q1*P1Q2/S**2-384*A12*MB*MT*P1Q2**2/S**2+
0975      &384*A12*P1P2*P1Q2**2/S**2-768*A1*A2*MB*MT*P1Q1*P2Q1/S**2+
0976      &768*A1*A2*P1P2*P1Q1*P2Q1/S**2-2688*A1*A2*MB*MT*P1Q2*P2Q1/S**2+
0977      &2688*A1*A2*P1P2*P1Q2*P2Q1/S**2-960*A12*P1Q1*P1Q2*P2Q1/S**2-
0978      &960*A1*A2*P1Q1*P1Q2*P2Q1/S**2+960*A12*P1Q2**2*P2Q1/S**2+
0979      &960*A1*A2*P1Q2**2*P2Q1/S**2-384*A2**2*MB*MT*P2Q1**2/S**2+
0980      &384*A2**2*P1P2*P2Q1**2/S**2-960*A1*A2*P1Q2*P2Q1**2/S**2-
0981      &960*A2**2*P1Q2*P2Q1**2/S**2-2688*A1*A2*MB*MT*P1Q1*P2Q2/S**2+
0982      &2688*A1*A2*P1P2*P1Q1*P2Q2/S**2+960*A12*P1Q1**2*P2Q2/S**2+
0983      &960*A1*A2*P1Q1**2*P2Q2/S**2-768*A1*A2*MB*MT*P1Q2*P2Q2/S**2+
0984      &768*A1*A2*P1P2*P1Q2*P2Q2/S**2-960*A12*P1Q1*P1Q2*P2Q2/S**2-
0985      &960*A1*A2*P1Q1*P1Q2*P2Q2/S**2-2688*A2**2*MB*MT*P2Q1*P2Q2/S**2+
0986      &2688*A2**2*P1P2*P2Q1*P2Q2/S**2+960*A1*A2*P1Q1*P2Q1*P2Q2/S**2+
0987      &960*A2**2*P1Q1*P2Q1*P2Q2/S**2+960*A1*A2*P1Q2*P2Q1*P2Q2/S**2
0988       A18=A18+960*A2**2*P1Q2*P2Q1*P2Q2/S**2-
0989      &384*A2**2*MB*MT*P2Q2**2/S**2+
0990      &384*A2**2*P1P2*P2Q2**2/S**2-960*A1*A2*P1Q1*P2Q2**2/S**2-
0991      &960*A2**2*P1Q1*P2Q2**2/S**2-96*A1*MB*MT/S-96*A2*MB*MT/S+
0992      &768*A2**2*MB**3*MT/S+768*A12*MB*MT**3/S-192*A1*P1P2/S-
0993      &192*A2*P1P2/S-768*A2**2*MB**2*P1P2/S+2304*A1*A2*MB*MT*P1P2/S-
0994      &768*A12*MT**2*P1P2/S-2304*A1*A2*P1P2**2/S+
0995      &96*A1*MB*MT**3/(P1Q1*S)+192*A2*MB*MT*P1P2/(P1Q1*S)-
0996      &96*A1*MT**2*P1P2/(P1Q1*S)-192*A2*P1P2**2/(P1Q1*S)-192*A1*P1Q1/S-
0997      &144*A2*P1Q1/S-384*A1*A2*MB**2*P1Q1/S-480*A2**2*MB**2*P1Q1/S+
0998      &480*A12*MB*MT*P1Q1/S-96*A1*A2*MB*MT*P1Q1/S-
0999      &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
[ Source navigation ] [ Diff markup ] [ Identifier search ] [ general search ]

This page was automatically generated by the 2.3.7 LXR engine.

The LXR team
Valid CSS 2.1! Valid HTML 4.01!