Pastebin launched a little side project called HostCabi.net, check it out ;-)Don't like ads? PRO users don't see any ads ;-)
Guest

Untitled

By: a guest on Nov 18th, 2010  |  syntax: C  |  size: 39.72 KB  |  hits: 45  |  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*Power(mH1,2) + 3*Power(g2,2)*Power(v1,2) - Power(g2,2)*Power(v2,2) +
  2.       4*Power(lambda,2)*Power(v2,2) +
  3.       Power(g1,2)*(3*Power(v1,2) - Power(v2,2)) +
  4.       4*Power(lambda,2)*Power(vevS,2))/4. +
  5.    (3*((-3*Power(2*v1*Power(Ybtm,2) -
  6.              (2*v1*Power(Ybtm,2)*
  7.                 (Power(v1,2)*Power(Ybtm,2) -
  8.                   Power(Ytop,2)*
  9.                    (Power(v2,2)*Power(Cos(phi2),2) +
  10.                      Power(v2,2)*Power(Sin(phi2),2))))/
  11.               Sqrt(Power(Power(v1,2)*Power(Ybtm,2) -
  12.                  Power(Ytop,2)*(Power(v2,2)*Power(Cos(phi2),2) +
  13.                     Power(v2,2)*Power(Sin(phi2),2)),2)),2))/2. -
  14.         (-1.5 + Log((Power(v1,2)*Power(Ybtm,2) +
  15.                Power(v2,2)*Power(Ytop,2)*Power(Cos(phi2),2) +
  16.                Power(v2,2)*Power(Ytop,2)*Power(Sin(phi2),2) -
  17.                Sqrt(Power(Power(v1,2)*Power(Ybtm,2) -
  18.                   Power(Ytop,2)*
  19.                    (Power(v2,2)*Power(Cos(phi2),2) +
  20.                      Power(v2,2)*Power(Sin(phi2),2)),2)))/(2.*Power(Rsc,2))))*
  21.          Power(2*v1*Power(Ybtm,2) -
  22.            (2*v1*Power(Ybtm,2)*(Power(v1,2)*Power(Ybtm,2) -
  23.                 Power(Ytop,2)*(Power(v2,2)*Power(Cos(phi2),2) +
  24.                    Power(v2,2)*Power(Sin(phi2),2))))/
  25.             Sqrt(Power(Power(v1,2)*Power(Ybtm,2) -
  26.                Power(Ytop,2)*(Power(v2,2)*Power(Cos(phi2),2) +
  27.                   Power(v2,2)*Power(Sin(phi2),2)),2)),2) -
  28.         (3*Power(2*v1*Power(Ybtm,2) +
  29.              (2*v1*Power(Ybtm,2)*
  30.                 (Power(v1,2)*Power(Ybtm,2) -
  31.                   Power(Ytop,2)*
  32.                    (Power(v2,2)*Power(Cos(phi2),2) +
  33.                      Power(v2,2)*Power(Sin(phi2),2))))/
  34.               Sqrt(Power(Power(v1,2)*Power(Ybtm,2) -
  35.                  Power(Ytop,2)*(Power(v2,2)*Power(Cos(phi2),2) +
  36.                     Power(v2,2)*Power(Sin(phi2),2)),2)),2))/2. -
  37.         (-1.5 + Log((Power(v1,2)*Power(Ybtm,2) +
  38.                Power(v2,2)*Power(Ytop,2)*Power(Cos(phi2),2) +
  39.                Power(v2,2)*Power(Ytop,2)*Power(Sin(phi2),2) +
  40.                Sqrt(Power(Power(v1,2)*Power(Ybtm,2) -
  41.                   Power(Ytop,2)*
  42.                    (Power(v2,2)*Power(Cos(phi2),2) +
  43.                      Power(v2,2)*Power(Sin(phi2),2)),2)))/(2.*Power(Rsc,2))))*
  44.          Power(2*v1*Power(Ybtm,2) +
  45.            (2*v1*Power(Ybtm,2)*(Power(v1,2)*Power(Ybtm,2) -
  46.                 Power(Ytop,2)*(Power(v2,2)*Power(Cos(phi2),2) +
  47.                    Power(v2,2)*Power(Sin(phi2),2))))/
  48.             Sqrt(Power(Power(v1,2)*Power(Ybtm,2) -
  49.                Power(Ytop,2)*(Power(v2,2)*Power(Cos(phi2),2) +
  50.                   Power(v2,2)*Power(Sin(phi2),2)),2)),2) -
  51.         ((2*Power(Ybtm,2) - (2*Power(Ybtm,2)*
  52.                 (Power(v1,2)*Power(Ybtm,2) -
  53.                   Power(Ytop,2)*
  54.                    (Power(v2,2)*Power(Cos(phi2),2) +
  55.                      Power(v2,2)*Power(Sin(phi2),2))))/
  56.               Sqrt(Power(Power(v1,2)*Power(Ybtm,2) -
  57.                  Power(Ytop,2)*(Power(v2,2)*Power(Cos(phi2),2) +
  58.                     Power(v2,2)*Power(Sin(phi2),2)),2)))*
  59.            (Power(v1,2)*Power(Ybtm,2) +
  60.              Power(v2,2)*Power(Ytop,2)*Power(Cos(phi2),2) +
  61.              Power(v2,2)*Power(Ytop,2)*Power(Sin(phi2),2) -
  62.              Sqrt(Power(Power(v1,2)*Power(Ybtm,2) -
  63.                 Power(Ytop,2)*(Power(v2,2)*Power(Cos(phi2),2) +
  64.                    Power(v2,2)*Power(Sin(phi2),2)),2))))/2. -
  65.         (-1.5 + Log((Power(v1,2)*Power(Ybtm,2) +
  66.                Power(v2,2)*Power(Ytop,2)*Power(Cos(phi2),2) +
  67.                Power(v2,2)*Power(Ytop,2)*Power(Sin(phi2),2) -
  68.                Sqrt(Power(Power(v1,2)*Power(Ybtm,2) -
  69.                   Power(Ytop,2)*
  70.                    (Power(v2,2)*Power(Cos(phi2),2) +
  71.                      Power(v2,2)*Power(Sin(phi2),2)),2)))/(2.*Power(Rsc,2))))*
  72.          (2*Power(Ybtm,2) - (2*Power(Ybtm,2)*
  73.               (Power(v1,2)*Power(Ybtm,2) -
  74.                 Power(Ytop,2)*(Power(v2,2)*Power(Cos(phi2),2) +
  75.                    Power(v2,2)*Power(Sin(phi2),2))))/
  76.             Sqrt(Power(Power(v1,2)*Power(Ybtm,2) -
  77.                Power(Ytop,2)*(Power(v2,2)*Power(Cos(phi2),2) +
  78.                   Power(v2,2)*Power(Sin(phi2),2)),2)))*
  79.          (Power(v1,2)*Power(Ybtm,2) +
  80.            Power(v2,2)*Power(Ytop,2)*Power(Cos(phi2),2) +
  81.            Power(v2,2)*Power(Ytop,2)*Power(Sin(phi2),2) -
  82.            Sqrt(Power(Power(v1,2)*Power(Ybtm,2) -
  83.               Power(Ytop,2)*(Power(v2,2)*Power(Cos(phi2),2) +
  84.                  Power(v2,2)*Power(Sin(phi2),2)),2))) -
  85.         ((2*Power(Ybtm,2) + (2*Power(Ybtm,2)*
  86.                 (Power(v1,2)*Power(Ybtm,2) -
  87.                   Power(Ytop,2)*
  88.                    (Power(v2,2)*Power(Cos(phi2),2) +
  89.                      Power(v2,2)*Power(Sin(phi2),2))))/
  90.               Sqrt(Power(Power(v1,2)*Power(Ybtm,2) -
  91.                  Power(Ytop,2)*(Power(v2,2)*Power(Cos(phi2),2) +
  92.                     Power(v2,2)*Power(Sin(phi2),2)),2)))*
  93.            (Power(v1,2)*Power(Ybtm,2) +
  94.              Power(v2,2)*Power(Ytop,2)*Power(Cos(phi2),2) +
  95.              Power(v2,2)*Power(Ytop,2)*Power(Sin(phi2),2) +
  96.              Sqrt(Power(Power(v1,2)*Power(Ybtm,2) -
  97.                 Power(Ytop,2)*(Power(v2,2)*Power(Cos(phi2),2) +
  98.                    Power(v2,2)*Power(Sin(phi2),2)),2))))/2. -
  99.         (-1.5 + Log((Power(v1,2)*Power(Ybtm,2) +
  100.                Power(v2,2)*Power(Ytop,2)*Power(Cos(phi2),2) +
  101.                Power(v2,2)*Power(Ytop,2)*Power(Sin(phi2),2) +
  102.                Sqrt(Power(Power(v1,2)*Power(Ybtm,2) -
  103.                   Power(Ytop,2)*
  104.                    (Power(v2,2)*Power(Cos(phi2),2) +
  105.                      Power(v2,2)*Power(Sin(phi2),2)),2)))/(2.*Power(Rsc,2))))*
  106.          (2*Power(Ybtm,2) + (2*Power(Ybtm,2)*
  107.               (Power(v1,2)*Power(Ybtm,2) -
  108.                 Power(Ytop,2)*(Power(v2,2)*Power(Cos(phi2),2) +
  109.                    Power(v2,2)*Power(Sin(phi2),2))))/
  110.             Sqrt(Power(Power(v1,2)*Power(Ybtm,2) -
  111.                Power(Ytop,2)*(Power(v2,2)*Power(Cos(phi2),2) +
  112.                   Power(v2,2)*Power(Sin(phi2),2)),2)))*
  113.          (Power(v1,2)*Power(Ybtm,2) +
  114.            Power(v2,2)*Power(Ytop,2)*Power(Cos(phi2),2) +
  115.            Power(v2,2)*Power(Ytop,2)*Power(Sin(phi2),2) +
  116.            Sqrt(Power(Power(v1,2)*Power(Ybtm,2) -
  117.               Power(Ytop,2)*(Power(v2,2)*Power(Cos(phi2),2) +
  118.                  Power(v2,2)*Power(Sin(phi2),2)),2))) +
  119.         (3*Power(Complex(0,4)*lambda*Power(Ytop,2)*
  120.               (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))*
  121.               (Atop*(v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) -
  122.                 Complex(0,1)*lambda*v1*
  123.                  (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS))) -
  124.              Complex(0,4)*lambda*Power(Ytop,2)*
  125.               (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS))*
  126.               (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) +
  127.                 Complex(0,1)*lambda*v1*
  128.                  (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))),2))/
  129.          (8.*(Power(Power(MQ33,2) - Power(MU33,2),2) +
  130.              4*Power(Ytop,2)*(Atop*
  131.                  (v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) -
  132.                 Complex(0,1)*lambda*v1*
  133.                  (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS)))*
  134.               (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) +
  135.                 Complex(0,1)*lambda*v1*
  136.                  (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))))) +
  137.         ((-1.5 + Log((Power(MQ33,2) + Power(MU33,2) +
  138.                  2*Power(Ytop,2)*
  139.                   (Power(v2,2)*Power(Cos(phi2),2) +
  140.                     Power(v2,2)*Power(Sin(phi2),2)) -
  141.                  Sqrt(Power(Power(MQ33,2) - Power(MU33,2),2) +
  142.                    4*Power(Ytop,2)*
  143.                     (Atop*(v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) -
  144.                       Complex(0,1)*lambda*v1*
  145.                        (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS)))*
  146.                     (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) +
  147.                       Complex(0,1)*lambda*v1*
  148.                        (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS)))))/
  149.                (2.*Power(Rsc,2))))*
  150.            Power(Complex(0,4)*lambda*Power(Ytop,2)*
  151.               (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))*
  152.               (Atop*(v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) -
  153.                 Complex(0,1)*lambda*v1*
  154.                  (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS))) -
  155.              Complex(0,4)*lambda*Power(Ytop,2)*
  156.               (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS))*
  157.               (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) +
  158.                 Complex(0,1)*lambda*v1*
  159.                  (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))),2))/
  160.          (8.*(Power(Power(MQ33,2) - Power(MU33,2),2) +
  161.              4*Power(Ytop,2)*(Atop*
  162.                  (v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) -
  163.                 Complex(0,1)*lambda*v1*
  164.                  (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS)))*
  165.               (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) +
  166.                 Complex(0,1)*lambda*v1*
  167.                  (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))))) +
  168.         ((-1.5 + Log((Power(MQ33,2) + Power(MU33,2) +
  169.                  2*Power(Ytop,2)*
  170.                   (Power(v2,2)*Power(Cos(phi2),2) +
  171.                     Power(v2,2)*Power(Sin(phi2),2)) +
  172.                  Sqrt(Power(Power(MQ33,2) - Power(MU33,2),2) +
  173.                    4*Power(Ytop,2)*
  174.                     (Atop*(v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) -
  175.                       Complex(0,1)*lambda*v1*
  176.                        (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS)))*
  177.                     (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) +
  178.                       Complex(0,1)*lambda*v1*
  179.                        (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS)))))/
  180.                (2.*Power(Rsc,2))))*
  181.            Power(Complex(0,4)*lambda*Power(Ytop,2)*
  182.               (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))*
  183.               (Atop*(v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) -
  184.                 Complex(0,1)*lambda*v1*
  185.                  (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS))) -
  186.              Complex(0,4)*lambda*Power(Ytop,2)*
  187.               (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS))*
  188.               (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) +
  189.                 Complex(0,1)*lambda*v1*
  190.                  (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))),2))/
  191.          (8.*(Power(Power(MQ33,2) - Power(MU33,2),2) +
  192.              4*Power(Ytop,2)*(Atop*
  193.                  (v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) -
  194.                 Complex(0,1)*lambda*v1*
  195.                  (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS)))*
  196.               (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) +
  197.                 Complex(0,1)*lambda*v1*
  198.                  (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))))) +
  199.         (Power(Complex(0,4)*lambda*Power(Ytop,2)*
  200.               (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))*
  201.               (Atop*(v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) -
  202.                 Complex(0,1)*lambda*v1*
  203.                  (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS))) -
  204.              Complex(0,4)*lambda*Power(Ytop,2)*
  205.               (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS))*
  206.               (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) +
  207.                 Complex(0,1)*lambda*v1*
  208.                  (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))),2)*
  209.            (Power(MQ33,2) + Power(MU33,2) +
  210.              2*Power(Ytop,2)*(Power(v2,2)*Power(Cos(phi2),2) +
  211.                 Power(v2,2)*Power(Sin(phi2),2)) -
  212.              Sqrt(Power(Power(MQ33,2) - Power(MU33,2),2) +
  213.                4*Power(Ytop,2)*(Atop*
  214.                    (v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) -
  215.                   Complex(0,1)*lambda*v1*
  216.                    (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS)))*
  217.                 (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) +
  218.                   Complex(0,1)*lambda*v1*
  219.                    (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))))))/
  220.          (16.*Power(Power(Power(MQ33,2) - Power(MU33,2),2) +
  221.              4*Power(Ytop,2)*(Atop*
  222.                  (v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) -
  223.                 Complex(0,1)*lambda*v1*
  224.                  (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS)))*
  225.               (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) +
  226.                 Complex(0,1)*lambda*v1*
  227.                  (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))),1.5)) +
  228.         ((-1.5 + Log((Power(MQ33,2) + Power(MU33,2) +
  229.                  2*Power(Ytop,2)*
  230.                   (Power(v2,2)*Power(Cos(phi2),2) +
  231.                     Power(v2,2)*Power(Sin(phi2),2)) -
  232.                  Sqrt(Power(Power(MQ33,2) - Power(MU33,2),2) +
  233.                    4*Power(Ytop,2)*
  234.                     (Atop*(v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) -
  235.                       Complex(0,1)*lambda*v1*
  236.                        (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS)))*
  237.                     (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) +
  238.                       Complex(0,1)*lambda*v1*
  239.                        (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS)))))/
  240.                (2.*Power(Rsc,2))))*
  241.            Power(Complex(0,4)*lambda*Power(Ytop,2)*
  242.               (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))*
  243.               (Atop*(v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) -
  244.                 Complex(0,1)*lambda*v1*
  245.                  (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS))) -
  246.              Complex(0,4)*lambda*Power(Ytop,2)*
  247.               (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS))*
  248.               (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) +
  249.                 Complex(0,1)*lambda*v1*
  250.                  (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))),2)*
  251.            (Power(MQ33,2) + Power(MU33,2) +
  252.              2*Power(Ytop,2)*(Power(v2,2)*Power(Cos(phi2),2) +
  253.                 Power(v2,2)*Power(Sin(phi2),2)) -
  254.              Sqrt(Power(Power(MQ33,2) - Power(MU33,2),2) +
  255.                4*Power(Ytop,2)*(Atop*
  256.                    (v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) -
  257.                   Complex(0,1)*lambda*v1*
  258.                    (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS)))*
  259.                 (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) +
  260.                   Complex(0,1)*lambda*v1*
  261.                    (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))))))/
  262.          (8.*Power(Power(Power(MQ33,2) - Power(MU33,2),2) +
  263.              4*Power(Ytop,2)*(Atop*
  264.                  (v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) -
  265.                 Complex(0,1)*lambda*v1*
  266.                  (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS)))*
  267.               (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) +
  268.                 Complex(0,1)*lambda*v1*
  269.                  (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))),1.5)) -
  270.         (Power(lambda,2)*Power(Ytop,2)*
  271.            (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS))*
  272.            (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))*
  273.            (Power(MQ33,2) + Power(MU33,2) +
  274.              2*Power(Ytop,2)*(Power(v2,2)*Power(Cos(phi2),2) +
  275.                 Power(v2,2)*Power(Sin(phi2),2)) -
  276.              Sqrt(Power(Power(MQ33,2) - Power(MU33,2),2) +
  277.                4*Power(Ytop,2)*(Atop*
  278.                    (v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) -
  279.                   Complex(0,1)*lambda*v1*
  280.                    (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS)))*
  281.                 (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) +
  282.                   Complex(0,1)*lambda*v1*
  283.                    (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))))))/
  284.          Sqrt(Power(Power(MQ33,2) - Power(MU33,2),2) +
  285.            4*Power(Ytop,2)*(Atop*(v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) -
  286.               Complex(0,1)*lambda*v1*
  287.                (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS)))*
  288.             (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) +
  289.               Complex(0,1)*lambda*v1*
  290.                (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS)))) -
  291.         (2*Power(lambda,2)*Power(Ytop,2)*
  292.            (-1.5 + Log((Power(MQ33,2) + Power(MU33,2) +
  293.                  2*Power(Ytop,2)*
  294.                   (Power(v2,2)*Power(Cos(phi2),2) +
  295.                     Power(v2,2)*Power(Sin(phi2),2)) -
  296.                  Sqrt(Power(Power(MQ33,2) - Power(MU33,2),2) +
  297.                    4*Power(Ytop,2)*
  298.                     (Atop*(v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) -
  299.                       Complex(0,1)*lambda*v1*
  300.                        (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS)))*
  301.                     (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) +
  302.                       Complex(0,1)*lambda*v1*
  303.                        (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS)))))/
  304.                (2.*Power(Rsc,2))))*
  305.            (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS))*
  306.            (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))*
  307.            (Power(MQ33,2) + Power(MU33,2) +
  308.              2*Power(Ytop,2)*(Power(v2,2)*Power(Cos(phi2),2) +
  309.                 Power(v2,2)*Power(Sin(phi2),2)) -
  310.              Sqrt(Power(Power(MQ33,2) - Power(MU33,2),2) +
  311.                4*Power(Ytop,2)*(Atop*
  312.                    (v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) -
  313.                   Complex(0,1)*lambda*v1*
  314.                    (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS)))*
  315.                 (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) +
  316.                   Complex(0,1)*lambda*v1*
  317.                    (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))))))/
  318.          Sqrt(Power(Power(MQ33,2) - Power(MU33,2),2) +
  319.            4*Power(Ytop,2)*(Atop*(v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) -
  320.               Complex(0,1)*lambda*v1*
  321.                (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS)))*
  322.             (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) +
  323.               Complex(0,1)*lambda*v1*
  324.                (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS)))) -
  325.         (Power(Complex(0,4)*lambda*Power(Ytop,2)*
  326.               (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))*
  327.               (Atop*(v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) -
  328.                 Complex(0,1)*lambda*v1*
  329.                  (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS))) -
  330.              Complex(0,4)*lambda*Power(Ytop,2)*
  331.               (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS))*
  332.               (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) +
  333.                 Complex(0,1)*lambda*v1*
  334.                  (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))),2)*
  335.            (Power(MQ33,2) + Power(MU33,2) +
  336.              2*Power(Ytop,2)*(Power(v2,2)*Power(Cos(phi2),2) +
  337.                 Power(v2,2)*Power(Sin(phi2),2)) +
  338.              Sqrt(Power(Power(MQ33,2) - Power(MU33,2),2) +
  339.                4*Power(Ytop,2)*(Atop*
  340.                    (v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) -
  341.                   Complex(0,1)*lambda*v1*
  342.                    (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS)))*
  343.                 (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) +
  344.                   Complex(0,1)*lambda*v1*
  345.                    (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))))))/
  346.          (16.*Power(Power(Power(MQ33,2) - Power(MU33,2),2) +
  347.              4*Power(Ytop,2)*(Atop*
  348.                  (v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) -
  349.                 Complex(0,1)*lambda*v1*
  350.                  (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS)))*
  351.               (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) +
  352.                 Complex(0,1)*lambda*v1*
  353.                  (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))),1.5)) -
  354.         ((-1.5 + Log((Power(MQ33,2) + Power(MU33,2) +
  355.                  2*Power(Ytop,2)*
  356.                   (Power(v2,2)*Power(Cos(phi2),2) +
  357.                     Power(v2,2)*Power(Sin(phi2),2)) +
  358.                  Sqrt(Power(Power(MQ33,2) - Power(MU33,2),2) +
  359.                    4*Power(Ytop,2)*
  360.                     (Atop*(v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) -
  361.                       Complex(0,1)*lambda*v1*
  362.                        (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS)))*
  363.                     (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) +
  364.                       Complex(0,1)*lambda*v1*
  365.                        (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS)))))/
  366.                (2.*Power(Rsc,2))))*
  367.            Power(Complex(0,4)*lambda*Power(Ytop,2)*
  368.               (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))*
  369.               (Atop*(v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) -
  370.                 Complex(0,1)*lambda*v1*
  371.                  (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS))) -
  372.              Complex(0,4)*lambda*Power(Ytop,2)*
  373.               (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS))*
  374.               (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) +
  375.                 Complex(0,1)*lambda*v1*
  376.                  (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))),2)*
  377.            (Power(MQ33,2) + Power(MU33,2) +
  378.              2*Power(Ytop,2)*(Power(v2,2)*Power(Cos(phi2),2) +
  379.                 Power(v2,2)*Power(Sin(phi2),2)) +
  380.              Sqrt(Power(Power(MQ33,2) - Power(MU33,2),2) +
  381.                4*Power(Ytop,2)*(Atop*
  382.                    (v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) -
  383.                   Complex(0,1)*lambda*v1*
  384.                    (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS)))*
  385.                 (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) +
  386.                   Complex(0,1)*lambda*v1*
  387.                    (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))))))/
  388.          (8.*Power(Power(Power(MQ33,2) - Power(MU33,2),2) +
  389.              4*Power(Ytop,2)*(Atop*
  390.                  (v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) -
  391.                 Complex(0,1)*lambda*v1*
  392.                  (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS)))*
  393.               (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) +
  394.                 Complex(0,1)*lambda*v1*
  395.                  (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))),1.5)) +
  396.         (Power(lambda,2)*Power(Ytop,2)*
  397.            (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS))*
  398.            (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))*
  399.            (Power(MQ33,2) + Power(MU33,2) +
  400.              2*Power(Ytop,2)*(Power(v2,2)*Power(Cos(phi2),2) +
  401.                 Power(v2,2)*Power(Sin(phi2),2)) +
  402.              Sqrt(Power(Power(MQ33,2) - Power(MU33,2),2) +
  403.                4*Power(Ytop,2)*(Atop*
  404.                    (v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) -
  405.                   Complex(0,1)*lambda*v1*
  406.                    (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS)))*
  407.                 (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) +
  408.                   Complex(0,1)*lambda*v1*
  409.                    (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))))))/
  410.          Sqrt(Power(Power(MQ33,2) - Power(MU33,2),2) +
  411.            4*Power(Ytop,2)*(Atop*(v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) -
  412.               Complex(0,1)*lambda*v1*
  413.                (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS)))*
  414.             (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) +
  415.               Complex(0,1)*lambda*v1*
  416.                (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS)))) +
  417.         (2*Power(lambda,2)*Power(Ytop,2)*
  418.            (-1.5 + Log((Power(MQ33,2) + Power(MU33,2) +
  419.                  2*Power(Ytop,2)*
  420.                   (Power(v2,2)*Power(Cos(phi2),2) +
  421.                     Power(v2,2)*Power(Sin(phi2),2)) +
  422.                  Sqrt(Power(Power(MQ33,2) - Power(MU33,2),2) +
  423.                    4*Power(Ytop,2)*
  424.                     (Atop*(v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) -
  425.                       Complex(0,1)*lambda*v1*
  426.                        (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS)))*
  427.                     (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) +
  428.                       Complex(0,1)*lambda*v1*
  429.                        (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS)))))/
  430.                (2.*Power(Rsc,2))))*
  431.            (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS))*
  432.            (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))*
  433.            (Power(MQ33,2) + Power(MU33,2) +
  434.              2*Power(Ytop,2)*(Power(v2,2)*Power(Cos(phi2),2) +
  435.                 Power(v2,2)*Power(Sin(phi2),2)) +
  436.              Sqrt(Power(Power(MQ33,2) - Power(MU33,2),2) +
  437.                4*Power(Ytop,2)*(Atop*
  438.                    (v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) -
  439.                   Complex(0,1)*lambda*v1*
  440.                    (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS)))*
  441.                 (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) +
  442.                   Complex(0,1)*lambda*v1*
  443.                    (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))))))/
  444.          Sqrt(Power(Power(MQ33,2) - Power(MU33,2),2) +
  445.            4*Power(Ytop,2)*(Atop*(v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) -
  446.               Complex(0,1)*lambda*v1*
  447.                (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS)))*
  448.             (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) +
  449.               Complex(0,1)*lambda*v1*
  450.                (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS)))) +
  451.         (3*Power(4*v1*Power(Ybtm,2) -
  452.              (2*Power(Ybtm,2)*(2*Power(Abtm,2)*v1 +
  453.                   2*Abtm*lambda*
  454.                    (-(v2*vevS*Cos(phi2)*Cos(phiS)) +
  455.                      v2*vevS*Sin(phi2)*Sin(phiS))))/
  456.               Sqrt(Power(MD33,4) - 2*Power(MD33,2)*Power(MQ33,2) +
  457.                 Power(MQ33,4) +
  458.                 4*Power(Ybtm,2)*
  459.                  (Power(Abtm,2)*Power(v1,2) +
  460.                    2*Abtm*lambda*
  461.                     (-(v1*v2*vevS*Cos(phi2)*Cos(phiS)) +
  462.                       v1*v2*vevS*Sin(phi2)*Sin(phiS)) +
  463.                    Power(lambda,2)*
  464.                     (Power(v2,2)*Power(Cos(phi2),2) +
  465.                       Power(v2,2)*Power(Sin(phi2),2))*
  466.                     (Power(vevS,2)*Power(Cos(phiS),2) +
  467.                       Power(vevS,2)*Power(Sin(phiS),2)))),2))/4. +
  468.         ((-1.5 + Log((Power(MD33,2) + Power(MQ33,2) +
  469.                  2*Power(v1,2)*Power(Ybtm,2) -
  470.                  Sqrt(Power(MD33,4) - 2*Power(MD33,2)*Power(MQ33,2) +
  471.                    Power(MQ33,4) +
  472.                    4*Power(Ybtm,2)*
  473.                     (Power(Abtm,2)*Power(v1,2) +
  474.                       2*Abtm*lambda*
  475.                        (-(v1*v2*vevS*Cos(phi2)*Cos(phiS)) +
  476.                          v1*v2*vevS*Sin(phi2)*Sin(phiS)) +
  477.                       Power(lambda,2)*
  478.                        (Power(v2,2)*Power(Cos(phi2),2) +
  479.                          Power(v2,2)*Power(Sin(phi2),2))*
  480.                        (Power(vevS,2)*Power(Cos(phiS),2) +
  481.                          Power(vevS,2)*Power(Sin(phiS),2)))))/
  482.                (2.*Power(Rsc,2))))*
  483.            Power(4*v1*Power(Ybtm,2) -
  484.              (2*Power(Ybtm,2)*(2*Power(Abtm,2)*v1 +
  485.                   2*Abtm*lambda*
  486.                    (-(v2*vevS*Cos(phi2)*Cos(phiS)) +
  487.                      v2*vevS*Sin(phi2)*Sin(phiS))))/
  488.               Sqrt(Power(MD33,4) - 2*Power(MD33,2)*Power(MQ33,2) +
  489.                 Power(MQ33,4) +
  490.                 4*Power(Ybtm,2)*
  491.                  (Power(Abtm,2)*Power(v1,2) +
  492.                    2*Abtm*lambda*
  493.                     (-(v1*v2*vevS*Cos(phi2)*Cos(phiS)) +
  494.                       v1*v2*vevS*Sin(phi2)*Sin(phiS)) +
  495.                    Power(lambda,2)*
  496.                     (Power(v2,2)*Power(Cos(phi2),2) +
  497.                       Power(v2,2)*Power(Sin(phi2),2))*
  498.                     (Power(vevS,2)*Power(Cos(phiS),2) +
  499.                       Power(vevS,2)*Power(Sin(phiS),2)))),2))/2. +
  500.         (3*Power(4*v1*Power(Ybtm,2) +
  501.              (2*Power(Ybtm,2)*(2*Power(Abtm,2)*v1 +
  502.                   2*Abtm*lambda*
  503.                    (-(v2*vevS*Cos(phi2)*Cos(phiS)) +
  504.                      v2*vevS*Sin(phi2)*Sin(phiS))))/
  505.               Sqrt(Power(MD33,4) - 2*Power(MD33,2)*Power(MQ33,2) +
  506.                 Power(MQ33,4) +
  507.                 4*Power(Ybtm,2)*
  508.                  (Power(Abtm,2)*Power(v1,2) +
  509.                    2*Abtm*lambda*
  510.                     (-(v1*v2*vevS*Cos(phi2)*Cos(phiS)) +
  511.                       v1*v2*vevS*Sin(phi2)*Sin(phiS)) +
  512.                    Power(lambda,2)*
  513.                     (Power(v2,2)*Power(Cos(phi2),2) +
  514.                       Power(v2,2)*Power(Sin(phi2),2))*
  515.                     (Power(vevS,2)*Power(Cos(phiS),2) +
  516.                       Power(vevS,2)*Power(Sin(phiS),2)))),2))/4. +
  517.         ((-1.5 + Log((Power(MD33,2) + Power(MQ33,2) +
  518.                  2*Power(v1,2)*Power(Ybtm,2) +
  519.                  Sqrt(Power(MD33,4) - 2*Power(MD33,2)*Power(MQ33,2) +
  520.                    Power(MQ33,4) +
  521.                    4*Power(Ybtm,2)*
  522.                     (Power(Abtm,2)*Power(v1,2) +
  523.                       2*Abtm*lambda*
  524.                        (-(v1*v2*vevS*Cos(phi2)*Cos(phiS)) +
  525.                          v1*v2*vevS*Sin(phi2)*Sin(phiS)) +
  526.                       Power(lambda,2)*
  527.                        (Power(v2,2)*Power(Cos(phi2),2) +
  528.                          Power(v2,2)*Power(Sin(phi2),2))*
  529.                        (Power(vevS,2)*Power(Cos(phiS),2) +
  530.                          Power(vevS,2)*Power(Sin(phiS),2)))))/
  531.                (2.*Power(Rsc,2))))*
  532.            Power(4*v1*Power(Ybtm,2) +
  533.              (2*Power(Ybtm,2)*(2*Power(Abtm,2)*v1 +
  534.                   2*Abtm*lambda*
  535.                    (-(v2*vevS*Cos(phi2)*Cos(phiS)) +
  536.                      v2*vevS*Sin(phi2)*Sin(phiS))))/
  537.               Sqrt(Power(MD33,4) - 2*Power(MD33,2)*Power(MQ33,2) +
  538.                 Power(MQ33,4) +
  539.                 4*Power(Ybtm,2)*
  540.                  (Power(Abtm,2)*Power(v1,2) +
  541.                    2*Abtm*lambda*
  542.                     (-(v1*v2*vevS*Cos(phi2)*Cos(phiS)) +
  543.                       v1*v2*vevS*Sin(phi2)*Sin(phiS)) +
  544.                    Power(lambda,2)*
  545.                     (Power(v2,2)*Power(Cos(phi2),2) +
  546.                       Power(v2,2)*Power(Sin(phi2),2))*
  547.                     (Power(vevS,2)*Power(Cos(phiS),2) +
  548.                       Power(vevS,2)*Power(Sin(phiS),2)))),2))/2. +
  549.         ((4*Power(Ybtm,2) + (4*Power(Ybtm,4)*
  550.                 Power(2*Power(Abtm,2)*v1 +
  551.                   2*Abtm*lambda*
  552.                    (-(v2*vevS*Cos(phi2)*Cos(phiS)) +
  553.                      v2*vevS*Sin(phi2)*Sin(phiS)),2))/
  554.               Power(Power(MD33,4) - 2*Power(MD33,2)*Power(MQ33,2) +
  555.                 Power(MQ33,4) +
  556.                 4*Power(Ybtm,2)*
  557.                  (Power(Abtm,2)*Power(v1,2) +
  558.                    2*Abtm*lambda*
  559.                     (-(v1*v2*vevS*Cos(phi2)*Cos(phiS)) +
  560.                       v1*v2*vevS*Sin(phi2)*Sin(phiS)) +
  561.                    Power(lambda,2)*
  562.                     (Power(v2,2)*Power(Cos(phi2),2) +
  563.                       Power(v2,2)*Power(Sin(phi2),2))*
  564.                     (Power(vevS,2)*Power(Cos(phiS),2) +
  565.                       Power(vevS,2)*Power(Sin(phiS),2))),1.5) -
  566.              (4*Power(Abtm,2)*Power(Ybtm,2))/
  567.               Sqrt(Power(MD33,4) - 2*Power(MD33,2)*Power(MQ33,2) +
  568.                 Power(MQ33,4) +
  569.                 4*Power(Ybtm,2)*
  570.                  (Power(Abtm,2)*Power(v1,2) +
  571.                    2*Abtm*lambda*
  572.                     (-(v1*v2*vevS*Cos(phi2)*Cos(phiS)) +
  573.                       v1*v2*vevS*Sin(phi2)*Sin(phiS)) +
  574.                    Power(lambda,2)*
  575.                     (Power(v2,2)*Power(Cos(phi2),2) +
  576.                       Power(v2,2)*Power(Sin(phi2),2))*
  577.                     (Power(vevS,2)*Power(Cos(phiS),2) +
  578.                       Power(vevS,2)*Power(Sin(phiS),2)))))*
  579.            (Power(MD33,2) + Power(MQ33,2) + 2*Power(v1,2)*Power(Ybtm,2) -
  580.              Sqrt(Power(MD33,4) - 2*Power(MD33,2)*Power(MQ33,2) +
  581.                Power(MQ33,4) + 4*Power(Ybtm,2)*
  582.                 (Power(Abtm,2)*Power(v1,2) +
  583.                   2*Abtm*lambda*
  584.                    (-(v1*v2*vevS*Cos(phi2)*Cos(phiS)) +
  585.                      v1*v2*vevS*Sin(phi2)*Sin(phiS)) +
  586.                   Power(lambda,2)*
  587.                    (Power(v2,2)*Power(Cos(phi2),2) +
  588.                      Power(v2,2)*Power(Sin(phi2),2))*
  589.                    (Power(vevS,2)*Power(Cos(phiS),2) +
  590.                      Power(vevS,2)*Power(Sin(phiS),2))))))/4. +
  591.         ((-1.5 + Log((Power(MD33,2) + Power(MQ33,2) +
  592.                  2*Power(v1,2)*Power(Ybtm,2) -
  593.                  Sqrt(Power(MD33,4) - 2*Power(MD33,2)*Power(MQ33,2) +
  594.                    Power(MQ33,4) +
  595.                    4*Power(Ybtm,2)*
  596.                     (Power(Abtm,2)*Power(v1,2) +
  597.                       2*Abtm*lambda*
  598.                        (-(v1*v2*vevS*Cos(phi2)*Cos(phiS)) +
  599.                          v1*v2*vevS*Sin(phi2)*Sin(phiS)) +
  600.                       Power(lambda,2)*
  601.                        (Power(v2,2)*Power(Cos(phi2),2) +
  602.                          Power(v2,2)*Power(Sin(phi2),2))*
  603.                        (Power(vevS,2)*Power(Cos(phiS),2) +
  604.                          Power(vevS,2)*Power(Sin(phiS),2)))))/
  605.                (2.*Power(Rsc,2))))*
  606.            (4*Power(Ybtm,2) + (4*Power(Ybtm,4)*
  607.                 Power(2*Power(Abtm,2)*v1 +
  608.                   2*Abtm*lambda*
  609.                    (-(v2*vevS*Cos(phi2)*Cos(phiS)) +
  610.                      v2*vevS*Sin(phi2)*Sin(phiS)),2))/
  611.               Power(Power(MD33,4) - 2*Power(MD33,2)*Power(MQ33,2) +
  612.                 Power(MQ33,4) +
  613.                 4*Power(Ybtm,2)*
  614.                  (Power(Abtm,2)*Power(v1,2) +
  615.                    2*Abtm*lambda*
  616.                     (-(v1*v2*vevS*Cos(phi2)*Cos(phiS)) +
  617.                       v1*v2*vevS*Sin(phi2)*Sin(phiS)) +
  618.                    Power(lambda,2)*
  619.                     (Power(v2,2)*Power(Cos(phi2),2) +
  620.                       Power(v2,2)*Power(Sin(phi2),2))*
  621.                     (Power(vevS,2)*Power(Cos(phiS),2) +
  622.                       Power(vevS,2)*Power(Sin(phiS),2))),1.5) -
  623.              (4*Power(Abtm,2)*Power(Ybtm,2))/
  624.               Sqrt(Power(MD33,4) - 2*Power(MD33,2)*Power(MQ33,2) +
  625.                 Power(MQ33,4) +
  626.                 4*Power(Ybtm,2)*
  627.                  (Power(Abtm,2)*Power(v1,2) +
  628.                    2*Abtm*lambda*
  629.                     (-(v1*v2*vevS*Cos(phi2)*Cos(phiS)) +
  630.                       v1*v2*vevS*Sin(phi2)*Sin(phiS)) +
  631.                    Power(lambda,2)*
  632.                     (Power(v2,2)*Power(Cos(phi2),2) +
  633.                       Power(v2,2)*Power(Sin(phi2),2))*
  634.                     (Power(vevS,2)*Power(Cos(phiS),2) +
  635.                       Power(vevS,2)*Power(Sin(phiS),2)))))*
  636.            (Power(MD33,2) + Power(MQ33,2) + 2*Power(v1,2)*Power(Ybtm,2) -
  637.              Sqrt(Power(MD33,4) - 2*Power(MD33,2)*Power(MQ33,2) +
  638.                Power(MQ33,4) + 4*Power(Ybtm,2)*
  639.                 (Power(Abtm,2)*Power(v1,2) +
  640.                   2*Abtm*lambda*
  641.                    (-(v1*v2*vevS*Cos(phi2)*Cos(phiS)) +
  642.                      v1*v2*vevS*Sin(phi2)*Sin(phiS)) +
  643.                   Power(lambda,2)*
  644.                    (Power(v2,2)*Power(Cos(phi2),2) +
  645.                      Power(v2,2)*Power(Sin(phi2),2))*
  646.                    (Power(vevS,2)*Power(Cos(phiS),2) +
  647.                      Power(vevS,2)*Power(Sin(phiS),2))))))/2. +
  648.         ((4*Power(Ybtm,2) - (4*Power(Ybtm,4)*
  649.                 Power(2*Power(Abtm,2)*v1 +
  650.                   2*Abtm*lambda*
  651.                    (-(v2*vevS*Cos(phi2)*Cos(phiS)) +
  652.                      v2*vevS*Sin(phi2)*Sin(phiS)),2))/
  653.               Power(Power(MD33,4) - 2*Power(MD33,2)*Power(MQ33,2) +
  654.                 Power(MQ33,4) +
  655.                 4*Power(Ybtm,2)*
  656.                  (Power(Abtm,2)*Power(v1,2) +
  657.                    2*Abtm*lambda*
  658.                     (-(v1*v2*vevS*Cos(phi2)*Cos(phiS)) +
  659.                       v1*v2*vevS*Sin(phi2)*Sin(phiS)) +
  660.                    Power(lambda,2)*
  661.                     (Power(v2,2)*Power(Cos(phi2),2) +
  662.                       Power(v2,2)*Power(Sin(phi2),2))*
  663.                     (Power(vevS,2)*Power(Cos(phiS),2) +
  664.                       Power(vevS,2)*Power(Sin(phiS),2))),1.5) +
  665.              (4*Power(Abtm,2)*Power(Ybtm,2))/
  666.               Sqrt(Power(MD33,4) - 2*Power(MD33,2)*Power(MQ33,2) +
  667.                 Power(MQ33,4) +
  668.                 4*Power(Ybtm,2)*
  669.                  (Power(Abtm,2)*Power(v1,2) +
  670.                    2*Abtm*lambda*
  671.                     (-(v1*v2*vevS*Cos(phi2)*Cos(phiS)) +
  672.                       v1*v2*vevS*Sin(phi2)*Sin(phiS)) +
  673.                    Power(lambda,2)*
  674.                     (Power(v2,2)*Power(Cos(phi2),2) +
  675.                       Power(v2,2)*Power(Sin(phi2),2))*
  676.                     (Power(vevS,2)*Power(Cos(phiS),2) +
  677.                       Power(vevS,2)*Power(Sin(phiS),2)))))*
  678.            (Power(MD33,2) + Power(MQ33,2) + 2*Power(v1,2)*Power(Ybtm,2) +
  679.              Sqrt(Power(MD33,4) - 2*Power(MD33,2)*Power(MQ33,2) +
  680.                Power(MQ33,4) + 4*Power(Ybtm,2)*
  681.                 (Power(Abtm,2)*Power(v1,2) +
  682.                   2*Abtm*lambda*
  683.                    (-(v1*v2*vevS*Cos(phi2)*Cos(phiS)) +
  684.                      v1*v2*vevS*Sin(phi2)*Sin(phiS)) +
  685.                   Power(lambda,2)*
  686.                    (Power(v2,2)*Power(Cos(phi2),2) +
  687.                      Power(v2,2)*Power(Sin(phi2),2))*
  688.                    (Power(vevS,2)*Power(Cos(phiS),2) +
  689.                      Power(vevS,2)*Power(Sin(phiS),2))))))/4. +
  690.         ((-1.5 + Log((Power(MD33,2) + Power(MQ33,2) +
  691.                  2*Power(v1,2)*Power(Ybtm,2) +
  692.                  Sqrt(Power(MD33,4) - 2*Power(MD33,2)*Power(MQ33,2) +
  693.                    Power(MQ33,4) +
  694.                    4*Power(Ybtm,2)*
  695.                     (Power(Abtm,2)*Power(v1,2) +
  696.                       2*Abtm*lambda*
  697.                        (-(v1*v2*vevS*Cos(phi2)*Cos(phiS)) +
  698.                          v1*v2*vevS*Sin(phi2)*Sin(phiS)) +
  699.                       Power(lambda,2)*
  700.                        (Power(v2,2)*Power(Cos(phi2),2) +
  701.                          Power(v2,2)*Power(Sin(phi2),2))*
  702.                        (Power(vevS,2)*Power(Cos(phiS),2) +
  703.                          Power(vevS,2)*Power(Sin(phiS),2)))))/
  704.                (2.*Power(Rsc,2))))*
  705.            (4*Power(Ybtm,2) - (4*Power(Ybtm,4)*
  706.                 Power(2*Power(Abtm,2)*v1 +
  707.                   2*Abtm*lambda*
  708.                    (-(v2*vevS*Cos(phi2)*Cos(phiS)) +
  709.                      v2*vevS*Sin(phi2)*Sin(phiS)),2))/
  710.               Power(Power(MD33,4) - 2*Power(MD33,2)*Power(MQ33,2) +
  711.                 Power(MQ33,4) +
  712.                 4*Power(Ybtm,2)*
  713.                  (Power(Abtm,2)*Power(v1,2) +
  714.                    2*Abtm*lambda*
  715.                     (-(v1*v2*vevS*Cos(phi2)*Cos(phiS)) +
  716.                       v1*v2*vevS*Sin(phi2)*Sin(phiS)) +
  717.                    Power(lambda,2)*
  718.                     (Power(v2,2)*Power(Cos(phi2),2) +
  719.                       Power(v2,2)*Power(Sin(phi2),2))*
  720.                     (Power(vevS,2)*Power(Cos(phiS),2) +
  721.                       Power(vevS,2)*Power(Sin(phiS),2))),1.5) +
  722.              (4*Power(Abtm,2)*Power(Ybtm,2))/
  723.               Sqrt(Power(MD33,4) - 2*Power(MD33,2)*Power(MQ33,2) +
  724.                 Power(MQ33,4) +
  725.                 4*Power(Ybtm,2)*
  726.                  (Power(Abtm,2)*Power(v1,2) +
  727.                    2*Abtm*lambda*
  728.                     (-(v1*v2*vevS*Cos(phi2)*Cos(phiS)) +
  729.                       v1*v2*vevS*Sin(phi2)*Sin(phiS)) +
  730.                    Power(lambda,2)*
  731.                     (Power(v2,2)*Power(Cos(phi2),2) +
  732.                       Power(v2,2)*Power(Sin(phi2),2))*
  733.                     (Power(vevS,2)*Power(Cos(phiS),2) +
  734.                       Power(vevS,2)*Power(Sin(phiS),2)))))*
  735.            (Power(MD33,2) + Power(MQ33,2) + 2*Power(v1,2)*Power(Ybtm,2) +
  736.              Sqrt(Power(MD33,4) - 2*Power(MD33,2)*Power(MQ33,2) +
  737.                Power(MQ33,4) + 4*Power(Ybtm,2)*
  738.                 (Power(Abtm,2)*Power(v1,2) +
  739.                   2*Abtm*lambda*
  740.                    (-(v1*v2*vevS*Cos(phi2)*Cos(phiS)) +
  741.                      v1*v2*vevS*Sin(phi2)*Sin(phiS)) +
  742.                   Power(lambda,2)*
  743.                    (Power(v2,2)*Power(Cos(phi2),2) +
  744.                      Power(v2,2)*Power(Sin(phi2),2))*
  745.                    (Power(vevS,2)*Power(Cos(phiS),2) +
  746.                      Power(vevS,2)*Power(Sin(phiS),2))))))/2.))/
  747.     (64.*Power(Pi,2))