PR121:=(mm1,mm2,mm)-> evalf( int( ln( mm2/mm*x^2 + mm1/mm*(1-x) -I*10^(-200)),x=0..1,`continuous`) ); PR11:=(mm1,mm2,mm)-> -evalf( int( ln( mm2/mm*x*(x-1) + mm1/mm -I*10^(-200)),x=0..1,`continuous`) ); f:=(mm1,mm2,a)-> evalf( 1/(mm1-mm2)^a); eP:=137.036; alpha:=1/(4*Pi*eP); mmW:=(80.419)^2; mmZ:=(91.188)^2; ########################################################################################## # # m_MS/M_pole = EXACT_1 # EXACT1_Z_B:=(mmH,mm)-> evalf( + mmH^2/mmW/mmZ * ( - 1/12 + 1/12*ln(mmH/mm) - 1/12*PR121(mmH,mmZ,mm) ) + mmZ/mmW * ( - 2/9 + 1/6*ln(mmZ/mm) + 1/12*PR11(mmW,mmZ,mm) - PR121(mmH,mmZ,mm) ) + mmZ^2/mmH/mmW * ( - 1/2 + 3/2*ln(mmZ/mm) ) + mmH/mmW * ( - 7/12 - 1/12*ln(mmZ/mm) + 1/2*ln(mmH/mm) + 1/3*PR121(mmH,mmZ,mm) ) + mmW/mmZ * ( 4/3*ln(mmW/mm) - 17/3*PR11(mmW,mmZ,mm) ) + mmW/mmH * ( - 1 + 3*ln(mmW/mm) ) + mmW^2/mmZ^2 * ( - 4*ln(mmW/mm) - 4*PR11(mmW,mmZ,mm) ) + ( - 13/18 + 1/6*ln(mmW/mm) + 4/3*PR11(mmW,mmZ,mm) ) ); EXACT1_W_B :=(mmH,mm)-> evalf( + mmZ^2/mmW^2 * ( - 1/12 + 1/12*ln(mmZ/mm) - 1/12*PR121(mmZ,mmW,mm) ) + mmH^2/mmW^2 * ( - 1/12 + 1/12*ln(mmH/mm) - 1/12*PR121(mmH,mmW,mm) ) + mmZ/mmW * ( - 3/4 - 1/12*ln(mmW/mm) + 2/3*ln(mmZ/mm) - 4/3*PR121(mmZ,mmW,mm) ) + mmZ^2/mmH/mmW * ( - 1/2 + 3/2*ln(mmZ/mm) ) + mmH/mmW * ( - 7/12 - 1/12*ln(mmW/mm) + 1/2*ln(mmH/mm) + 1/3*PR121(mmH,mmW,mm) ) + mmW/mmZ * ( 8 - 4*ln(mmW/mm) + 4*PR121(mmZ,mmW,mm) ) + mmW/mmH * ( - 1 + 3*ln(mmW/mm) ) + ( - 73/9 + 3*ln(mmW/mm) - 2*ln(mmZ/mm) + 17/3*PR121(mmZ,mmW,mm) - PR121(mmH,mmW,mm) ) ); EXACT1_1_H :=(mmH,mm)-> evalf( + mmZ/mmW * ( - 1/4 + 1/4*ln(mmZ/mm) - 1/2*PR11(mmZ,mmH,mm) ) + mmZ^2/mmH/mmW * ( - 3/2 + 3/2*ln(mmZ/mm) + 3/2*PR11(mmZ,mmH,mm) ) + mmH/mmW * ( 3/2 - 9/8*Pi/sqrt(3) - 3/8*ln(mmH/mm) + 1/4*PR11(mmW,mmH,mm) + 1/8*PR11(mmZ,mmH,mm) ) + mmW/mmH * ( - 3 + 3*ln(mmW/mm) + 3*PR11(mmW,mmH,mm) ) - 1/2 + 1/2*ln(mmW/mm) - PR11(mmW,mmH,mm) );