/
dr3.in
10 lines (9 loc) · 41.8 KB
/
dr3.in
1
2
3
4
5
6
7
8
9
double dr_SMNNLO::res3fff(void) const
{
double MZs = sqr(ival->get(MZ)),
MWs = sqr(ival->get(MW)),
MTs = sqr(ival->get(MT)),
AL = ival->get(al),
deltaAlpha = ival->get(Delal);
return(real(-(AL*sqr(deltaAlpha)*(128*MWs*(MWs - MZs)*MZs + 3*MTs*(-32*MWs*MZs + 128*sqr(MWs) - 9*sqr(MZs))))/(144.*MWs*(MWs - MZs)*MZs*Pi) + trip(deltaAlpha) + (deltaAlpha*powint(MWs - MZs,-3)*sqr(AL)*(768*MTs*MWs*MZs*sqr(MWs - MZs)*(-32*MWs*MZs + 128*sqr(MWs) - 9*sqr(MZs)) + 9*sqr(MTs)*(-24576*MZs*powint(MWs,4) + 16384*powint(MWs,5) - 495*MWs*powint(MZs,4) - 81*powint(MZs,5) + 9216*sqr(MZs)*trip(MWs) - 448*sqr(MWs)*trip(MZs)) + 4*MWs*sqr(MZs)*(4096*powint(MWs,4) + 95481*powint(MZs,4)*sqr(Pi) + 24*sqr(MWs)*sqr(MZs)*(512 + 9423*sqr(Pi)) - 768*MZs*(16 + 75*sqr(Pi))*trip(MWs) - 4*MWs*(1024 + 63639*sqr(Pi))*trip(MZs))))/(62208.*sqr(MWs)*sqr(MZs)*sqr(Pi)) + sqr(real(B0fin(MWs,0,0)))*((27*deltaAlpha*powint(MWs - MZs,-4)*sqr(AL)*sqr(MZs)*(-4*MWs*MZs + 5*sqr(MWs) + sqr(MZs)))/(16.*sqr(Pi)) - (27*(MTs - MWs)*(MTs + 2*MWs)*powint(MWs - MZs,-6)*powint(Pi,-3)*real(B0fin(MWs,MTs,0))*trip(AL)*(-14*MZs*sqr(MWs) + 6*MWs*sqr(MZs) + 14*trip(MWs) - trip(MZs))*trip(MZs))/(128.*sqr(MWs)) + (3*powint(MWs - MZs,-6)*powint(Pi,-3)*real(B0fin(MZs,0,0))*trip(AL)*trip(MZs)*(812*MZs*sqr(MWs) - 344*MWs*sqr(MZs) - 840*trip(MWs) + 57*trip(MZs)))/128. - (3*MZs*powint(MWs - MZs,-6)*powint(Pi,-3)*real(B0fin(MZs,MTs,MTs))*trip(AL)*(MZs*(160*powint(MWs,4) + 17*powint(MZs,4) + 252*sqr(MWs)*sqr(MZs) - 280*MZs*trip(MWs) - 104*MWs*trip(MZs)) + MTs*(320*powint(MWs,4) + 7*powint(MZs,4) + 252*sqr(MWs)*sqr(MZs) - 560*MZs*trip(MWs) - 64*MWs*trip(MZs))))/128. + (MZs*powint(MWs - MZs,-5)*powint(Pi,-3)*trip(AL)*(-128*MWs*MZs*(-9*MZs*sqr(MWs) + 5*MWs*sqr(MZs) + 5*trip(MWs) - trip(MZs)) - 3*MTs*(640*powint(MWs,4) - 9*powint(MZs,4) + 150*sqr(MWs)*sqr(MZs) - 480*MZs*trip(MWs) + 4*MWs*trip(MZs))))/(256.*MWs)) + sqr(real(B0fin(MWs,MTs,0)))*((3*deltaAlpha*powint(MWs,-4)*powint(MWs - MZs,-4)*sqr(AL)*sqr(MZs)*(-4*MWs*MZs + 5*sqr(MWs) + sqr(MZs))*sqr(MTs*MWs + pow(MTs,2) - 2*pow(MWs,2)))/(64.*sqr(Pi)) - (powint(MWs,-4)*powint(MWs - MZs,-6)*powint(Pi,-3)*real(B0fin(MZs,0,0))*sqr(MTs - MWs)*sqr(MTs + 2*MWs)*trip(AL)*(-812*MZs*sqr(MWs) + 344*MWs*sqr(MZs) + 840*trip(MWs) - 57*trip(MZs))*trip(MZs))/1536. - (MZs*powint(MWs,-4)*powint(MWs - MZs,-6)*powint(Pi,-3)*real(B0fin(MZs,MTs,MTs))*sqr(MTs - MWs)*sqr(MTs + 2*MWs)*trip(AL)*(MZs*(160*powint(MWs,4) + 17*powint(MZs,4) + 252*sqr(MWs)*sqr(MZs) - 280*MZs*trip(MWs) - 104*MWs*trip(MZs)) + MTs*(320*powint(MWs,4) + 7*powint(MZs,4) + 252*sqr(MWs)*sqr(MZs) - 560*MZs*trip(MWs) - 64*MWs*trip(MZs))))/1536. - (MZs*powint(MWs,-5)*powint(MWs - MZs,-5)*powint(Pi,-3)*sqr(MTs*MWs + pow(MTs,2) - 2*pow(MWs,2))*trip(AL)*(128*MWs*MZs*(-9*MZs*sqr(MWs) + 5*MWs*sqr(MZs) + 5*trip(MWs) - trip(MZs)) + 3*MTs*(640*powint(MWs,4) - 9*powint(MZs,4) + 150*sqr(MWs)*sqr(MZs) - 480*MZs*trip(MWs) + 4*MWs*trip(MZs))))/9216.) + sqr(real(B0fin(MZs,0,0)))*((deltaAlpha*powint(MWs - MZs,-4)*sqr(AL)*sqr(MZs)*(-740*MWs*MZs + 1000*sqr(MWs) + 181*sqr(MZs)))/(96.*sqr(Pi)) - (MZs*powint(MWs - MZs,-6)*powint(Pi,-3)*real(B0fin(MZs,MTs,MTs))*trip(AL)*(MZs*(64000*powint(MWs,4) + 6229*powint(MZs,4) + 96384*sqr(MWs)*sqr(MZs) - 108160*MZs*trip(MWs) - 38608*MWs*trip(MZs)) + MTs*(128000*powint(MWs,4) + 3323*powint(MZs,4) + 91968*sqr(MWs)*sqr(MZs) - 216320*MZs*trip(MWs) - 26816*MWs*trip(MZs))))/13824. + (MZs*powint(MWs - MZs,-5)*powint(Pi,-3)*trip(AL)*(-256*MWs*MZs*(-1740*MZs*sqr(MWs) + 921*MWs*sqr(MZs) + 1000*trip(MWs) - 181*trip(MZs)) - 3*MTs*(256000*powint(MWs,4) - 3255*powint(MZs,4) + 57696*sqr(MWs)*sqr(MZs) - 176640*MZs*trip(MWs) + 704*MWs*trip(MZs))))/(82944.*MWs)) + sqr(A0fin(MTs))*((deltaAlpha*powint(MWs,-4)*powint(MWs - MZs,-4)*sqr(AL)*(1024*powint(MWs,4)*powint(MWs - MZs,4)*sqr(MZs) + 81*powint(MTs,4)*powint(MZs,4)*(-4*MWs*MZs + 5*sqr(MWs) + sqr(MZs)) + 64*MTs*MZs*(-16*MWs*MZs + 64*sqr(MWs) + 9*sqr(MZs))*trip(MWs)*trip(MWs - MZs) + 18*MWs*sqr(MZs)*trip(MTs)*(128*powint(MWs,4) + 9*powint(MZs,4) + 106*sqr(MWs)*sqr(MZs) - 200*MZs*trip(MWs) - 43*MWs*trip(MZs)) + sqr(MTs)*sqr(MWs)*sqr(MWs - MZs)*(4096*powint(MWs,4) - 495*powint(MZs,4) + 1984*sqr(MWs)*sqr(MZs) - 2048*MZs*trip(MWs) + 864*MWs*trip(MZs))))/(1728.*sqr(MTs)*sqr(MZs)*sqr(Pi)) + (powint(MWs,-4)*powint(MWs - MZs,-6)*powint(Pi,-3)*real(B0fin(MWs,0,0))*trip(AL)*(1024*(2*MWs - MZs)*powint(MWs,4)*powint(MWs - MZs,4)*sqr(MZs) + 64*MTs*MZs*trip(MWs)*trip(MWs - MZs)*(-72*MZs*sqr(MWs) + 34*MWs*sqr(MZs) + 128*trip(MWs) - 9*trip(MZs)) - 81*powint(MTs,4)*powint(MZs,4)*(14*MZs*sqr(MWs) - 6*MWs*sqr(MZs) - 14*trip(MWs) + trip(MZs)) + 18*MWs*sqr(MZs)*trip(MTs)*(-560*MZs*powint(MWs,4) + 320*powint(MWs,5) + 61*MWs*powint(MZs,4) - 9*powint(MZs,5) + 378*sqr(MZs)*trip(MWs) - 190*sqr(MWs)*trip(MZs)) + sqr(MTs)*sqr(MWs)*sqr(MWs - MZs)*(-5120*MZs*powint(MWs,4) + 8192*powint(MWs,5) - 1854*MWs*powint(MZs,4) + 495*powint(MZs,5) + 5248*sqr(MZs)*trip(MWs) + 800*sqr(MWs)*trip(MZs))))/(2304.*MZs*sqr(MTs)) - ((MTs - MWs)*(MTs + 2*MWs)*powint(MWs,-6)*powint(MWs - MZs,-6)*powint(Pi,-3)*real(B0fin(MWs,MTs,0))*trip(AL)*(1024*(2*MWs - MZs)*powint(MWs,4)*powint(MWs - MZs,4)*sqr(MZs) + 64*MTs*MZs*trip(MWs)*trip(MWs - MZs)*(-72*MZs*sqr(MWs) + 34*MWs*sqr(MZs) + 128*trip(MWs) - 9*trip(MZs)) - 81*powint(MTs,4)*powint(MZs,4)*(14*MZs*sqr(MWs) - 6*MWs*sqr(MZs) - 14*trip(MWs) + trip(MZs)) + 18*MWs*sqr(MZs)*trip(MTs)*(-560*MZs*powint(MWs,4) + 320*powint(MWs,5) + 61*MWs*powint(MZs,4) - 9*powint(MZs,5) + 378*sqr(MZs)*trip(MWs) - 190*sqr(MWs)*trip(MZs)) + sqr(MTs)*sqr(MWs)*sqr(MWs - MZs)*(-5120*MZs*powint(MWs,4) + 8192*powint(MWs,5) - 1854*MWs*powint(MZs,4) + 495*powint(MZs,5) + 5248*sqr(MZs)*trip(MWs) + 800*sqr(MWs)*trip(MZs))))/(13824.*MZs*sqr(MTs)) - (powint(MWs,-4)*powint(MWs - MZs,-6)*powint(Pi,-3)*real(B0fin(MZs,0,0))*trip(AL)*(1024*(40*MWs - 19*MZs)*powint(MWs,4)*powint(MWs - MZs,4)*sqr(MZs) + 64*MTs*MZs*trip(MWs)*trip(MWs - MZs)*(-1376*MZs*sqr(MWs) + 688*MWs*sqr(MZs) + 2560*trip(MWs) - 171*trip(MZs)) + 27*powint(MTs,4)*powint(MZs,4)*(-812*MZs*sqr(MWs) + 344*MWs*sqr(MZs) + 840*trip(MWs) - 57*trip(MZs)) + 18*MWs*sqr(MZs)*trip(MTs)*(-11008*MZs*powint(MWs,4) + 6400*powint(MWs,5) + 1187*MWs*powint(MZs,4) - 171*powint(MZs,5) + 7336*sqr(MZs)*trip(MWs) - 3744*sqr(MWs)*trip(MZs)) + sqr(MTs)*sqr(MWs)*sqr(MWs - MZs)*(-98304*MZs*powint(MWs,4) + 163840*powint(MWs,5) - 35352*MWs*powint(MZs,4) + 9405*powint(MZs,5) + 105984*sqr(MZs)*trip(MWs) + 17408*sqr(MWs)*trip(MZs))))/(41472.*MZs*sqr(MTs)) - (powint(MWs,-4)*powint(MWs - MZs,-6)*powint(MZs,-3)*powint(Pi,-3)*real(B0fin(MZs,MTs,MTs))*trip(AL)*(1024*powint(MWs,4)*powint(MWs - MZs,4)*(-40*MWs*MZs + 32*sqr(MWs) + 17*sqr(MZs))*trip(MZs) + 192*MTs*sqr(MZs)*trip(MWs)*trip(MWs - MZs)*(1024*powint(MWs,4) + 51*powint(MZs,4) + 1232*sqr(MWs)*sqr(MZs) - 1792*MZs*trip(MWs) - 272*MWs*trip(MZs)) + 81*powint(MTs,5)*powint(MZs,4)*(320*powint(MWs,4) + 7*powint(MZs,4) + 252*sqr(MWs)*sqr(MZs) - 560*MZs*trip(MWs) - 64*MWs*trip(MZs)) + 9*powint(MTs,4)*sqr(MZs)*(-43008*MZs*powint(MWs,6) + 16384*powint(MWs,7) - 810*MWs*powint(MZs,6) + 153*powint(MZs,7) - 434*powint(MZs,5)*sqr(MWs) + 41600*powint(MWs,5)*sqr(MZs) + 10392*powint(MZs,4)*trip(MWs) - 23872*powint(MWs,4)*trip(MZs)) + MWs*trip(MTs)*(-983040*MZs*powint(MWs,8) + 262144*powint(MWs,9) + 632256*powint(MWs,5)*powint(MZs,4) + 98544*powint(MWs,4)*powint(MZs,5) - 23643*MWs*powint(MZs,8) + 2754*powint(MZs,9) + 95778*powint(MZs,7)*sqr(MWs) + 1572864*powint(MWs,7)*sqr(MZs) - 248633*powint(MZs,6)*trip(MWs) - 1409024*powint(MWs,6)*trip(MZs)) + 3*MZs*sqr(MTs)*sqr(MWs)*sqr(MWs - MZs)*(-294912*MZs*powint(MWs,5) + 131072*powint(MWs,6) + 9288*MWs*powint(MZs,5) - 2805*powint(MZs,6) + 35232*powint(MZs,4)*sqr(MWs) + 276480*powint(MWs,4)*sqr(MZs) - 131072*trip(MWs)*trip(MZs))))/(124416.*sqr(MTs)) - (powint(MWs,-5)*powint(MWs - MZs,-5)*powint(MZs,-3)*powint(Pi,-3)*trip(AL)*(1024*MTs*powint(MWs,4)*powint(MWs - MZs,4)*sqr(MZs)*(-224*MWs*MZs + 896*sqr(MWs) + 45*sqr(MZs)) + 131072*powint(MWs,5)*powint(MWs - MZs,5)*trip(MZs) + 243*powint(MTs,5)*powint(MZs,4)*(640*powint(MWs,4) - 9*powint(MZs,4) + 150*sqr(MWs)*sqr(MZs) - 480*MZs*trip(MWs) + 4*MWs*trip(MZs)) + 64*MZs*sqr(MTs)*trip(MWs)*trip(MWs - MZs)*(32768*powint(MWs,4) - 1233*powint(MZs,4) + 10688*sqr(MWs)*sqr(MZs) - 16384*MZs*trip(MWs) + 1296*MWs*trip(MZs)) + 54*MWs*powint(MTs,4)*sqr(MZs)*(-26624*MZs*powint(MWs,5) + 16384*powint(MWs,6) + 1059*MWs*powint(MZs,5) - 273*powint(MZs,6) + 254*powint(MZs,4)*sqr(MWs) + 17280*powint(MWs,4)*sqr(MZs) - 8080*trip(MWs)*trip(MZs)) + 3*sqr(MWs)*sqr(MWs - MZs)*trip(MTs)*(-393216*MZs*powint(MWs,5) + 524288*powint(MWs,6) + 34176*MWs*powint(MZs,5) - 2457*powint(MZs,6) - 120960*powint(MZs,4)*sqr(MWs) + 430080*powint(MWs,4)*sqr(MZs) + 126976*trip(MWs)*trip(MZs))))/(746496.*sqr(MTs))) + (powint(MWs,-3)*powint(MWs - MZs,-4)*powint(MZs,-3)*powint(Pi,-3)*trip(AL)*(314928*powint(MZs,4)*sqr(MWs)*sqr(Pi)*(-8*MZs*sqr(MWs)*(6*MWs*MZs + 32*sqr(MWs) - 11*sqr(MZs)) + 3*sqr(MTs)*(-160*MZs*sqr(MWs) + 14*MWs*sqr(MZs) + 256*trip(MWs) + 9*trip(MZs)) - MTs*MWs*(-736*MZs*sqr(MWs) + 102*MWs*sqr(MZs) + 768*trip(MWs) + 61*trip(MZs))) + (-MTs + MWs)*(-MWs + MZs)*(-147456*MTs*sqr(MWs)*sqr(MWs - MZs)*(-32*MWs*MZs + 128*sqr(MWs) - 9*sqr(MZs))*sqr(MZs) - 2097152*trip(MWs)*trip(MWs - MZs)*trip(MZs) - 3456*MWs*MZs*sqr(MTs)*(-24576*MZs*powint(MWs,4) + 16384*powint(MWs,5) - 495*MWs*powint(MZs,4) - 81*powint(MZs,5) + 9216*sqr(MZs)*trip(MWs) - 448*sqr(MWs)*trip(MZs)) - 27*trip(MTs)*(-1572864*MZs*powint(MWs,5) + 2097152*powint(MWs,6) - 7776*MWs*powint(MZs,5) - 729*powint(MZs,6) - 38016*powint(MZs,4)*sqr(MWs) + 835584*powint(MWs,4)*sqr(MZs) - 32768*trip(MWs)*trip(MZs))) - (108*MWs*(-MTs + MWs)*sqr(MZs)*sqr(Pi)*(MTs*sqr(MZs)*(72089600*MZs*powint(MWs,5) - 29491200*powint(MWs,6) - 143701936*MWs*powint(MZs,5) + 29249013*powint(MZs,6) + 282442336*powint(MZs,4)*sqr(MWs) + 38978560*powint(MWs,4)*sqr(MZs) - 244339200*trip(MWs)*trip(MZs)) - 2*MZs*sqr(MTs)*(15564800*MZs*powint(MWs,5) - 9830400*powint(MWs,6) - 129022096*MWs*powint(MZs,5) + 27212085*powint(MZs,6) + 247517056*powint(MZs,4)*sqr(MWs) + 73507840*powint(MWs,4)*sqr(MZs) - 224294400*trip(MWs)*trip(MZs)) - 8*trip(MZs)*(1587200*MZs*powint(MWs,5) - 614400*powint(MWs,6) - 2360554*MWs*powint(MZs,5) + 477405*powint(MZs,6) + 4628392*powint(MZs,4)*sqr(MWs) + 236800*powint(MWs,4)*sqr(MZs) - 3847680*trip(MWs)*trip(MZs)) + 12*trip(MTs)*(-24166400*MZs*powint(MWs,5) + 8192000*powint(MWs,6) + 2992160*MWs*powint(MZs,5) - 222789*powint(MZs,6) - 8301376*powint(MZs,4)*sqr(MWs) + 17566720*powint(MWs,4)*sqr(MZs) + 2014720*trip(MWs)*trip(MZs))))/sqr(-4*MTs + MZs)))/(8.0621568e7*(MTs - MWs)) + sqr(real(B0fin(MZs,MTs,MTs)))*((deltaAlpha*powint(MWs - MZs,-4)*sqr(AL)*(sqr(MTs)*(2048*powint(MWs,4) + 569*powint(MZs,4) + 3936*sqr(MWs)*sqr(MZs) - 5120*MZs*trip(MWs) - 1352*MWs*trip(MZs)) + 2*MTs*MZs*(1024*powint(MWs,4) + 109*powint(MZs,4) + 2400*sqr(MWs)*sqr(MZs) - 2560*MZs*trip(MWs) - 1054*MWs*trip(MZs)) + sqr(MZs)*(512*powint(MWs,4) + 149*powint(MZs,4) + 1416*sqr(MWs)*sqr(MZs) - 1280*MZs*trip(MWs) - 716*MWs*trip(MZs))))/(864.*sqr(MZs)*sqr(Pi)) + (powint(MWs - MZs,-5)*powint(MZs,-3)*powint(Pi,-3)*trip(AL)*(-256*MWs*trip(MZs)*(-1792*MZs*powint(MWs,4) + 512*powint(MWs,5) + 865*MWs*powint(MZs,4) - 149*powint(MZs,5) + 2696*sqr(MZs)*trip(MWs) - 2132*sqr(MWs)*trip(MZs)) - 3*trip(MTs)*(-1441792*MZs*powint(MWs,5) + 524288*powint(MWs,6) - 29504*MWs*powint(MZs,5) - 6975*powint(MZs,6) + 425728*powint(MZs,4)*sqr(MWs) + 1445888*powint(MWs,4)*sqr(MZs) - 892928*trip(MWs)*trip(MZs)) + 2*MZs*sqr(MTs)*(3080192*MZs*powint(MWs,5) - 1048576*powint(MWs,6) + 137824*MWs*powint(MZs,5) + 4995*powint(MZs,6) - 705632*powint(MZs,4)*sqr(MWs) - 3659776*powint(MWs,4)*sqr(MZs) + 2265088*trip(MWs)*trip(MZs)) + MTs*sqr(MZs)*(2916352*MZs*powint(MWs,5) - 917504*powint(MWs,6) + 63680*MWs*powint(MZs,5) + 7965*powint(MZs,6) - 875872*powint(MZs,4)*sqr(MWs) - 3955712*powint(MWs,4)*sqr(MZs) + 2686976*trip(MWs)*trip(MZs))))/(746496.*MWs)) + real(B0fin(MZs,0,0))*(-(AL*(40*MWs - 19*MZs)*MZs*sqr(deltaAlpha))/(8.*Pi*sqr(MWs - MZs)) + (deltaAlpha*powint(MWs - MZs,-3)*sqr(AL)*(128*MWs*MZs*(-59*MWs*MZs + 40*sqr(MWs) + 19*sqr(MZs)) + 3*MTs*(-2752*MZs*sqr(MWs) + 296*MWs*sqr(MZs) + 5120*trip(MWs) + 171*trip(MZs))))/(1728.*MWs*sqr(Pi)) + (powint(MWs - MZs,-6)*powint(Pi,-3)*sqr(real(B0fin(MZs,MTs,MTs)))*trip(AL)*(sqr(MZs)*(106496*MZs*powint(MWs,4) - 40960*powint(MWs,5) - 36944*MWs*powint(MZs,4) + 5719*powint(MZs,5) - 139648*sqr(MZs)*trip(MWs) + 96832*sqr(MWs)*trip(MZs)) + 2*MTs*MZs*(212992*MZs*powint(MWs,4) - 81920*powint(MWs,5) - 34144*MWs*powint(MZs,4) + 4121*powint(MZs,5) - 227456*sqr(MZs)*trip(MWs) + 134912*sqr(MWs)*trip(MZs)) + sqr(MTs)*(425984*MZs*powint(MWs,4) - 163840*powint(MWs,5) - 86000*MWs*powint(MZs,4) + 14263*powint(MZs,5) - 351232*sqr(MZs)*trip(MWs) + 152320*sqr(MWs)*trip(MZs))))/(41472.*MZs) - (powint(MWs - MZs,-5)*powint(Pi,-3)*trip(AL)*(768*MTs*MWs*MZs*sqr(MWs - MZs)*(-2752*MZs*sqr(MWs) + 296*MWs*sqr(MZs) + 5120*trip(MWs) + 171*trip(MZs)) + 4*MWs*sqr(MZs)*(163840*powint(MWs,5) - 1874085*powint(MZs,5)*sqr(Pi) - 2048*MZs*powint(MWs,4)*(278 + 975*sqr(Pi)) + 4*MWs*powint(MZs,4)*(19456 + 2076261*sqr(Pi)) + 96*sqr(MZs)*(7552 + 92785*sqr(Pi))*trip(MWs) - 32*sqr(MWs)*(12416 + 406557*sqr(Pi))*trip(MZs)) + 9*sqr(MTs)*(-1048576*MZs*powint(MWs,5) + 655360*powint(MWs,6) + 6741*MWs*powint(MZs,5) + 1539*powint(MZs,6) - 13400*powint(MZs,4)*sqr(MWs) + 540672*powint(MWs,4)*sqr(MZs) - 142336*trip(MWs)*trip(MZs))))/(1.492992e6*MZs*sqr(MWs)) + real(B0fin(MZs,MTs,MTs))*((deltaAlpha*powint(MWs - MZs,-4)*sqr(AL)*(MZs*(-984*MZs*sqr(MWs) + 696*MWs*sqr(MZs) + 640*trip(MWs) - 163*trip(MZs)) + MTs*(-1968*MZs*sqr(MWs) + 582*MWs*sqr(MZs) + 1280*trip(MWs) - 83*trip(MZs))))/(144.*sqr(Pi)) - (powint(MWs - MZs,-5)*powint(Pi,-3)*trip(AL)*(256*MWs*sqr(MZs)*(640*powint(MWs,4) + 163*powint(MZs,4) + 1680*sqr(MWs)*sqr(MZs) - 1624*MZs*trip(MWs) - 859*MWs*trip(MZs)) + MTs*MZs*(-1617920*MZs*powint(MWs,4) + 819200*powint(MWs,5) + 22592*MWs*powint(MZs,4) + 8775*powint(MZs,5) + 1300224*sqr(MZs)*trip(MWs) - 359936*sqr(MWs)*trip(MZs)) + 3*sqr(MTs)*(-524288*MZs*powint(MWs,4) + 327680*powint(MWs,5) + 12128*MWs*powint(MZs,4) + 1395*powint(MZs,5) + 224256*sqr(MZs)*trip(MWs) - 98816*sqr(MWs)*trip(MZs))))/(124416.*MWs*MZs))) + real(B0fin(MWs,MTs,0))*((-3*AL*(MTs - MWs)*(MTs + 2*MWs)*(2*MWs - MZs)*MZs*sqr(deltaAlpha))/(8.*Pi*sqr(MWs)*sqr(MWs - MZs)) - (35*(MTs - MWs)*(MTs + 2*MWs)*powint(MWs - MZs,-6)*powint(Pi,-3)*sqr(real(B0fin(MZs,0,0)))*trip(AL)*(-448*MZs*sqr(MWs) + 188*MWs*sqr(MZs) + 480*trip(MWs) - 31*trip(MZs))*trip(MZs))/(4608.*sqr(MWs)) + (deltaAlpha*(MTs - MWs)*(MTs + 2*MWs)*powint(MWs,-3)*powint(MWs - MZs,-3)*sqr(AL)*(128*MWs*MZs*(-3*MWs*MZs + 2*sqr(MWs) + sqr(MZs)) + 3*MTs*(-144*MZs*sqr(MWs) + 14*MWs*sqr(MZs) + 256*trip(MWs) + 9*trip(MZs))))/(576.*sqr(Pi)) + (powint(MWs - MZs,-6)*powint(Pi,-3)*sqr(real(B0fin(MZs,MTs,MTs)))*trip(AL)*(sqr(MTs)*(-51200*MZs*powint(MWs,6) + 16384*powint(MWs,7) - 1514*MWs*powint(MZs,6) + 295*powint(MZs,7) - 54*powint(MZs,5)*sqr(MWs) + 55296*powint(MWs,5)*sqr(MZs) + 14968*powint(MZs,4)*trip(MWs) - 32960*powint(MWs,4)*trip(MZs)) + 2*sqr(MWs)*sqr(MZs)*(-5376*MZs*powint(MWs,4) + 2048*powint(MWs,5) + 1884*MWs*powint(MZs,4) - 295*powint(MZs,5) + 7072*sqr(MZs)*trip(MWs) - 4928*sqr(MWs)*trip(MZs)) + powint(MTs,4)*(21504*MZs*powint(MWs,4) - 8192*powint(MWs,5) - 4188*MWs*powint(MZs,4) + 775*powint(MZs,5) - 17920*sqr(MZs)*trip(MWs) + 7616*sqr(MWs)*trip(MZs)) + MTs*MWs*MZs*(-45056*MZs*powint(MWs,5) + 16384*powint(MWs,6) - 2624*MWs*powint(MZs,5) + 295*powint(MZs,6) + 11792*powint(MZs,4)*sqr(MWs) + 51584*powint(MWs,4)*sqr(MZs) - 34400*trip(MWs)*trip(MZs)) + trip(MTs)*(13312*MZs*powint(MWs,5) - 8192*powint(MWs,6) - 2657*MWs*powint(MZs,5) + 370*powint(MZs,6) + 9476*powint(MZs,4)*sqr(MWs) + 3584*powint(MWs,4)*sqr(MZs) - 15488*trip(MWs)*trip(MZs))))/(13824.*MZs*sqr(MWs)) + real(B0fin(MZs,MTs,MTs))*((deltaAlpha*powint(MWs - MZs,-4)*sqr(AL)*(MTs*MWs + sqr(MTs) - 2*sqr(MWs))*(MZs*(-100*MZs*sqr(MWs) + 71*MWs*sqr(MZs) + 64*trip(MWs) - 17*trip(MZs)) + MTs*(-200*MZs*sqr(MWs) + 61*MWs*sqr(MZs) + 128*trip(MWs) - 7*trip(MZs))))/(96.*sqr(MWs)*sqr(Pi)) - (powint(MWs,-3)*powint(MWs - MZs,-5)*powint(Pi,-3)*(MTs*MWs + sqr(MTs) - 2*sqr(MWs))*trip(AL)*(128*MWs*sqr(MZs)*(64*powint(MWs,4) + 17*powint(MZs,4) + 171*sqr(MWs)*sqr(MZs) - 164*MZs*trip(MWs) - 88*MWs*trip(MZs)) + MTs*MZs*(-81920*MZs*powint(MWs,4) + 40960*powint(MWs,5) + 872*MWs*powint(MZs,4) + 459*powint(MZs,5) + 66240*sqr(MZs)*trip(MWs) - 18376*sqr(MWs)*trip(MZs)) + 3*sqr(MTs)*(-26624*MZs*powint(MWs,4) + 16384*powint(MWs,5) + 632*MWs*powint(MZs,4) + 63*powint(MZs,5) + 11520*sqr(MZs)*trip(MWs) - 4720*sqr(MWs)*trip(MZs))))/(41472.*MZs)) + real(B0fin(MZs,0,0))*((deltaAlpha*powint(MWs - MZs,-4)*sqr(AL)*(MTs*MWs + sqr(MTs) - 2*sqr(MWs))*sqr(MZs)*(-77*MWs*MZs + 100*sqr(MWs) + 19*sqr(MZs)))/(32.*sqr(MWs)*sqr(Pi)) - (MZs*powint(MWs,-3)*powint(MWs - MZs,-5)*powint(Pi,-3)*(MTs*MWs + sqr(MTs) - 2*sqr(MWs))*trip(AL)*(128*MWs*MZs*(-177*MZs*sqr(MWs) + 96*MWs*sqr(MZs) + 100*trip(MWs) - 19*trip(MZs)) + 3*MTs*(12800*powint(MWs,4) - 171*powint(MZs,4) + 2936*sqr(MWs)*sqr(MZs) - 9216*MZs*trip(MWs) + 56*MWs*trip(MZs))))/13824. + (MZs*powint(MWs - MZs,-6)*powint(Pi,-3)*real(B0fin(MZs,MTs,MTs))*trip(AL)*(2*MZs*sqr(MWs)*(3200*powint(MWs,4) + 325*powint(MZs,4) + 4928*sqr(MWs)*sqr(MZs) - 5504*MZs*trip(MWs) - 2004*MWs*trip(MZs)) + trip(MTs)*(-6400*powint(MWs,4) - 155*powint(MZs,4) - 4816*sqr(MWs)*sqr(MZs) + 11008*MZs*trip(MWs) + 1308*MWs*trip(MZs)) + MTs*MWs*(-25216*MZs*powint(MWs,4) + 12800*powint(MWs,5) + 2314*MWs*powint(MZs,4) - 325*powint(MZs,5) + 15136*sqr(MZs)*trip(MWs) - 7544*sqr(MWs)*trip(MZs)) + sqr(MTs)*(7808*MZs*powint(MWs,4) - 6400*powint(MWs,5) + 1849*MWs*powint(MZs,4) - 325*powint(MZs,5) + 688*sqr(MZs)*trip(MWs) - 3620*sqr(MWs)*trip(MZs))))/(2304.*sqr(MWs))) + (powint(MWs,-4)*powint(MWs - MZs,-5)*powint(Pi,-3)*trip(AL)*(34992*powint(MZs,4)*sqr(MWs)*sqr(Pi)*(-3*MTs*sqr(MWs)*(-9*MWs*MZs + 12*sqr(MWs) + 2*sqr(MZs)) + (-12*MWs*MZs + 14*sqr(MWs) + 3*sqr(MZs))*trip(MTs) + 2*(-12*MWs*MZs + 14*sqr(MWs) + 3*sqr(MZs))*trip(MWs)) + 12*MWs*(MTs + 2*MWs)*sqr(MTs - MWs)*sqr(Pi)*trip(MZs)*(32000*powint(MWs,4) + 31827*powint(MZs,4) + 249360*sqr(MWs)*sqr(MZs) - 185760*MZs*trip(MWs) - 147272*MWs*trip(MZs)) - (MTs + 2*MWs)*(MWs - MZs)*sqr(MTs - MWs)*(16384*(2*MWs - MZs)*sqr(MWs)*sqr(MWs - MZs)*sqr(MZs) + 768*MTs*MWs*MZs*(256*powint(MWs,4) - 9*powint(MZs,4) + 158*sqr(MWs)*sqr(MZs) - 400*MZs*trip(MWs) - 5*MWs*trip(MZs)) + 9*sqr(MTs)*(-20480*MZs*powint(MWs,4) + 32768*powint(MWs,5) - 414*MWs*powint(MZs,4) - 81*powint(MZs,5) + 7168*sqr(MZs)*trip(MWs) + 320*sqr(MWs)*trip(MZs)))))/(497664.*(MTs - MWs)*MZs)) + real(B0fin(MWs,0,0))*((9*AL*(2*MWs - MZs)*MZs*sqr(deltaAlpha))/(4.*Pi*sqr(MWs - MZs)) + (35*powint(MWs - MZs,-6)*powint(Pi,-3)*sqr(real(B0fin(MZs,0,0)))*trip(AL)*(-448*MZs*sqr(MWs) + 188*MWs*sqr(MZs) + 480*trip(MWs) - 31*trip(MZs))*trip(MZs))/768. + (9*powint(MWs,-4)*powint(MWs - MZs,-6)*powint(Pi,-3)*sqr(MTs - MWs)*sqr(MTs + 2*MWs)*sqr(real(B0fin(MWs,MTs,0)))*trip(AL)*(-14*MZs*sqr(MWs) + 6*MWs*sqr(MZs) + 14*trip(MWs) - trip(MZs))*trip(MZs))/256. - (deltaAlpha*powint(MWs - MZs,-3)*sqr(AL)*(128*MWs*MZs*(-3*MWs*MZs + 2*sqr(MWs) + sqr(MZs)) + 3*MTs*(-144*MZs*sqr(MWs) + 14*MWs*sqr(MZs) + 256*trip(MWs) + 9*trip(MZs))))/(96.*MWs*sqr(Pi)) + (powint(MWs - MZs,-6)*powint(Pi,-3)*sqr(real(B0fin(MZs,MTs,MTs)))*trip(AL)*(sqr(MTs)*(-21504*MZs*powint(MWs,4) + 8192*powint(MWs,5) + 4188*MWs*powint(MZs,4) - 775*powint(MZs,5) + 17920*sqr(MZs)*trip(MWs) - 7616*sqr(MWs)*trip(MZs)) + 2*MTs*MZs*(-10752*MZs*powint(MWs,4) + 4096*powint(MWs,5) + 1716*MWs*powint(MZs,4) - 185*powint(MZs,5) + 11552*sqr(MZs)*trip(MWs) - 6832*sqr(MWs)*trip(MZs)) + sqr(MZs)*(-5376*MZs*powint(MWs,4) + 2048*powint(MWs,5) + 1884*MWs*powint(MZs,4) - 295*powint(MZs,5) + 7072*sqr(MZs)*trip(MWs) - 4928*sqr(MWs)*trip(MZs))))/(2304.*MZs) + (powint(MWs - MZs,-5)*powint(Pi,-3)*trip(AL)*(768*MTs*MWs*MZs*sqr(MWs - MZs)*(-144*MZs*sqr(MWs) + 14*MWs*sqr(MZs) + 256*trip(MWs) + 9*trip(MZs)) + 4*MWs*sqr(MZs)*(8192*powint(MWs,5) - 95481*powint(MZs,5)*sqr(Pi) - 256*MZs*powint(MWs,4)*(112 + 375*sqr(Pi)) + 4*MWs*powint(MZs,4)*(1024 + 103893*sqr(Pi)) + 72*sqr(MZs)*(512 + 6039*sqr(Pi))*trip(MWs) - 32*sqr(MWs)*(640 + 20097*sqr(Pi))*trip(MZs)) + 9*sqr(MTs)*(-53248*MZs*powint(MWs,5) + 32768*powint(MWs,6) + 333*MWs*powint(MZs,5) + 81*powint(MZs,6) - 734*powint(MZs,4)*sqr(MWs) + 27648*powint(MWs,4)*sqr(MZs) - 6848*trip(MWs)*trip(MZs))))/(82944.*MZs*sqr(MWs)) + real(B0fin(MZs,0,0))*((-3*deltaAlpha*powint(MWs - MZs,-4)*sqr(AL)*sqr(MZs)*(-77*MWs*MZs + 100*sqr(MWs) + 19*sqr(MZs)))/(16.*sqr(Pi)) + (MZs*powint(MWs - MZs,-6)*powint(Pi,-3)*real(B0fin(MZs,MTs,MTs))*trip(AL)*(MZs*(3200*powint(MWs,4) + 325*powint(MZs,4) + 4928*sqr(MWs)*sqr(MZs) - 5504*MZs*trip(MWs) - 2004*MWs*trip(MZs)) + MTs*(6400*powint(MWs,4) + 155*powint(MZs,4) + 4816*sqr(MWs)*sqr(MZs) - 11008*MZs*trip(MWs) - 1308*MWs*trip(MZs))))/384. + (MZs*powint(MWs - MZs,-5)*powint(Pi,-3)*trip(AL)*(128*MWs*MZs*(-177*MZs*sqr(MWs) + 96*MWs*sqr(MZs) + 100*trip(MWs) - 19*trip(MZs)) + 3*MTs*(12800*powint(MWs,4) - 171*powint(MZs,4) + 2936*sqr(MWs)*sqr(MZs) - 9216*MZs*trip(MWs) + 56*MWs*trip(MZs))))/(2304.*MWs)) + real(B0fin(MZs,MTs,MTs))*((deltaAlpha*powint(MWs - MZs,-4)*sqr(AL)*(MTs*(200*MZs*sqr(MWs) - 61*MWs*sqr(MZs) - 128*trip(MWs) + 7*trip(MZs)) + MZs*(100*MZs*sqr(MWs) - 71*MWs*sqr(MZs) - 64*trip(MWs) + 17*trip(MZs))))/(16.*sqr(Pi)) + (powint(MWs - MZs,-5)*powint(Pi,-3)*trip(AL)*(128*MWs*sqr(MZs)*(64*powint(MWs,4) + 17*powint(MZs,4) + 171*sqr(MWs)*sqr(MZs) - 164*MZs*trip(MWs) - 88*MWs*trip(MZs)) + MTs*MZs*(-81920*MZs*powint(MWs,4) + 40960*powint(MWs,5) + 872*MWs*powint(MZs,4) + 459*powint(MZs,5) + 66240*sqr(MZs)*trip(MWs) - 18376*sqr(MWs)*trip(MZs)) + 3*sqr(MTs)*(-26624*MZs*powint(MWs,4) + 16384*powint(MWs,5) + 632*MWs*powint(MZs,4) + 63*powint(MZs,5) + 11520*sqr(MZs)*trip(MWs) - 4720*sqr(MWs)*trip(MZs))))/(6912.*MWs*MZs)) + real(B0fin(MWs,MTs,0))*((-9*deltaAlpha*powint(MWs - MZs,-4)*sqr(AL)*(MTs*MWs + sqr(MTs) - 2*sqr(MWs))*sqr(MZs)*(-4*MWs*MZs + 5*sqr(MWs) + sqr(MZs)))/(16.*sqr(MWs)*sqr(Pi)) + ((MTs - MWs)*(MTs + 2*MWs)*powint(MWs - MZs,-6)*powint(Pi,-3)*real(B0fin(MZs,0,0))*trip(AL)*(-812*MZs*sqr(MWs) + 344*MWs*sqr(MZs) + 840*trip(MWs) - 57*trip(MZs))*trip(MZs))/(128.*sqr(MWs)) + (MZs*powint(MWs,-3)*powint(MWs - MZs,-5)*powint(Pi,-3)*(MTs*MWs + sqr(MTs) - 2*sqr(MWs))*trip(AL)*(128*MWs*MZs*(-9*MZs*sqr(MWs) + 5*MWs*sqr(MZs) + 5*trip(MWs) - trip(MZs)) + 3*MTs*(640*powint(MWs,4) - 9*powint(MZs,4) + 150*sqr(MWs)*sqr(MZs) - 480*MZs*trip(MWs) + 4*MWs*trip(MZs))))/768. + (MZs*powint(MWs - MZs,-6)*powint(Pi,-3)*real(B0fin(MZs,MTs,MTs))*trip(AL)*(-2*MZs*sqr(MWs)*(160*powint(MWs,4) + 17*powint(MZs,4) + 252*sqr(MWs)*sqr(MZs) - 280*MZs*trip(MWs) - 104*MWs*trip(MZs)) + trip(MTs)*(320*powint(MWs,4) + 7*powint(MZs,4) + 252*sqr(MWs)*sqr(MZs) - 560*MZs*trip(MWs) - 64*MWs*trip(MZs)) + sqr(MTs)*(-400*MZs*powint(MWs,4) + 320*powint(MWs,5) - 97*MWs*powint(MZs,4) + 17*powint(MZs,5) - 28*sqr(MZs)*trip(MWs) + 188*sqr(MWs)*trip(MZs)) + MTs*MWs*(1280*MZs*powint(MWs,4) - 640*powint(MWs,5) - 118*MWs*powint(MZs,4) + 17*powint(MZs,5) - 784*sqr(MZs)*trip(MWs) + 380*sqr(MWs)*trip(MZs))))/(128.*sqr(MWs)))) + real(B0fin(MZs,MTs,MTs))*(-(AL*sqr(deltaAlpha)*(MTs*(-80*MWs*MZs + 64*sqr(MWs) + 7*sqr(MZs)) + MZs*(-40*MWs*MZs + 32*sqr(MWs) + 17*sqr(MZs))))/(24.*MZs*Pi*sqr(MWs - MZs)) + (deltaAlpha*powint(MWs - MZs,-3)*sqr(AL)*(128*MWs*sqr(MZs)*(-72*MZs*sqr(MWs) + 57*MWs*sqr(MZs) + 32*trip(MWs) - 17*trip(MZs)) + MTs*MZs*(20480*powint(MWs,4) - 459*powint(MZs,4) + 22368*sqr(MWs)*sqr(MZs) - 36864*MZs*trip(MWs) - 1880*MWs*trip(MZs)) + 3*sqr(MTs)*(8192*powint(MWs,4) - 63*powint(MZs,4) + 4032*sqr(MWs)*sqr(MZs) - 12288*MZs*trip(MWs) - 1088*MWs*trip(MZs))))/(5184.*MWs*sqr(MZs)*sqr(Pi)) + (powint(MWs - MZs,-5)*powint(MZs,-3)*powint(Pi,-3)*trip(AL)*(104976*powint(MZs,4)*sqr(Pi)*(MZs*(-104*MZs*sqr(MWs) + 72*MWs*sqr(MZs) + 64*trip(MWs) - 17*trip(MZs)) + MTs*(-208*MZs*sqr(MWs) + 72*MWs*sqr(MZs) + 128*trip(MWs) - 7*trip(MZs))) - (36*sqr(MZs)*sqr(Pi)*(2*MZs*sqr(MTs)*(-10854400*MZs*powint(MWs,6) + 2457600*powint(MWs,7) + 1810431*MWs*powint(MZs,6) - 201571*powint(MZs,7) - 6033672*powint(MZs,5)*sqr(MWs) + 19768320*powint(MWs,5)*sqr(MZs) + 13251872*powint(MZs,4)*trip(MWs) - 19484160*powint(MWs,4)*trip(MZs)) - 2*trip(MTs)*(-32358400*MZs*powint(MWs,6) + 8192000*powint(MWs,7) - 2690669*MWs*powint(MZs,6) + 222789*powint(MZs,7) + 9784416*powint(MZs,5)*sqr(MWs) + 41733120*powint(MWs,5)*sqr(MZs) - 8547776*powint(MZs,4)*trip(MWs) - 15859200*powint(MWs,4)*trip(MZs)) + MWs*powint(MZs,4)*(1297920*MZs*powint(MWs,4) - 204800*powint(MWs,5) - 1077408*MWs*powint(MZs,4) + 200335*powint(MZs,5) - 2512000*sqr(MZs)*trip(MWs) + 2355488*sqr(MWs)*trip(MZs)) + 3*MTs*trip(MZs)*(-2595840*MZs*powint(MWs,5) + 409600*powint(MWs,6) - 591632*MWs*powint(MZs,5) + 31827*powint(MZs,6) + 2529936*powint(MZs,4)*sqr(MWs) + 5081600*powint(MWs,4)*sqr(MZs) - 5044096*trip(MWs)*trip(MZs))))/(sqr(MWs)*sqr(-4*MTs + MZs)) + ((-MWs + MZs)*(16384*sqr(MWs)*sqr(MWs - MZs)*(-40*MWs*MZs + 32*sqr(MWs) + 17*sqr(MZs))*trip(MZs) + 256*MTs*MWs*sqr(MZs)*(-44032*MZs*powint(MWs,4) + 16384*powint(MWs,5) + 973*MWs*powint(MZs,4) + 459*powint(MZs,5) + 44448*sqr(MZs)*trip(MWs) - 18232*sqr(MWs)*trip(MZs)) + 3*MZs*sqr(MTs)*(-7995392*MZs*powint(MWs,5) + 3670016*powint(MWs,6) + 40968*MWs*powint(MZs,5) + 4131*powint(MZs,6) + 280928*powint(MZs,4)*sqr(MWs) + 6537216*powint(MWs,4)*sqr(MZs) - 2017280*trip(MWs)*trip(MZs)) + 9*trip(MTs)*(-1835008*MZs*powint(MWs,5) + 1048576*powint(MWs,6) + 16560*MWs*powint(MZs,5) + 567*powint(MZs,6) + 46912*powint(MZs,4)*sqr(MWs) + 1130496*powint(MWs,4)*sqr(MZs) - 581632*trip(MWs)*trip(MZs))))/sqr(MWs)))/4.478976e6) + A0fin(MTs)*((AL*sqr(deltaAlpha)*(32*MZs*sqr(MWs)*sqr(MWs - MZs) + 9*(2*MWs - MZs)*sqr(MTs)*sqr(MZs) + MTs*MWs*(-80*MZs*sqr(MWs) + 25*MWs*sqr(MZs) + 64*trip(MWs) - 9*trip(MZs))))/(24.*MTs*MZs*Pi*sqr(MWs)*sqr(MWs - MZs)) + (MZs*powint(MWs - MZs,-6)*powint(Pi,-3)*sqr(real(B0fin(MZs,0,0)))*trip(AL)*(64*MZs*sqr(MWs)*sqr(MWs - MZs)*(-740*MWs*MZs + 1000*sqr(MWs) + 181*sqr(MZs)) + 105*sqr(MTs)*sqr(MZs)*(-448*MZs*sqr(MWs) + 188*MWs*sqr(MZs) + 480*trip(MWs) - 31*trip(MZs)) + MTs*MWs*(-216320*MZs*powint(MWs,4) + 128000*powint(MWs,5) + 23063*MWs*powint(MZs,4) - 3255*powint(MZs,5) + 142368*sqr(MZs)*trip(MWs) - 73856*sqr(MWs)*trip(MZs))))/(13824.*MTs*sqr(MWs)) + (powint(MWs - MZs,-6)*powint(MZs,-3)*powint(Pi,-3)*sqr(real(B0fin(MZs,MTs,MTs)))*trip(AL)*(64*sqr(MWs)*sqr(MWs - MZs)*trip(MZs)*(512*powint(MWs,4) + 149*powint(MZs,4) + 1416*sqr(MWs)*sqr(MZs) - 1280*MZs*trip(MWs) - 716*MWs*trip(MZs)) + 3*MTs*MWs*sqr(MZs)*(-278528*MZs*powint(MWs,6) + 65536*powint(MWs,7) + 12309*MWs*powint(MZs,6) - 885*powint(MZs,7) - 89152*powint(MZs,5)*sqr(MWs) + 507904*powint(MWs,5)*sqr(MZs) + 287648*powint(MZs,4)*trip(MWs) - 504832*powint(MWs,4)*trip(MZs)) + trip(MTs)*(-983040*MZs*powint(MWs,7) + 262144*powint(MWs,8) + 707712*powint(MWs,4)*powint(MZs,4) + 23913*MWs*powint(MZs,7) - 3330*powint(MZs,8) - 71873*powint(MZs,6)*sqr(MWs) + 1499136*powint(MWs,6)*sqr(MZs) - 152736*powint(MZs,5)*trip(MWs) - 1289216*powint(MWs,5)*trip(MZs)) + 3*MZs*sqr(MTs)*(-524288*MZs*powint(MWs,7) + 131072*powint(MWs,8) + 453632*powint(MWs,4)*powint(MZs,4) + 4542*MWs*powint(MZs,7) - 885*powint(MZs,8) + 23254*powint(MZs,6)*sqr(MWs) + 882688*powint(MWs,6)*sqr(MZs) - 150112*powint(MZs,5)*trip(MWs) - 818688*powint(MWs,5)*trip(MZs)) + 9*powint(MTs,4)*sqr(MZs)*(-21504*MZs*powint(MWs,4) + 8192*powint(MWs,5) + 4188*MWs*powint(MZs,4) - 775*powint(MZs,5) + 17920*sqr(MZs)*trip(MWs) - 7616*sqr(MWs)*trip(MZs))))/(124416.*MTs*sqr(MWs)) - (deltaAlpha*powint(MWs,-3)*powint(MWs - MZs,-3)*sqr(AL)*(32*MTs*MZs*sqr(MWs)*sqr(MWs - MZs)*(-160*MWs*MZs + 640*sqr(MWs) + 9*sqr(MZs)) + 4096*sqr(MZs)*trip(MWs)*trip(MWs - MZs) + 27*sqr(MZs)*trip(MTs)*(-144*MZs*sqr(MWs) + 14*MWs*sqr(MZs) + 256*trip(MWs) + 9*trip(MZs)) + 3*MWs*sqr(MTs)*(-12288*MZs*powint(MWs,4) + 8192*powint(MWs,5) - 1089*MWs*powint(MZs,4) + 465*powint(MZs,5) + 6336*sqr(MZs)*trip(MWs) - 1616*sqr(MWs)*trip(MZs))))/(5184.*MTs*sqr(MZs)*sqr(Pi)) + (3*MZs*powint(MWs - MZs,-6)*powint(Pi,-3)*sqr(real(B0fin(MWs,0,0)))*trip(AL)*(32*MZs*sqr(MWs)*sqr(MWs - MZs)*(-4*MWs*MZs + 5*sqr(MWs) + sqr(MZs)) - 9*sqr(MTs)*sqr(MZs)*(14*MZs*sqr(MWs) - 6*MWs*sqr(MZs) - 14*trip(MWs) + trip(MZs)) + MTs*MWs*(-560*MZs*powint(MWs,4) + 320*powint(MWs,5) + 61*MWs*powint(MZs,4) - 9*powint(MZs,5) + 378*sqr(MZs)*trip(MWs) - 190*sqr(MWs)*trip(MZs))))/(128.*MTs*sqr(MWs)) + (MZs*powint(MWs,-6)*powint(MWs - MZs,-6)*powint(Pi,-3)*sqr(MTs*MWs + pow(MTs,2) - 2*pow(MWs,2))*sqr(real(B0fin(MWs,MTs,0)))*trip(AL)*(32*MZs*sqr(MWs)*sqr(MWs - MZs)*(-4*MWs*MZs + 5*sqr(MWs) + sqr(MZs)) - 9*sqr(MTs)*sqr(MZs)*(14*MZs*sqr(MWs) - 6*MWs*sqr(MZs) - 14*trip(MWs) + trip(MZs)) + MTs*MWs*(-560*MZs*powint(MWs,4) + 320*powint(MWs,5) + 61*MWs*powint(MZs,4) - 9*powint(MZs,5) + 378*sqr(MZs)*trip(MWs) - 190*sqr(MWs)*trip(MZs))))/(1536.*MTs) + (powint(MWs,-4)*powint(MWs - MZs,-5)*powint(MZs,-3)*powint(Pi,-3)*trip(AL)*((104976*powint(MZs,4)*sqr(MWs)*sqr(Pi)*(-3*sqr(MZs)*(-12*MWs*MZs + 14*sqr(MWs) + 3*sqr(MZs))*trip(MTs) + 32*(2*MWs - MZs)*MZs*sqr(MWs - MZs)*trip(MWs) + sqr(MTs)*sqr(MWs)*(208*MZs*sqr(MWs) - 72*MWs*sqr(MZs) - 128*trip(MWs) + 7*trip(MZs)) + 16*MTs*sqr(MWs)*(8*powint(MWs,4) + 2*powint(MZs,4) + 16*sqr(MWs)*sqr(MZs) - 17*MZs*trip(MWs) - 9*MWs*trip(MZs))))/(MTs - MWs) + (36*MWs*sqr(MZs)*sqr(Pi)*(-32*powint(MZs,4)*sqr(MWs)*sqr(MWs - MZs)*(-30960*MZs*sqr(MWs) + 31200*MWs*sqr(MZs) + 6400*trip(MWs) - 10609*trip(MZs)) - 48*powint(MTs,4)*trip(MZs)*(32000*powint(MWs,4) + 31827*powint(MZs,4) + 249360*sqr(MWs)*sqr(MZs) - 185760*MZs*trip(MWs) - 147272*MWs*trip(MZs)) + MZs*sqr(MTs)*(-21708800*MZs*powint(MWs,7) + 4915200*powint(MWs,8) + 37493824*powint(MWs,4)*powint(MZs,4) + 802522*MWs*powint(MZs,7) - 95481*powint(MZs,8) + 1564702*powint(MZs,6)*sqr(MWs) + 41175040*powint(MWs,6)*sqr(MZs) - 15124912*powint(MZs,5)*trip(MWs) - 48962560*powint(MWs,5)*trip(MZs)) + trip(MTs)*(64716800*MZs*powint(MWs,7) - 16384000*powint(MWs,8) + 26780032*powint(MWs,4)*powint(MZs,4) - 5507802*MWs*powint(MZs,7) + 763848*powint(MZs,8) + 18435034*powint(MZs,6)*sqr(MWs) - 83466240*powint(MWs,6)*sqr(MZs) - 35996352*powint(MZs,5)*trip(MWs) + 30182400*powint(MWs,5)*trip(MZs)) + 32*MTs*sqr(MWs)*trip(MZs)*(-268960*MZs*powint(MWs,4) + 38400*powint(MWs,5) + 363796*MWs*powint(MZs,4) - 76447*powint(MZs,5) + 641560*sqr(MZs)*trip(MWs) - 698349*sqr(MWs)*trip(MZs))))/sqr(-4*MTs + MZs) + (MWs - MZs)*(8192*MTs*(-128*MWs*MZs + 512*sqr(MWs) - 9*sqr(MZs))*sqr(MZs)*trip(MWs)*trip(MWs - MZs) + 524288*powint(MWs,4)*powint(MWs - MZs,4)*trip(MZs) + 96*MZs*sqr(MTs)*sqr(MWs)*sqr(MWs - MZs)*(114688*powint(MWs,4) - 1941*powint(MZs,4) + 20992*sqr(MWs)*sqr(MZs) - 57344*MZs*trip(MWs) + 3648*MWs*trip(MZs)) + 9*MWs*trip(MTs)*(-1835008*MZs*powint(MWs,6) + 1048576*powint(MWs,7) - 6999*MWs*powint(MZs,6) - 7641*powint(MZs,7) + 140784*powint(MZs,5)*sqr(MWs) + 1425408*powint(MWs,5)*sqr(MZs) - 195776*powint(MZs,4)*trip(MWs) - 569344*powint(MWs,4)*trip(MZs)) + 81*powint(MTs,4)*sqr(MZs)*(-20480*MZs*powint(MWs,4) + 32768*powint(MWs,5) - 414*MWs*powint(MZs,4) - 81*powint(MZs,5) + 7168*sqr(MZs)*trip(MWs) + 320*sqr(MWs)*trip(MZs)))))/(4.478976e6*MTs) + real(B0fin(MZs,0,0))*(-(deltaAlpha*powint(MWs - MZs,-4)*sqr(AL)*(32*(40*MWs - 19*MZs)*MZs*sqr(MWs)*sqr(MWs - MZs) + 9*sqr(MTs)*sqr(MZs)*(-77*MWs*MZs + 100*sqr(MWs) + 19*sqr(MZs)) + MTs*MWs*(2560*powint(MWs,4) + 171*powint(MZs,4) + 2064*sqr(MWs)*sqr(MZs) - 3936*MZs*trip(MWs) - 859*MWs*trip(MZs))))/(288.*MTs*sqr(MWs)*sqr(Pi)) + (powint(MWs - MZs,-6)*powint(Pi,-3)*real(B0fin(MZs,MTs,MTs))*trip(AL)*(64*sqr(MWs)*sqr(MWs - MZs)*sqr(MZs)*(-984*MZs*sqr(MWs) + 696*MWs*sqr(MZs) + 640*trip(MWs) - 163*trip(MZs)) + 9*sqr(MZs)*trip(MTs)*(6400*powint(MWs,4) + 155*powint(MZs,4) + 4816*sqr(MWs)*sqr(MZs) - 11008*MZs*trip(MWs) - 1308*MWs*trip(MZs)) + sqr(MTs)*(-425984*MZs*powint(MWs,6) + 163840*powint(MWs,7) - 16641*MWs*powint(MZs,6) + 2925*powint(MZs,7) + 18317*powint(MZs,5)*sqr(MWs) + 408832*powint(MWs,5)*sqr(MZs) + 79808*powint(MZs,4)*trip(MWs) - 222592*powint(MWs,4)*trip(MZs)) + MTs*MWs*MZs*(-502784*MZs*powint(MWs,5) + 163840*powint(MWs,6) - 27469*MWs*powint(MZs,5) + 2925*powint(MZs,6) + 126368*powint(MZs,4)*sqr(MWs) + 627328*powint(MWs,4)*sqr(MZs) - 390208*trip(MWs)*trip(MZs))))/(20736.*MTs*MZs*sqr(MWs)) + (powint(MWs,-3)*powint(MWs - MZs,-5)*powint(Pi,-3)*trip(AL)*(4096*(40*MWs - 19*MZs)*sqr(MZs)*trip(MWs)*trip(MWs - MZs) + 32*MTs*MZs*sqr(MWs)*sqr(MWs - MZs)*(-13760*MZs*sqr(MWs) + 3640*MWs*sqr(MZs) + 25600*trip(MWs) - 171*trip(MZs)) + 27*sqr(MZs)*trip(MTs)*(12800*powint(MWs,4) - 171*powint(MZs,4) + 2936*sqr(MWs)*sqr(MZs) - 9216*MZs*trip(MWs) + 56*MWs*trip(MZs)) + 3*MWs*sqr(MTs)*(-524288*MZs*powint(MWs,5) + 327680*powint(MWs,6) + 38763*MWs*powint(MZs,5) - 8835*powint(MZs,6) - 29416*powint(MZs,4)*sqr(MWs) + 339456*powint(MWs,4)*sqr(MZs) - 143360*trip(MWs)*trip(MZs))))/(124416.*MTs*MZs)) + real(B0fin(MWs,0,0))*((deltaAlpha*powint(MWs - MZs,-4)*sqr(AL)*(32*(2*MWs - MZs)*MZs*sqr(MWs)*sqr(MWs - MZs) + 9*sqr(MTs)*sqr(MZs)*(-4*MWs*MZs + 5*sqr(MWs) + sqr(MZs)) + MTs*MWs*(128*powint(MWs,4) + 9*powint(MZs,4) + 106*sqr(MWs)*sqr(MZs) - 200*MZs*trip(MWs) - 43*MWs*trip(MZs))))/(16.*MTs*sqr(MWs)*sqr(Pi)) + (MZs*powint(MWs,-4)*powint(MWs - MZs,-6)*powint(Pi,-3)*real(B0fin(MWs,MTs,0))*trip(AL)*(64*MZs*powint(MWs,4)*sqr(MWs - MZs)*(-4*MWs*MZs + 5*sqr(MWs) + sqr(MZs)) + 9*powint(MTs,4)*sqr(MZs)*(14*MZs*sqr(MWs) - 6*MWs*sqr(MZs) - 14*trip(MWs) + trip(MZs)) + sqr(MTs)*sqr(MWs)*(400*MZs*powint(MWs,4) - 320*powint(MWs,5) + 239*MWs*powint(MZs,4) - 41*powint(MZs,5) + 322*sqr(MZs)*trip(MWs) - 510*sqr(MWs)*trip(MZs)) + 2*MTs*trip(MWs)*(-640*MZs*powint(MWs,4) + 320*powint(MWs,5) + 157*MWs*powint(MZs,4) - 25*powint(MZs,5) + 602*sqr(MZs)*trip(MWs) - 414*sqr(MWs)*trip(MZs)) + MWs*trip(MTs)*(560*MZs*powint(MWs,4) - 320*powint(MWs,5) - 115*MWs*powint(MZs,4) + 18*powint(MZs,5) - 504*sqr(MZs)*trip(MWs) + 316*sqr(MWs)*trip(MZs))))/(128.*MTs) + (MZs*powint(MWs - MZs,-6)*powint(Pi,-3)*real(B0fin(MZs,0,0))*trip(AL)*(-32*MZs*sqr(MWs)*sqr(MWs - MZs)*(-77*MWs*MZs + 100*sqr(MWs) + 19*sqr(MZs)) + 3*sqr(MTs)*sqr(MZs)*(812*MZs*sqr(MWs) - 344*MWs*sqr(MZs) - 840*trip(MWs) + 57*trip(MZs)) + MTs*MWs*(11008*MZs*powint(MWs,4) - 6400*powint(MWs,5) - 1187*MWs*powint(MZs,4) + 171*powint(MZs,5) - 7336*sqr(MZs)*trip(MWs) + 3744*sqr(MWs)*trip(MZs))))/(384.*MTs*sqr(MWs)) - (powint(MWs - MZs,-6)*powint(Pi,-3)*real(B0fin(MZs,MTs,MTs))*trip(AL)*(32*sqr(MWs)*sqr(MWs - MZs)*sqr(MZs)*(-100*MZs*sqr(MWs) + 71*MWs*sqr(MZs) + 64*trip(MWs) - 17*trip(MZs)) + 9*sqr(MZs)*trip(MTs)*(320*powint(MWs,4) + 7*powint(MZs,4) + 252*sqr(MWs)*sqr(MZs) - 560*MZs*trip(MWs) - 64*MWs*trip(MZs)) + sqr(MTs)*(-21504*MZs*powint(MWs,6) + 8192*powint(MWs,7) - 873*MWs*powint(MZs,6) + 153*powint(MZs,7) + 917*powint(MZs,5)*sqr(MWs) + 20800*powint(MWs,5)*sqr(MZs) + 3936*powint(MZs,4)*trip(MWs) - 11216*powint(MWs,4)*trip(MZs)) + MTs*MWs*MZs*(-25344*MZs*powint(MWs,5) + 8192*powint(MWs,6) - 1345*MWs*powint(MZs,5) + 153*powint(MZs,6) + 6384*powint(MZs,4)*sqr(MWs) + 31840*powint(MWs,4)*sqr(MZs) - 19880*trip(MWs)*trip(MZs))))/(1152.*MTs*MZs*sqr(MWs)) - (powint(MWs,-3)*powint(MWs - MZs,-5)*powint(Pi,-3)*trip(AL)*(4096*(2*MWs - MZs)*sqr(MZs)*trip(MWs)*trip(MWs - MZs) + 32*MTs*MZs*sqr(MWs)*sqr(MWs - MZs)*(-720*MZs*sqr(MWs) + 178*MWs*sqr(MZs) + 1280*trip(MWs) - 9*trip(MZs)) + 27*sqr(MZs)*trip(MTs)*(640*powint(MWs,4) - 9*powint(MZs,4) + 150*sqr(MWs)*sqr(MZs) - 480*MZs*trip(MWs) + 4*MWs*trip(MZs)) + 3*MWs*sqr(MTs)*(-26624*MZs*powint(MWs,5) + 16384*powint(MWs,6) + 2019*MWs*powint(MZs,5) - 465*powint(MZs,6) - 1474*powint(MZs,4)*sqr(MWs) + 17280*powint(MWs,4)*sqr(MZs) - 7120*trip(MWs)*trip(MZs))))/(6912.*MTs*MZs)) + real(B0fin(MWs,MTs,0))*(-(deltaAlpha*powint(MWs,-4)*powint(MWs - MZs,-4)*sqr(AL)*(MTs*MWs + sqr(MTs) - 2*sqr(MWs))*(32*(2*MWs - MZs)*MZs*sqr(MWs)*sqr(MWs - MZs) + 9*sqr(MTs)*sqr(MZs)*(-4*MWs*MZs + 5*sqr(MWs) + sqr(MZs)) + MTs*MWs*(128*powint(MWs,4) + 9*powint(MZs,4) + 106*sqr(MWs)*sqr(MZs) - 200*MZs*trip(MWs) - 43*MWs*trip(MZs))))/(96.*MTs*sqr(Pi)) + (MZs*powint(MWs,-4)*powint(MWs - MZs,-6)*powint(Pi,-3)*real(B0fin(MZs,0,0))*trip(AL)*(-64*MZs*powint(MWs,4)*sqr(MWs - MZs)*(-77*MWs*MZs + 100*sqr(MWs) + 19*sqr(MZs)) + 3*powint(MTs,4)*sqr(MZs)*(-812*MZs*sqr(MWs) + 344*MWs*sqr(MZs) + 840*trip(MWs) - 57*trip(MZs)) - 2*MTs*trip(MWs)*(-12608*MZs*powint(MWs,4) + 6400*powint(MWs,5) + 3027*MWs*powint(MZs,4) - 475*powint(MZs,5) + 11768*sqr(MZs)*trip(MWs) - 8112*sqr(MWs)*trip(MZs)) + MWs*trip(MTs)*(-11008*MZs*powint(MWs,4) + 6400*powint(MWs,5) + 2219*MWs*powint(MZs,4) - 342*powint(MZs,5) + 9856*sqr(MZs)*trip(MWs) - 6180*sqr(MWs)*trip(MZs)) + sqr(MTs)*sqr(MWs)*(-7808*MZs*powint(MWs,4) + 6400*powint(MWs,5) - 4557*MWs*powint(MZs,4) + 779*powint(MZs,5) - 6568*sqr(MZs)*trip(MWs) + 9864*sqr(MWs)*trip(MZs))))/(2304.*MTs) + (powint(MWs,-4)*powint(MWs - MZs,-6)*powint(Pi,-3)*real(B0fin(MZs,MTs,MTs))*trip(AL)*(-64*powint(MWs,4)*sqr(MWs - MZs)*sqr(MZs)*(-100*MZs*sqr(MWs) + 71*MWs*sqr(MZs) + 64*trip(MWs) - 17*trip(MZs)) + 9*powint(MTs,5)*sqr(MZs)*(320*powint(MWs,4) + 7*powint(MZs,4) + 252*sqr(MWs)*sqr(MZs) - 560*MZs*trip(MWs) - 64*MWs*trip(MZs)) - sqr(MTs)*sqr(MWs)*(-51200*MZs*powint(MWs,6) + 16384*powint(MWs,7) - 3761*MWs*powint(MZs,6) + 697*powint(MZs,7) + 3738*powint(MZs,5)*sqr(MWs) + 64896*powint(MWs,5)*sqr(MZs) + 17032*powint(MZs,4)*trip(MWs) - 46976*powint(MWs,4)*trip(MZs)) + powint(MTs,4)*(-21504*MZs*powint(MWs,6) + 8192*powint(MWs,7) - 810*MWs*powint(MZs,6) + 153*powint(MZs,7) + 341*powint(MZs,5)*sqr(MWs) + 23680*powint(MWs,5)*sqr(MZs) + 6204*powint(MZs,4)*trip(MWs) - 16256*powint(MWs,4)*trip(MZs)) + MWs*trip(MTs)*(-13312*MZs*powint(MWs,6) + 8192*powint(MWs,7) - 2344*MWs*powint(MZs,6) + 306*powint(MZs,7) + 8453*powint(MZs,5)*sqr(MWs) - 10304*powint(MWs,5)*sqr(MZs) - 20480*powint(MZs,4)*trip(MWs) + 30704*powint(MWs,4)*trip(MZs)) - 2*MTs*MZs*trip(MWs)*(-26368*MZs*powint(MWs,5) + 8192*powint(MWs,6) - 3025*MWs*powint(MZs,5) + 425*powint(MZs,6) + 10528*powint(MZs,4)*sqr(MWs) + 35488*powint(MWs,4)*sqr(MZs) - 25240*trip(MWs)*trip(MZs))))/(6912.*MTs*MZs) + (powint(MWs,-5)*powint(MWs - MZs,-5)*powint(Pi,-3)*(MTs*MWs + sqr(MTs) - 2*sqr(MWs))*trip(AL)*(4096*(2*MWs - MZs)*sqr(MZs)*trip(MWs)*trip(MWs - MZs) + 32*MTs*MZs*sqr(MWs)*sqr(MWs - MZs)*(-720*MZs*sqr(MWs) + 178*MWs*sqr(MZs) + 1280*trip(MWs) - 9*trip(MZs)) + 27*sqr(MZs)*trip(MTs)*(640*powint(MWs,4) - 9*powint(MZs,4) + 150*sqr(MWs)*sqr(MZs) - 480*MZs*trip(MWs) + 4*MWs*trip(MZs)) + 3*MWs*sqr(MTs)*(-26624*MZs*powint(MWs,5) + 16384*powint(MWs,6) + 2019*MWs*powint(MZs,5) - 465*powint(MZs,6) - 1474*powint(MZs,4)*sqr(MWs) + 17280*powint(MWs,4)*sqr(MZs) - 7120*trip(MWs)*trip(MZs))))/(41472.*MTs*MZs)) + real(B0fin(MZs,MTs,MTs))*((powint(MWs,-3)*powint(MWs - MZs,-5)*powint(MZs,-3)*powint(Pi,-3)*trip(AL)*(4096*(-40*MWs*MZs + 32*sqr(MWs) + 17*sqr(MZs))*trip(MWs)*trip(MWs - MZs)*trip(MZs) + 32*MTs*sqr(MWs)*sqr(MWs - MZs)*sqr(MZs)*(28672*powint(MWs,4) + 153*powint(MZs,4) + 31584*sqr(MWs)*sqr(MZs) - 49152*MZs*trip(MWs) - 4696*MWs*trip(MZs)) + MWs*MZs*sqr(MTs)*(-6750208*MZs*powint(MWs,6) + 2097152*powint(MWs,7) - 104571*MWs*powint(MZs,6) + 23715*powint(MZs,7) - 271400*powint(MZs,5)*sqr(MWs) + 8720384*powint(MWs,5)*sqr(MZs) + 1948672*powint(MZs,4)*trip(MWs) - 5663744*powint(MWs,4)*trip(MZs)) + 3*trip(MTs)*(-1441792*MZs*powint(MWs,7) + 524288*powint(MWs,8) + 283648*powint(MWs,4)*powint(MZs,4) + 3183*MWs*powint(MZs,7) + 1377*powint(MZs,8) - 56415*powint(MZs,6)*sqr(MWs) + 1593344*powint(MWs,6)*sqr(MZs) + 126736*powint(MZs,5)*trip(MWs) - 1009664*powint(MWs,5)*trip(MZs)) + 27*powint(MTs,4)*sqr(MZs)*(-26624*MZs*powint(MWs,4) + 16384*powint(MWs,5) + 632*MWs*powint(MZs,4) + 63*powint(MZs,5) + 11520*sqr(MZs)*trip(MWs) - 4720*sqr(MWs)*trip(MZs))))/(373248.*MTs) - (deltaAlpha*powint(MWs - MZs,-4)*sqr(AL)*(32*sqr(MWs)*sqr(MWs - MZs)*sqr(MZs)*(-40*MWs*MZs + 32*sqr(MWs) + 17*sqr(MZs)) + 9*sqr(MZs)*trip(MTs)*(-200*MZs*sqr(MWs) + 61*MWs*sqr(MZs) + 128*trip(MWs) - 7*trip(MZs)) + MTs*MWs*MZs*(-11776*MZs*powint(MWs,4) + 4096*powint(MWs,5) + 1081*MWs*powint(MZs,4) - 153*powint(MZs,5) + 12768*sqr(MZs)*trip(MWs) - 6016*sqr(MWs)*trip(MZs)) + sqr(MTs)*(-10240*MZs*powint(MWs,5) + 4096*powint(MWs,6) + 576*MWs*powint(MZs,5) - 153*powint(MZs,6) + 787*powint(MZs,4)*sqr(MWs) + 9024*powint(MWs,4)*sqr(MZs) - 3928*trip(MWs)*trip(MZs))))/(864.*MTs*sqr(MWs)*sqr(MZs)*sqr(Pi)))) + (powint(MTs,-3)*powint(MWs,-6)*powint(MWs - MZs,-6)*powint(MZs,-3)*powint(Pi,-3)*trip(AL)*(3072*MTs*powint(MWs,5)*powint(MWs - MZs,5)*sqr(MZs)*(-16*MWs*MZs + 64*sqr(MWs) + 9*sqr(MZs)) + 32768*powint(MWs,6)*powint(MWs - MZs,6)*trip(MZs) - 729*powint(MTs,6)*powint(MZs,6)*(14*MZs*sqr(MWs) - 6*MWs*sqr(MZs) - 14*trip(MWs) + trip(MZs)) + 96*MZs*powint(MWs,4)*powint(MWs - MZs,4)*sqr(MTs)*(4096*powint(MWs,4) - 207*powint(MZs,4) + 1984*sqr(MWs)*sqr(MZs) - 2048*MZs*trip(MWs) + 288*MWs*trip(MZs)) + 27*powint(MTs,4)*sqr(MWs)*sqr(MWs - MZs)*sqr(MZs)*(-5120*MZs*powint(MWs,4) + 8192*powint(MWs,5) - 702*MWs*powint(MZs,4) + 207*powint(MZs,5) + 5248*sqr(MZs)*trip(MWs) - 640*sqr(MWs)*trip(MZs)) + 243*MWs*powint(MTs,5)*powint(MZs,4)*(-560*MZs*powint(MWs,4) + 320*powint(MWs,5) + 61*MWs*powint(MZs,4) - 9*powint(MZs,5) + 378*sqr(MZs)*trip(MWs) - 190*sqr(MWs)*trip(MZs)) + trip(MTs)*trip(MWs)*trip(MWs - MZs)*(-196608*MZs*powint(MWs,5) + 262144*powint(MWs,6) + 54864*MWs*powint(MZs,5) - 14823*powint(MZs,6) - 91584*powint(MZs,4)*sqr(MWs) + 270336*powint(MWs,4)*sqr(MZs) + 134144*trip(MWs)*trip(MZs)))*trip(A0fin(MTs)))/373248. - (27*powint(MWs - MZs,-6)*powint(Pi,-3)*trip(AL)*trip(MZs)*(14*MZs*sqr(MWs) - 6*MWs*sqr(MZs) - 14*trip(MWs) + trip(MZs))*trip(real(B0fin(MWs,0,0))))/64. + (powint(MWs,-6)*powint(MWs - MZs,-6)*powint(Pi,-3)*trip(AL)*trip(-MTs + MWs)*trip(MTs + 2*MWs)*(-14*MZs*sqr(MWs) + 6*MWs*sqr(MZs) + 14*trip(MWs) - trip(MZs))*trip(MZs)*trip(real(B0fin(MWs,MTs,0))))/512. + (35*powint(MWs - MZs,-6)*powint(Pi,-3)*trip(AL)*trip(MZs)*(2880*MZs*sqr(MWs) - 1200*MWs*sqr(MZs) - 3200*trip(MWs) + 197*trip(MZs))*trip(real(B0fin(MZs,0,0))))/13824. - (powint(MWs - MZs,-6)*powint(MZs,-3)*powint(Pi,-3)*trip(AL)*(trip(MTs)*(-983040*MZs*powint(MWs,5) + 262144*powint(MWs,6) - 292128*MWs*powint(MZs,5) + 13411*powint(MZs,6) + 739968*powint(MZs,4)*sqr(MWs) + 1425408*powint(MWs,4)*sqr(MZs) - 1169408*trip(MWs)*trip(MZs)) + 3*MZs*sqr(MTs)*(-491520*MZs*powint(MWs,5) + 131072*powint(MWs,6) - 115824*MWs*powint(MZs,5) + 23405*powint(MZs,6) + 369984*powint(MZs,4)*sqr(MWs) + 768000*powint(MWs,4)*sqr(MZs) - 681472*trip(MWs)*trip(MZs)) + 3*MTs*sqr(MZs)*(-245760*MZs*powint(MWs,5) + 65536*powint(MWs,6) - 60288*MWs*powint(MZs,5) + 6019*powint(MZs,6) + 208320*powint(MZs,4)*sqr(MWs) + 411648*powint(MWs,4)*sqr(MZs) - 389120*trip(MWs)*trip(MZs)) + trip(MZs)*(-122880*MZs*powint(MWs,5) + 32768*powint(MWs,6) - 40080*MWs*powint(MZs,5) + 5453*powint(MZs,6) + 127488*powint(MZs,4)*sqr(MWs) + 219648*powint(MWs,4)*sqr(MZs) - 218752*trip(MWs)*trip(MZs)))*trip(real(B0fin(MZs,MTs,MTs))))/373248.));
}