Pastebin launched a little side project called VERYVIRAL.com, check it out ;-) Want more features on Pastebin? Sign Up, it's FREE!
Guest

Untitled

By: a guest on Nov 18th, 2010  |  syntax: Fortran  |  size: 11.03 KB  |  views: 36  |  expires: Never
download  |  raw  |  embed  |  report abuse  |  print
Text below is selected. Please press Ctrl+C to copy to your clipboard. (⌘+C on Mac)
  1.        MX11=(4.d0*mH1**2+3.d0*g2**2*v1**2-g2**2*v2**2+4.d0*lambda**2*
  2.      . v2**2+g1**2*(3.d0*v1**2-v2**2)+4.d0*lambda**2*vevS**2)/4.d0+
  3.      . (3.d0*(-28.d0*v1**2*Ybtm**4+(24.d0*lambda**2*vevS**2*Ytop**4*
  4.      . (lambda*v1*vevS-Atop*v2*DCOS(phi2+phiS))**2)/(MQ33**4-2.d0*
  5.      . MQ33**2*MU33**2+MU33**4+4.d0*Atop**2*v2**2*Ytop**2+4.d0*
  6.      . lambda**2*v1**2*vevS**2*Ytop**2-8.d0*Atop*lambda*v1*v2*vevS*
  7.      . Ytop**2*DCOS(phi2+phiS))-12.d0*v1**2*Ybtm**4*(-3.d0+2.d0*
  8.      . DLOG((v1**2*Ybtm**2)/Rsc**2))+(4.d0*lambda**2*vevS**2*Ytop**4*
  9.      . (lambda*v1*vevS-Atop*v2*DCOS(phi2+phiS))**2*(-3.d0+2.d0*
  10.      . DLOG((MQ33**2+MU33**2+2.d0*v2**2*Ytop**2-DSQRT(MQ33**4-2.d0*
  11.      . MQ33**2*MU33**2+MU33**4+4.d0*Atop**2*v2**2*Ytop**2+4.d0*
  12.      . lambda**2*v1**2*vevS**2*Ytop**2-8.d0*Atop*lambda*v1*v2*vevS*
  13.      . Ytop**2*DCOS(phi2+phiS)))/(2.d0*Rsc**2))))/(MQ33**4-2.d0*MQ33**2
  14.      . *MU33**2+MU33**4+4.d0*Atop**2*v2**2*Ytop**2+4.d0*lambda**2*v1**2
  15.      . *vevS**2*Ytop**2-8.d0*Atop*lambda*v1*v2*vevS*Ytop**2*DCOS(phi2+
  16.      . phiS))+(4.d0*lambda**2*vevS**2*Ytop**4*(lambda*v1*vevS-Atop*v2*
  17.      . DCOS(phi2+phiS))**2*(-3.d0+2.d0*DLOG((MQ33**2+MU33**2+2.d0*v2**2
  18.      . *Ytop**2+DSQRT(MQ33**4-2.d0*MQ33**2*MU33**2+MU33**4+4.d0*Atop**2
  19.      . *v2**2*Ytop**2+4.d0*lambda**2*v1**2*vevS**2*Ytop**2-8.d0*Atop*
  20.      . lambda*v1*v2*vevS*Ytop**2*DCOS(phi2+phiS)))/(2.d0*Rsc**2))))/
  21.      . (MQ33**4-2.d0*MQ33**2*MU33**2+MU33**4+4.d0*Atop**2*v2**2*Ytop**2
  22.      . +4.d0*lambda**2*v1**2*vevS**2*Ytop**2-8.d0*Atop*lambda*v1*v2*
  23.      . vevS*Ytop**2*DCOS(phi2+phiS))+12.d0*Ybtm**4*(v1+(Abtm*(Abtm*v1-
  24.      . lambda*v2*vevS*DCOS(phi2+phiS)))/DSQRT(MD33**4-2.d0*MD33**2*
  25.      . MQ33**2+MQ33**4+4.d0*Abtm**2*v1**2*Ybtm**2+4.d0*lambda**2*v2**2*
  26.      . vevS**2*Ybtm**2-8.d0*Abtm*lambda*v1*v2*vevS*Ybtm**2*DCOS(phi2+
  27.      . phiS)))**2+8.d0*Ybtm**4*(-1.5d0+DLOG((MD33**2+MQ33**2+2.d0*v1**2
  28.      . *Ybtm**2+DSQRT(MD33**4-2.d0*MD33**2*MQ33**2+MQ33**4+4.d0*Abtm**2
  29.      . *v1**2*Ybtm**2+4.d0*lambda**2*v2**2*vevS**2*Ybtm**2-8.d0*Abtm*
  30.      . lambda*v1*v2*vevS*Ybtm**2*DCOS(phi2+phiS)))/(2.d0*Rsc**2)))*(v1+
  31.      . (Abtm*(Abtm*v1-lambda*v2*vevS*DCOS(phi2+phiS)))/DSQRT(MD33**4-
  32.      . 2.d0*MD33**2*MQ33**2+MQ33**4+4.d0*Abtm**2*v1**2*Ybtm**2+4.d0*
  33.      . lambda**2*v2**2*vevS**2*Ybtm**2-8.d0*Abtm*lambda*v1*v2*vevS*
  34.      . Ybtm**2*DCOS(phi2+phiS)))**2+12.d0*Ybtm**4*(v1+(Abtm*(-(Abtm*v1)
  35.      . +lambda*v2*vevS*DCOS(phi2+phiS)))/DSQRT(MD33**4-2.d0*MD33**2*
  36.      . MQ33**2+MQ33**4+4.d0*Abtm**2*v1**2*Ybtm**2+4.d0*lambda**2*v2**2*
  37.      . vevS**2*Ybtm**2-8.d0*Abtm*lambda*v1*v2*vevS*Ybtm**2*DCOS(phi2+
  38.      . phiS)))**2+8.d0*Ybtm**4*(-1.5d0+DLOG((MD33**2+MQ33**2+2.d0*v1**2
  39.      . *Ybtm**2-DSQRT(MD33**4-2.d0*MD33**2*MQ33**2+MQ33**4+4.d0*Abtm**2
  40.      . *v1**2*Ybtm**2+4.d0*lambda**2*v2**2*vevS**2*Ybtm**2-8.d0*Abtm*
  41.      . lambda*v1*v2*vevS*Ybtm**2*DCOS(phi2+phiS)))/(2.d0*Rsc**2)))*(v1+
  42.      . (Abtm*(-(Abtm*v1)+lambda*v2*vevS*DCOS(phi2+phiS)))/DSQRT(MD33**4
  43.      . -2.d0*MD33**2*MQ33**2+MQ33**4+4.d0*Abtm**2*v1**2*Ybtm**2+4.d0*
  44.      . lambda**2*v2**2*vevS**2*Ybtm**2-8.d0*Abtm*lambda*v1*v2*vevS*
  45.      . Ybtm**2*DCOS(phi2+phiS)))**2+Ybtm**2*(1.d0+(4.d0*Abtm**2*Ybtm**2
  46.      . *(Abtm*v1-lambda*v2*vevS*DCOS(phi2+phiS))**2)/(MD33**4-2.d0*
  47.      . MD33**2*MQ33**2+MQ33**4+4.d0*Abtm**2*v1**2*Ybtm**2+4.d0*
  48.      . lambda**2*v2**2*vevS**2*Ybtm**2-8.d0*Abtm*lambda*v1*v2*vevS*
  49.      . Ybtm**2*DCOS(phi2+phiS))**1.5-Abtm**2/DSQRT(MD33**4-2.d0*
  50.      . MD33**2*MQ33**2+MQ33**4+4.d0*Abtm**2*v1**2*Ybtm**2+4.d0*
  51.      . lambda**2*v2**2*vevS**2*Ybtm**2-8.d0*Abtm*lambda*v1*v2*vevS*
  52.      . Ybtm**2*DCOS(phi2+phiS)))*(MD33**2+MQ33**2+2.d0*v1**2*Ybtm**2-
  53.      . DSQRT(MD33**4-2.d0*MD33**2*MQ33**2+MQ33**4+4.d0*Abtm**2*v1**2*
  54.      . Ybtm**2+4.d0*lambda**2*v2**2*vevS**2*Ybtm**2-8.d0*Abtm*lambda*v1
  55.      . *v2*vevS*Ybtm**2*DCOS(phi2+phiS)))+2.d0*Ybtm**2*(-1.5d0+
  56.      . DLOG((MD33**2+MQ33**2+2.d0*v1**2*Ybtm**2-DSQRT(MD33**4-2.d0*
  57.      . MD33**2*MQ33**2+MQ33**4+4.d0*Abtm**2*v1**2*Ybtm**2+4.d0*
  58.      . lambda**2*v2**2*vevS**2*Ybtm**2-8.d0*Abtm*lambda*v1*v2*vevS*
  59.      . Ybtm**2*DCOS(phi2+phiS)))/(2.d0*Rsc**2)))*(1.d0+(4.d0*Abtm**2*
  60.      . Ybtm**2*(Abtm*v1-lambda*v2*vevS*DCOS(phi2+phiS))**2)/(MD33**4-
  61.      . 2.d0*MD33**2*MQ33**2+MQ33**4+4.d0*Abtm**2*v1**2*Ybtm**2+4.d0*
  62.      . lambda**2*v2**2*vevS**2*Ybtm**2-8.d0*Abtm*lambda*v1*v2*vevS*
  63.      . Ybtm**2*DCOS(phi2+phiS))**1.5-Abtm**2/DSQRT(MD33**4-2.d0*
  64.      . MD33**2*MQ33**2+MQ33**4+4.d0*Abtm**2*v1**2*Ybtm**2+4.d0*
  65.      . lambda**2*v2**2*vevS**2*Ybtm**2-8.d0*Abtm*lambda*v1*v2*vevS*
  66.      . Ybtm**2*DCOS(phi2+phiS)))*(MD33**2+MQ33**2+2.d0*v1**2*Ybtm**2-
  67.      . DSQRT(MD33**4-2.d0*MD33**2*MQ33**2+MQ33**4+4.d0*Abtm**2*v1**2*
  68.      . Ybtm**2+4.d0*lambda**2*v2**2*vevS**2*Ybtm**2-8.d0*Abtm*lambda*v1
  69.      . *v2*vevS*Ybtm**2*DCOS(phi2+phiS)))+Ybtm**2*(1.d0-(4.d0*Abtm**2*
  70.      . Ybtm**2*(Abtm*v1-lambda*v2*vevS*DCOS(phi2+phiS))**2)/(MD33**4-
  71.      . 2.d0*MD33**2*MQ33**2+MQ33**4+4.d0*Abtm**2*v1**2*Ybtm**2+4.d0*
  72.      . lambda**2*v2**2*vevS**2*Ybtm**2-8.d0*Abtm*lambda*v1*v2*vevS*
  73.      . Ybtm**2*DCOS(phi2+phiS))**1.5+Abtm**2/DSQRT(MD33**4-2.d0*
  74.      . MD33**2*MQ33**2+MQ33**4+4.d0*Abtm**2*v1**2*Ybtm**2+4.d0*
  75.      . lambda**2*v2**2*vevS**2*Ybtm**2-8.d0*Abtm*lambda*v1*v2*vevS*
  76.      . Ybtm**2*DCOS(phi2+phiS)))*(MD33**2+MQ33**2+2.d0*v1**2*Ybtm**2+
  77.      . DSQRT(MD33**4-2.d0*MD33**2*MQ33**2+MQ33**4+4.d0*Abtm**2*v1**2*
  78.      . Ybtm**2+4.d0*lambda**2*v2**2*vevS**2*Ybtm**2-8.d0*Abtm*lambda*v1
  79.      . *v2*vevS*Ybtm**2*DCOS(phi2+phiS)))+2.d0*Ybtm**2*(-1.5d0+
  80.      . DLOG((MD33**2+MQ33**2+2.d0*v1**2*Ybtm**2+DSQRT(MD33**4-2.d0*
  81.      . MD33**2*MQ33**2+MQ33**4+4.d0*Abtm**2*v1**2*Ybtm**2+4.d0*
  82.      . lambda**2*v2**2*vevS**2*Ybtm**2-8.d0*Abtm*lambda*v1*v2*vevS*
  83.      . Ybtm**2*DCOS(phi2+phiS)))/(2.d0*Rsc**2)))*(1.d0-(4.d0*Abtm**2*
  84.      . Ybtm**2*(Abtm*v1-lambda*v2*vevS*DCOS(phi2+phiS))**2)/(MD33**4-
  85.      . 2.d0*MD33**2*MQ33**2+MQ33**4+4.d0*Abtm**2*v1**2*Ybtm**2+4.d0*
  86.      . lambda**2*v2**2*vevS**2*Ybtm**2-8.d0*Abtm*lambda*v1*v2*vevS*
  87.      . Ybtm**2*DCOS(phi2+phiS))**1.5+Abtm**2/DSQRT(MD33**4-2.d0*
  88.      . MD33**2*MQ33**2+MQ33**4+4.d0*Abtm**2*v1**2*Ybtm**2+4.d0*
  89.      . lambda**2*v2**2*vevS**2*Ybtm**2-8.d0*Abtm*lambda*v1*v2*vevS*
  90.      . Ybtm**2*DCOS(phi2+phiS)))*(MD33**2+MQ33**2+2.d0*v1**2*Ybtm**2+
  91.      . DSQRT(MD33**4-2.d0*MD33**2*MQ33**2+MQ33**4+4.d0*Abtm**2*v1**2*
  92.      . Ybtm**2+4.d0*lambda**2*v2**2*vevS**2*Ybtm**2-8.d0*Abtm*lambda*v1
  93.      . *v2*vevS*Ybtm**2*DCOS(phi2+phiS)))+(4.d0*lambda**2*vevS**2*
  94.      . Ytop**4*(lambda*v1*vevS-Atop*v2*DCOS(phi2+phiS))**2*(MQ33**2+
  95.      . MU33**2+2.d0*v2**2*Ytop**2-DSQRT(MQ33**4-2.d0*MQ33**2*MU33**2+
  96.      . MU33**4+4.d0*Atop**2*v2**2*Ytop**2+4.d0*lambda**2*v1**2*vevS**2*
  97.      . Ytop**2-8.d0*Atop*lambda*v1*v2*vevS*Ytop**2*DCOS(phi2+phiS))))/
  98.      . (MQ33**4-2.d0*MQ33**2*MU33**2+MU33**4+4.d0*Atop**2*v2**2*Ytop**2
  99.      . +4.d0*lambda**2*v1**2*vevS**2*Ytop**2-8.d0*Atop*lambda*v1*v2*
  100.      . vevS*Ytop**2*DCOS(phi2+phiS))**1.5+(4.d0*lambda**2*vevS**2*
  101.      . Ytop**4*(lambda*v1*vevS-Atop*v2*DCOS(phi2+phiS))**2*(-3.d0+2.d0*
  102.      . DLOG((MQ33**2+MU33**2+2.d0*v2**2*Ytop**2-DSQRT(MQ33**4-2.d0*
  103.      . MQ33**2*MU33**2+MU33**4+4.d0*Atop**2*v2**2*Ytop**2+4.d0*
  104.      . lambda**2*v1**2*vevS**2*Ytop**2-8.d0*Atop*lambda*v1*v2*vevS*
  105.      . Ytop**2*DCOS(phi2+phiS)))/(2.d0*Rsc**2)))*(MQ33**2+MU33**2+2.d0*
  106.      . v2**2*Ytop**2-DSQRT(MQ33**4-2.d0*MQ33**2*MU33**2+MU33**4+4.d0*
  107.      . Atop**2*v2**2*Ytop**2+4.d0*lambda**2*v1**2*vevS**2*Ytop**2-8.d0*
  108.      . Atop*lambda*v1*v2*vevS*Ytop**2*DCOS(phi2+phiS))))/(MQ33**4-2.d0*
  109.      . MQ33**2*MU33**2+MU33**4+4.d0*Atop**2*v2**2*Ytop**2+4.d0*
  110.      . lambda**2*v1**2*vevS**2*Ytop**2-8.d0*Atop*lambda*v1*v2*vevS*
  111.      . Ytop**2*DCOS(phi2+phiS))**1.5-(lambda**2*vevS**2*Ytop**2*
  112.      . (MQ33**2+MU33**2+2.d0*v2**2*Ytop**2-DSQRT(MQ33**4-2.d0*MQ33**2*
  113.      . MU33**2+MU33**4+4.d0*Atop**2*v2**2*Ytop**2+4.d0*lambda**2*v1**2*
  114.      . vevS**2*Ytop**2-8.d0*Atop*lambda*v1*v2*vevS*Ytop**2*DCOS(phi2+
  115.      . phiS))))/DSQRT(MQ33**4-2.d0*MQ33**2*MU33**2+MU33**4+4.d0*Atop**2
  116.      . *v2**2*Ytop**2+4.d0*lambda**2*v1**2*vevS**2*Ytop**2-8.d0*Atop*
  117.      . lambda*v1*v2*vevS*Ytop**2*DCOS(phi2+phiS))-(lambda**2*vevS**2*
  118.      . Ytop**2*(-3.d0+2.d0*DLOG((MQ33**2+MU33**2+2.d0*v2**2*Ytop**2-
  119.      . DSQRT(MQ33**4-2.d0*MQ33**2*MU33**2+MU33**4+4.d0*Atop**2*v2**2*
  120.      . Ytop**2+4.d0*lambda**2*v1**2*vevS**2*Ytop**2-8.d0*Atop*lambda*v1
  121.      . *v2*vevS*Ytop**2*DCOS(phi2+phiS)))/(2.d0*Rsc**2)))*(MQ33**2+
  122.      . MU33**2+2.d0*v2**2*Ytop**2-DSQRT(MQ33**4-2.d0*MQ33**2*MU33**2+
  123.      . MU33**4+4.d0*Atop**2*v2**2*Ytop**2+4.d0*lambda**2*v1**2*vevS**2*
  124.      . Ytop**2-8.d0*Atop*lambda*v1*v2*vevS*Ytop**2*DCOS(phi2+phiS))))/
  125.      . DSQRT(MQ33**4-2.d0*MQ33**2*MU33**2+MU33**4+4.d0*Atop**2*v2**2*
  126.      . Ytop**2+4.d0*lambda**2*v1**2*vevS**2*Ytop**2-8.d0*Atop*lambda*v1
  127.      . *v2*vevS*Ytop**2*DCOS(phi2+phiS))-(4.d0*lambda**2*vevS**2*
  128.      . Ytop**4*(lambda*v1*vevS-Atop*v2*DCOS(phi2+phiS))**2*(MQ33**2+
  129.      . MU33**2+2.d0*v2**2*Ytop**2+DSQRT(MQ33**4-2.d0*MQ33**2*MU33**2+
  130.      . MU33**4+4.d0*Atop**2*v2**2*Ytop**2+4.d0*lambda**2*v1**2*vevS**2*
  131.      . Ytop**2-8.d0*Atop*lambda*v1*v2*vevS*Ytop**2*DCOS(phi2+phiS))))/
  132.      . (MQ33**4-2.d0*MQ33**2*MU33**2+MU33**4+4.d0*Atop**2*v2**2*Ytop**2
  133.      . +4.d0*lambda**2*v1**2*vevS**2*Ytop**2-8.d0*Atop*lambda*v1*v2*
  134.      . vevS*Ytop**2*DCOS(phi2+phiS))**1.5-(4.d0*lambda**2*vevS**2*
  135.      . Ytop**4*(lambda*v1*vevS-Atop*v2*DCOS(phi2+phiS))**2*(-3.d0+2.d0*
  136.      . DLOG((MQ33**2+MU33**2+2.d0*v2**2*Ytop**2+DSQRT(MQ33**4-2.d0*
  137.      . MQ33**2*MU33**2+MU33**4+4.d0*Atop**2*v2**2*Ytop**2+4.d0*
  138.      . lambda**2*v1**2*vevS**2*Ytop**2-8.d0*Atop*lambda*v1*v2*vevS*
  139.      . Ytop**2*DCOS(phi2+phiS)))/(2.d0*Rsc**2)))*(MQ33**2+MU33**2+2.d0*
  140.      . v2**2*Ytop**2+DSQRT(MQ33**4-2.d0*MQ33**2*MU33**2+MU33**4+4.d0*
  141.      . Atop**2*v2**2*Ytop**2+4.d0*lambda**2*v1**2*vevS**2*Ytop**2-8.d0*
  142.      . Atop*lambda*v1*v2*vevS*Ytop**2*DCOS(phi2+phiS))))/(MQ33**4-2.d0*
  143.      . MQ33**2*MU33**2+MU33**4+4.d0*Atop**2*v2**2*Ytop**2+4.d0*
  144.      . lambda**2*v1**2*vevS**2*Ytop**2-8.d0*Atop*lambda*v1*v2*vevS*
  145.      . Ytop**2*DCOS(phi2+phiS))**1.5+(lambda**2*vevS**2*Ytop**2*
  146.      . (MQ33**2+MU33**2+2.d0*v2**2*Ytop**2+DSQRT(MQ33**4-2.d0*MQ33**2*
  147.      . MU33**2+MU33**4+4.d0*Atop**2*v2**2*Ytop**2+4.d0*lambda**2*v1**2*
  148.      . vevS**2*Ytop**2-8.d0*Atop*lambda*v1*v2*vevS*Ytop**2*DCOS(phi2+
  149.      . phiS))))/DSQRT(MQ33**4-2.d0*MQ33**2*MU33**2+MU33**4+4.d0*Atop**2
  150.      . *v2**2*Ytop**2+4.d0*lambda**2*v1**2*vevS**2*Ytop**2-8.d0*Atop*
  151.      . lambda*v1*v2*vevS*Ytop**2*DCOS(phi2+phiS))+(lambda**2*vevS**2*
  152.      . Ytop**2*(-3.d0+2.d0*DLOG((MQ33**2+MU33**2+2.d0*v2**2*Ytop**2+
  153.      . DSQRT(MQ33**4-2.d0*MQ33**2*MU33**2+MU33**4+4.d0*Atop**2*v2**2*
  154.      . Ytop**2+4.d0*lambda**2*v1**2*vevS**2*Ytop**2-8.d0*Atop*lambda*v1
  155.      . *v2*vevS*Ytop**2*DCOS(phi2+phiS)))/(2.d0*Rsc**2)))*(MQ33**2+
  156.      . MU33**2+2.d0*v2**2*Ytop**2+DSQRT(MQ33**4-2.d0*MQ33**2*MU33**2+
  157.      . MU33**4+4.d0*Atop**2*v2**2*Ytop**2+4.d0*lambda**2*v1**2*vevS**2*
  158.      . Ytop**2-8.d0*Atop*lambda*v1*v2*vevS*Ytop**2*DCOS(phi2+phiS))))/
  159.      . DSQRT(MQ33**4-2.d0*MQ33**2*MU33**2+MU33**4+4.d0*Atop**2*v2**2*
  160.      . Ytop**2+4.d0*lambda**2*v1**2*vevS**2*Ytop**2-8.d0*Atop*lambda*v1
  161.      . *v2*vevS*Ytop**2*DCOS(phi2+phiS))))/(64.d0*Pi**2)