MX11 = (4*Power(mH1,2) + 3*Power(g2,2)*Power(v1,2) - Power(g2,2)*Power(v2,2) + 
      4*Power(lambda,2)*Power(v2,2) + 
      Power(g1,2)*(3*Power(v1,2) - Power(v2,2)) + 
      4*Power(lambda,2)*Power(vevS,2))/4. + 
   (3*((-3*Power(2*v1*Power(Ybtm,2) - 
             (2*v1*Power(Ybtm,2)*
                (Power(v1,2)*Power(Ybtm,2) - 
                  Power(Ytop,2)*
                   (Power(v2,2)*Power(Cos(phi2),2) + 
                     Power(v2,2)*Power(Sin(phi2),2))))/
              Sqrt(Power(Power(v1,2)*Power(Ybtm,2) - 
                 Power(Ytop,2)*(Power(v2,2)*Power(Cos(phi2),2) + 
                    Power(v2,2)*Power(Sin(phi2),2)),2)),2))/2. - 
        (-1.5 + Log((Power(v1,2)*Power(Ybtm,2) + 
               Power(v2,2)*Power(Ytop,2)*Power(Cos(phi2),2) + 
               Power(v2,2)*Power(Ytop,2)*Power(Sin(phi2),2) - 
               Sqrt(Power(Power(v1,2)*Power(Ybtm,2) - 
                  Power(Ytop,2)*
                   (Power(v2,2)*Power(Cos(phi2),2) + 
                     Power(v2,2)*Power(Sin(phi2),2)),2)))/(2.*Power(Rsc,2))))*
         Power(2*v1*Power(Ybtm,2) - 
           (2*v1*Power(Ybtm,2)*(Power(v1,2)*Power(Ybtm,2) - 
                Power(Ytop,2)*(Power(v2,2)*Power(Cos(phi2),2) + 
                   Power(v2,2)*Power(Sin(phi2),2))))/
            Sqrt(Power(Power(v1,2)*Power(Ybtm,2) - 
               Power(Ytop,2)*(Power(v2,2)*Power(Cos(phi2),2) + 
                  Power(v2,2)*Power(Sin(phi2),2)),2)),2) - 
        (3*Power(2*v1*Power(Ybtm,2) + 
             (2*v1*Power(Ybtm,2)*
                (Power(v1,2)*Power(Ybtm,2) - 
                  Power(Ytop,2)*
                   (Power(v2,2)*Power(Cos(phi2),2) + 
                     Power(v2,2)*Power(Sin(phi2),2))))/
              Sqrt(Power(Power(v1,2)*Power(Ybtm,2) - 
                 Power(Ytop,2)*(Power(v2,2)*Power(Cos(phi2),2) + 
                    Power(v2,2)*Power(Sin(phi2),2)),2)),2))/2. - 
        (-1.5 + Log((Power(v1,2)*Power(Ybtm,2) + 
               Power(v2,2)*Power(Ytop,2)*Power(Cos(phi2),2) + 
               Power(v2,2)*Power(Ytop,2)*Power(Sin(phi2),2) + 
               Sqrt(Power(Power(v1,2)*Power(Ybtm,2) - 
                  Power(Ytop,2)*
                   (Power(v2,2)*Power(Cos(phi2),2) + 
                     Power(v2,2)*Power(Sin(phi2),2)),2)))/(2.*Power(Rsc,2))))*
         Power(2*v1*Power(Ybtm,2) + 
           (2*v1*Power(Ybtm,2)*(Power(v1,2)*Power(Ybtm,2) - 
                Power(Ytop,2)*(Power(v2,2)*Power(Cos(phi2),2) + 
                   Power(v2,2)*Power(Sin(phi2),2))))/
            Sqrt(Power(Power(v1,2)*Power(Ybtm,2) - 
               Power(Ytop,2)*(Power(v2,2)*Power(Cos(phi2),2) + 
                  Power(v2,2)*Power(Sin(phi2),2)),2)),2) - 
        ((2*Power(Ybtm,2) - (2*Power(Ybtm,2)*
                (Power(v1,2)*Power(Ybtm,2) - 
                  Power(Ytop,2)*
                   (Power(v2,2)*Power(Cos(phi2),2) + 
                     Power(v2,2)*Power(Sin(phi2),2))))/
              Sqrt(Power(Power(v1,2)*Power(Ybtm,2) - 
                 Power(Ytop,2)*(Power(v2,2)*Power(Cos(phi2),2) + 
                    Power(v2,2)*Power(Sin(phi2),2)),2)))*
           (Power(v1,2)*Power(Ybtm,2) + 
             Power(v2,2)*Power(Ytop,2)*Power(Cos(phi2),2) + 
             Power(v2,2)*Power(Ytop,2)*Power(Sin(phi2),2) - 
             Sqrt(Power(Power(v1,2)*Power(Ybtm,2) - 
                Power(Ytop,2)*(Power(v2,2)*Power(Cos(phi2),2) + 
                   Power(v2,2)*Power(Sin(phi2),2)),2))))/2. - 
        (-1.5 + Log((Power(v1,2)*Power(Ybtm,2) + 
               Power(v2,2)*Power(Ytop,2)*Power(Cos(phi2),2) + 
               Power(v2,2)*Power(Ytop,2)*Power(Sin(phi2),2) - 
               Sqrt(Power(Power(v1,2)*Power(Ybtm,2) - 
                  Power(Ytop,2)*
                   (Power(v2,2)*Power(Cos(phi2),2) + 
                     Power(v2,2)*Power(Sin(phi2),2)),2)))/(2.*Power(Rsc,2))))*
         (2*Power(Ybtm,2) - (2*Power(Ybtm,2)*
              (Power(v1,2)*Power(Ybtm,2) - 
                Power(Ytop,2)*(Power(v2,2)*Power(Cos(phi2),2) + 
                   Power(v2,2)*Power(Sin(phi2),2))))/
            Sqrt(Power(Power(v1,2)*Power(Ybtm,2) - 
               Power(Ytop,2)*(Power(v2,2)*Power(Cos(phi2),2) + 
                  Power(v2,2)*Power(Sin(phi2),2)),2)))*
         (Power(v1,2)*Power(Ybtm,2) + 
           Power(v2,2)*Power(Ytop,2)*Power(Cos(phi2),2) + 
           Power(v2,2)*Power(Ytop,2)*Power(Sin(phi2),2) - 
           Sqrt(Power(Power(v1,2)*Power(Ybtm,2) - 
              Power(Ytop,2)*(Power(v2,2)*Power(Cos(phi2),2) + 
                 Power(v2,2)*Power(Sin(phi2),2)),2))) - 
        ((2*Power(Ybtm,2) + (2*Power(Ybtm,2)*
                (Power(v1,2)*Power(Ybtm,2) - 
                  Power(Ytop,2)*
                   (Power(v2,2)*Power(Cos(phi2),2) + 
                     Power(v2,2)*Power(Sin(phi2),2))))/
              Sqrt(Power(Power(v1,2)*Power(Ybtm,2) - 
                 Power(Ytop,2)*(Power(v2,2)*Power(Cos(phi2),2) + 
                    Power(v2,2)*Power(Sin(phi2),2)),2)))*
           (Power(v1,2)*Power(Ybtm,2) + 
             Power(v2,2)*Power(Ytop,2)*Power(Cos(phi2),2) + 
             Power(v2,2)*Power(Ytop,2)*Power(Sin(phi2),2) + 
             Sqrt(Power(Power(v1,2)*Power(Ybtm,2) - 
                Power(Ytop,2)*(Power(v2,2)*Power(Cos(phi2),2) + 
                   Power(v2,2)*Power(Sin(phi2),2)),2))))/2. - 
        (-1.5 + Log((Power(v1,2)*Power(Ybtm,2) + 
               Power(v2,2)*Power(Ytop,2)*Power(Cos(phi2),2) + 
               Power(v2,2)*Power(Ytop,2)*Power(Sin(phi2),2) + 
               Sqrt(Power(Power(v1,2)*Power(Ybtm,2) - 
                  Power(Ytop,2)*
                   (Power(v2,2)*Power(Cos(phi2),2) + 
                     Power(v2,2)*Power(Sin(phi2),2)),2)))/(2.*Power(Rsc,2))))*
         (2*Power(Ybtm,2) + (2*Power(Ybtm,2)*
              (Power(v1,2)*Power(Ybtm,2) - 
                Power(Ytop,2)*(Power(v2,2)*Power(Cos(phi2),2) + 
                   Power(v2,2)*Power(Sin(phi2),2))))/
            Sqrt(Power(Power(v1,2)*Power(Ybtm,2) - 
               Power(Ytop,2)*(Power(v2,2)*Power(Cos(phi2),2) + 
                  Power(v2,2)*Power(Sin(phi2),2)),2)))*
         (Power(v1,2)*Power(Ybtm,2) + 
           Power(v2,2)*Power(Ytop,2)*Power(Cos(phi2),2) + 
           Power(v2,2)*Power(Ytop,2)*Power(Sin(phi2),2) + 
           Sqrt(Power(Power(v1,2)*Power(Ybtm,2) - 
              Power(Ytop,2)*(Power(v2,2)*Power(Cos(phi2),2) + 
                 Power(v2,2)*Power(Sin(phi2),2)),2))) + 
        (3*Power(Complex(0,4)*lambda*Power(Ytop,2)*
              (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))*
              (Atop*(v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) - 
                Complex(0,1)*lambda*v1*
                 (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS))) - 
             Complex(0,4)*lambda*Power(Ytop,2)*
              (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS))*
              (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) + 
                Complex(0,1)*lambda*v1*
                 (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))),2))/
         (8.*(Power(Power(MQ33,2) - Power(MU33,2),2) + 
             4*Power(Ytop,2)*(Atop*
                 (v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) - 
                Complex(0,1)*lambda*v1*
                 (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS)))*
              (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) + 
                Complex(0,1)*lambda*v1*
                 (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))))) + 
        ((-1.5 + Log((Power(MQ33,2) + Power(MU33,2) + 
                 2*Power(Ytop,2)*
                  (Power(v2,2)*Power(Cos(phi2),2) + 
                    Power(v2,2)*Power(Sin(phi2),2)) - 
                 Sqrt(Power(Power(MQ33,2) - Power(MU33,2),2) + 
                   4*Power(Ytop,2)*
                    (Atop*(v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) - 
                      Complex(0,1)*lambda*v1*
                       (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS)))*
                    (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) + 
                      Complex(0,1)*lambda*v1*
                       (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS)))))/
               (2.*Power(Rsc,2))))*
           Power(Complex(0,4)*lambda*Power(Ytop,2)*
              (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))*
              (Atop*(v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) - 
                Complex(0,1)*lambda*v1*
                 (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS))) - 
             Complex(0,4)*lambda*Power(Ytop,2)*
              (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS))*
              (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) + 
                Complex(0,1)*lambda*v1*
                 (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))),2))/
         (8.*(Power(Power(MQ33,2) - Power(MU33,2),2) + 
             4*Power(Ytop,2)*(Atop*
                 (v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) - 
                Complex(0,1)*lambda*v1*
                 (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS)))*
              (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) + 
                Complex(0,1)*lambda*v1*
                 (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))))) + 
        ((-1.5 + Log((Power(MQ33,2) + Power(MU33,2) + 
                 2*Power(Ytop,2)*
                  (Power(v2,2)*Power(Cos(phi2),2) + 
                    Power(v2,2)*Power(Sin(phi2),2)) + 
                 Sqrt(Power(Power(MQ33,2) - Power(MU33,2),2) + 
                   4*Power(Ytop,2)*
                    (Atop*(v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) - 
                      Complex(0,1)*lambda*v1*
                       (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS)))*
                    (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) + 
                      Complex(0,1)*lambda*v1*
                       (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS)))))/
               (2.*Power(Rsc,2))))*
           Power(Complex(0,4)*lambda*Power(Ytop,2)*
              (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))*
              (Atop*(v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) - 
                Complex(0,1)*lambda*v1*
                 (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS))) - 
             Complex(0,4)*lambda*Power(Ytop,2)*
              (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS))*
              (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) + 
                Complex(0,1)*lambda*v1*
                 (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))),2))/
         (8.*(Power(Power(MQ33,2) - Power(MU33,2),2) + 
             4*Power(Ytop,2)*(Atop*
                 (v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) - 
                Complex(0,1)*lambda*v1*
                 (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS)))*
              (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) + 
                Complex(0,1)*lambda*v1*
                 (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))))) + 
        (Power(Complex(0,4)*lambda*Power(Ytop,2)*
              (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))*
              (Atop*(v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) - 
                Complex(0,1)*lambda*v1*
                 (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS))) - 
             Complex(0,4)*lambda*Power(Ytop,2)*
              (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS))*
              (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) + 
                Complex(0,1)*lambda*v1*
                 (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))),2)*
           (Power(MQ33,2) + Power(MU33,2) + 
             2*Power(Ytop,2)*(Power(v2,2)*Power(Cos(phi2),2) + 
                Power(v2,2)*Power(Sin(phi2),2)) - 
             Sqrt(Power(Power(MQ33,2) - Power(MU33,2),2) + 
               4*Power(Ytop,2)*(Atop*
                   (v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) - 
                  Complex(0,1)*lambda*v1*
                   (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS)))*
                (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) + 
                  Complex(0,1)*lambda*v1*
                   (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))))))/
         (16.*Power(Power(Power(MQ33,2) - Power(MU33,2),2) + 
             4*Power(Ytop,2)*(Atop*
                 (v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) - 
                Complex(0,1)*lambda*v1*
                 (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS)))*
              (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) + 
                Complex(0,1)*lambda*v1*
                 (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))),1.5)) + 
        ((-1.5 + Log((Power(MQ33,2) + Power(MU33,2) + 
                 2*Power(Ytop,2)*
                  (Power(v2,2)*Power(Cos(phi2),2) + 
                    Power(v2,2)*Power(Sin(phi2),2)) - 
                 Sqrt(Power(Power(MQ33,2) - Power(MU33,2),2) + 
                   4*Power(Ytop,2)*
                    (Atop*(v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) - 
                      Complex(0,1)*lambda*v1*
                       (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS)))*
                    (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) + 
                      Complex(0,1)*lambda*v1*
                       (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS)))))/
               (2.*Power(Rsc,2))))*
           Power(Complex(0,4)*lambda*Power(Ytop,2)*
              (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))*
              (Atop*(v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) - 
                Complex(0,1)*lambda*v1*
                 (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS))) - 
             Complex(0,4)*lambda*Power(Ytop,2)*
              (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS))*
              (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) + 
                Complex(0,1)*lambda*v1*
                 (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))),2)*
           (Power(MQ33,2) + Power(MU33,2) + 
             2*Power(Ytop,2)*(Power(v2,2)*Power(Cos(phi2),2) + 
                Power(v2,2)*Power(Sin(phi2),2)) - 
             Sqrt(Power(Power(MQ33,2) - Power(MU33,2),2) + 
               4*Power(Ytop,2)*(Atop*
                   (v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) - 
                  Complex(0,1)*lambda*v1*
                   (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS)))*
                (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) + 
                  Complex(0,1)*lambda*v1*
                   (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))))))/
         (8.*Power(Power(Power(MQ33,2) - Power(MU33,2),2) + 
             4*Power(Ytop,2)*(Atop*
                 (v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) - 
                Complex(0,1)*lambda*v1*
                 (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS)))*
              (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) + 
                Complex(0,1)*lambda*v1*
                 (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))),1.5)) - 
        (Power(lambda,2)*Power(Ytop,2)*
           (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS))*
           (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))*
           (Power(MQ33,2) + Power(MU33,2) + 
             2*Power(Ytop,2)*(Power(v2,2)*Power(Cos(phi2),2) + 
                Power(v2,2)*Power(Sin(phi2),2)) - 
             Sqrt(Power(Power(MQ33,2) - Power(MU33,2),2) + 
               4*Power(Ytop,2)*(Atop*
                   (v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) - 
                  Complex(0,1)*lambda*v1*
                   (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS)))*
                (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) + 
                  Complex(0,1)*lambda*v1*
                   (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))))))/
         Sqrt(Power(Power(MQ33,2) - Power(MU33,2),2) + 
           4*Power(Ytop,2)*(Atop*(v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) - 
              Complex(0,1)*lambda*v1*
               (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS)))*
            (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) + 
              Complex(0,1)*lambda*v1*
               (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS)))) - 
        (2*Power(lambda,2)*Power(Ytop,2)*
           (-1.5 + Log((Power(MQ33,2) + Power(MU33,2) + 
                 2*Power(Ytop,2)*
                  (Power(v2,2)*Power(Cos(phi2),2) + 
                    Power(v2,2)*Power(Sin(phi2),2)) - 
                 Sqrt(Power(Power(MQ33,2) - Power(MU33,2),2) + 
                   4*Power(Ytop,2)*
                    (Atop*(v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) - 
                      Complex(0,1)*lambda*v1*
                       (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS)))*
                    (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) + 
                      Complex(0,1)*lambda*v1*
                       (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS)))))/
               (2.*Power(Rsc,2))))*
           (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS))*
           (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))*
           (Power(MQ33,2) + Power(MU33,2) + 
             2*Power(Ytop,2)*(Power(v2,2)*Power(Cos(phi2),2) + 
                Power(v2,2)*Power(Sin(phi2),2)) - 
             Sqrt(Power(Power(MQ33,2) - Power(MU33,2),2) + 
               4*Power(Ytop,2)*(Atop*
                   (v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) - 
                  Complex(0,1)*lambda*v1*
                   (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS)))*
                (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) + 
                  Complex(0,1)*lambda*v1*
                   (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))))))/
         Sqrt(Power(Power(MQ33,2) - Power(MU33,2),2) + 
           4*Power(Ytop,2)*(Atop*(v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) - 
              Complex(0,1)*lambda*v1*
               (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS)))*
            (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) + 
              Complex(0,1)*lambda*v1*
               (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS)))) - 
        (Power(Complex(0,4)*lambda*Power(Ytop,2)*
              (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))*
              (Atop*(v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) - 
                Complex(0,1)*lambda*v1*
                 (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS))) - 
             Complex(0,4)*lambda*Power(Ytop,2)*
              (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS))*
              (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) + 
                Complex(0,1)*lambda*v1*
                 (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))),2)*
           (Power(MQ33,2) + Power(MU33,2) + 
             2*Power(Ytop,2)*(Power(v2,2)*Power(Cos(phi2),2) + 
                Power(v2,2)*Power(Sin(phi2),2)) + 
             Sqrt(Power(Power(MQ33,2) - Power(MU33,2),2) + 
               4*Power(Ytop,2)*(Atop*
                   (v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) - 
                  Complex(0,1)*lambda*v1*
                   (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS)))*
                (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) + 
                  Complex(0,1)*lambda*v1*
                   (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))))))/
         (16.*Power(Power(Power(MQ33,2) - Power(MU33,2),2) + 
             4*Power(Ytop,2)*(Atop*
                 (v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) - 
                Complex(0,1)*lambda*v1*
                 (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS)))*
              (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) + 
                Complex(0,1)*lambda*v1*
                 (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))),1.5)) - 
        ((-1.5 + Log((Power(MQ33,2) + Power(MU33,2) + 
                 2*Power(Ytop,2)*
                  (Power(v2,2)*Power(Cos(phi2),2) + 
                    Power(v2,2)*Power(Sin(phi2),2)) + 
                 Sqrt(Power(Power(MQ33,2) - Power(MU33,2),2) + 
                   4*Power(Ytop,2)*
                    (Atop*(v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) - 
                      Complex(0,1)*lambda*v1*
                       (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS)))*
                    (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) + 
                      Complex(0,1)*lambda*v1*
                       (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS)))))/
               (2.*Power(Rsc,2))))*
           Power(Complex(0,4)*lambda*Power(Ytop,2)*
              (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))*
              (Atop*(v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) - 
                Complex(0,1)*lambda*v1*
                 (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS))) - 
             Complex(0,4)*lambda*Power(Ytop,2)*
              (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS))*
              (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) + 
                Complex(0,1)*lambda*v1*
                 (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))),2)*
           (Power(MQ33,2) + Power(MU33,2) + 
             2*Power(Ytop,2)*(Power(v2,2)*Power(Cos(phi2),2) + 
                Power(v2,2)*Power(Sin(phi2),2)) + 
             Sqrt(Power(Power(MQ33,2) - Power(MU33,2),2) + 
               4*Power(Ytop,2)*(Atop*
                   (v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) - 
                  Complex(0,1)*lambda*v1*
                   (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS)))*
                (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) + 
                  Complex(0,1)*lambda*v1*
                   (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))))))/
         (8.*Power(Power(Power(MQ33,2) - Power(MU33,2),2) + 
             4*Power(Ytop,2)*(Atop*
                 (v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) - 
                Complex(0,1)*lambda*v1*
                 (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS)))*
              (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) + 
                Complex(0,1)*lambda*v1*
                 (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))),1.5)) + 
        (Power(lambda,2)*Power(Ytop,2)*
           (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS))*
           (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))*
           (Power(MQ33,2) + Power(MU33,2) + 
             2*Power(Ytop,2)*(Power(v2,2)*Power(Cos(phi2),2) + 
                Power(v2,2)*Power(Sin(phi2),2)) + 
             Sqrt(Power(Power(MQ33,2) - Power(MU33,2),2) + 
               4*Power(Ytop,2)*(Atop*
                   (v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) - 
                  Complex(0,1)*lambda*v1*
                   (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS)))*
                (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) + 
                  Complex(0,1)*lambda*v1*
                   (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))))))/
         Sqrt(Power(Power(MQ33,2) - Power(MU33,2),2) + 
           4*Power(Ytop,2)*(Atop*(v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) - 
              Complex(0,1)*lambda*v1*
               (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS)))*
            (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) + 
              Complex(0,1)*lambda*v1*
               (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS)))) + 
        (2*Power(lambda,2)*Power(Ytop,2)*
           (-1.5 + Log((Power(MQ33,2) + Power(MU33,2) + 
                 2*Power(Ytop,2)*
                  (Power(v2,2)*Power(Cos(phi2),2) + 
                    Power(v2,2)*Power(Sin(phi2),2)) + 
                 Sqrt(Power(Power(MQ33,2) - Power(MU33,2),2) + 
                   4*Power(Ytop,2)*
                    (Atop*(v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) - 
                      Complex(0,1)*lambda*v1*
                       (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS)))*
                    (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) + 
                      Complex(0,1)*lambda*v1*
                       (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS)))))/
               (2.*Power(Rsc,2))))*
           (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS))*
           (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))*
           (Power(MQ33,2) + Power(MU33,2) + 
             2*Power(Ytop,2)*(Power(v2,2)*Power(Cos(phi2),2) + 
                Power(v2,2)*Power(Sin(phi2),2)) + 
             Sqrt(Power(Power(MQ33,2) - Power(MU33,2),2) + 
               4*Power(Ytop,2)*(Atop*
                   (v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) - 
                  Complex(0,1)*lambda*v1*
                   (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS)))*
                (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) + 
                  Complex(0,1)*lambda*v1*
                   (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS))))))/
         Sqrt(Power(Power(MQ33,2) - Power(MU33,2),2) + 
           4*Power(Ytop,2)*(Atop*(v2*Cos(phi2) - Complex(0,1)*v2*Sin(phi2)) - 
              Complex(0,1)*lambda*v1*
               (Complex(0,-1)*vevS*Cos(phiS) + vevS*Sin(phiS)))*
            (Atop*(v2*Cos(phi2) + Complex(0,1)*v2*Sin(phi2)) + 
              Complex(0,1)*lambda*v1*
               (Complex(0,1)*vevS*Cos(phiS) + vevS*Sin(phiS)))) + 
        (3*Power(4*v1*Power(Ybtm,2) - 
             (2*Power(Ybtm,2)*(2*Power(Abtm,2)*v1 + 
                  2*Abtm*lambda*
                   (-(v2*vevS*Cos(phi2)*Cos(phiS)) + 
                     v2*vevS*Sin(phi2)*Sin(phiS))))/
              Sqrt(Power(MD33,4) - 2*Power(MD33,2)*Power(MQ33,2) + 
                Power(MQ33,4) + 
                4*Power(Ybtm,2)*
                 (Power(Abtm,2)*Power(v1,2) + 
                   2*Abtm*lambda*
                    (-(v1*v2*vevS*Cos(phi2)*Cos(phiS)) + 
                      v1*v2*vevS*Sin(phi2)*Sin(phiS)) + 
                   Power(lambda,2)*
                    (Power(v2,2)*Power(Cos(phi2),2) + 
                      Power(v2,2)*Power(Sin(phi2),2))*
                    (Power(vevS,2)*Power(Cos(phiS),2) + 
                      Power(vevS,2)*Power(Sin(phiS),2)))),2))/4. + 
        ((-1.5 + Log((Power(MD33,2) + Power(MQ33,2) + 
                 2*Power(v1,2)*Power(Ybtm,2) - 
                 Sqrt(Power(MD33,4) - 2*Power(MD33,2)*Power(MQ33,2) + 
                   Power(MQ33,4) + 
                   4*Power(Ybtm,2)*
                    (Power(Abtm,2)*Power(v1,2) + 
                      2*Abtm*lambda*
                       (-(v1*v2*vevS*Cos(phi2)*Cos(phiS)) + 
                         v1*v2*vevS*Sin(phi2)*Sin(phiS)) + 
                      Power(lambda,2)*
                       (Power(v2,2)*Power(Cos(phi2),2) + 
                         Power(v2,2)*Power(Sin(phi2),2))*
                       (Power(vevS,2)*Power(Cos(phiS),2) + 
                         Power(vevS,2)*Power(Sin(phiS),2)))))/
               (2.*Power(Rsc,2))))*
           Power(4*v1*Power(Ybtm,2) - 
             (2*Power(Ybtm,2)*(2*Power(Abtm,2)*v1 + 
                  2*Abtm*lambda*
                   (-(v2*vevS*Cos(phi2)*Cos(phiS)) + 
                     v2*vevS*Sin(phi2)*Sin(phiS))))/
              Sqrt(Power(MD33,4) - 2*Power(MD33,2)*Power(MQ33,2) + 
                Power(MQ33,4) + 
                4*Power(Ybtm,2)*
                 (Power(Abtm,2)*Power(v1,2) + 
                   2*Abtm*lambda*
                    (-(v1*v2*vevS*Cos(phi2)*Cos(phiS)) + 
                      v1*v2*vevS*Sin(phi2)*Sin(phiS)) + 
                   Power(lambda,2)*
                    (Power(v2,2)*Power(Cos(phi2),2) + 
                      Power(v2,2)*Power(Sin(phi2),2))*
                    (Power(vevS,2)*Power(Cos(phiS),2) + 
                      Power(vevS,2)*Power(Sin(phiS),2)))),2))/2. + 
        (3*Power(4*v1*Power(Ybtm,2) + 
             (2*Power(Ybtm,2)*(2*Power(Abtm,2)*v1 + 
                  2*Abtm*lambda*
                   (-(v2*vevS*Cos(phi2)*Cos(phiS)) + 
                     v2*vevS*Sin(phi2)*Sin(phiS))))/
              Sqrt(Power(MD33,4) - 2*Power(MD33,2)*Power(MQ33,2) + 
                Power(MQ33,4) + 
                4*Power(Ybtm,2)*
                 (Power(Abtm,2)*Power(v1,2) + 
                   2*Abtm*lambda*
                    (-(v1*v2*vevS*Cos(phi2)*Cos(phiS)) + 
                      v1*v2*vevS*Sin(phi2)*Sin(phiS)) + 
                   Power(lambda,2)*
                    (Power(v2,2)*Power(Cos(phi2),2) + 
                      Power(v2,2)*Power(Sin(phi2),2))*
                    (Power(vevS,2)*Power(Cos(phiS),2) + 
                      Power(vevS,2)*Power(Sin(phiS),2)))),2))/4. + 
        ((-1.5 + Log((Power(MD33,2) + Power(MQ33,2) + 
                 2*Power(v1,2)*Power(Ybtm,2) + 
                 Sqrt(Power(MD33,4) - 2*Power(MD33,2)*Power(MQ33,2) + 
                   Power(MQ33,4) + 
                   4*Power(Ybtm,2)*
                    (Power(Abtm,2)*Power(v1,2) + 
                      2*Abtm*lambda*
                       (-(v1*v2*vevS*Cos(phi2)*Cos(phiS)) + 
                         v1*v2*vevS*Sin(phi2)*Sin(phiS)) + 
                      Power(lambda,2)*
                       (Power(v2,2)*Power(Cos(phi2),2) + 
                         Power(v2,2)*Power(Sin(phi2),2))*
                       (Power(vevS,2)*Power(Cos(phiS),2) + 
                         Power(vevS,2)*Power(Sin(phiS),2)))))/
               (2.*Power(Rsc,2))))*
           Power(4*v1*Power(Ybtm,2) + 
             (2*Power(Ybtm,2)*(2*Power(Abtm,2)*v1 + 
                  2*Abtm*lambda*
                   (-(v2*vevS*Cos(phi2)*Cos(phiS)) + 
                     v2*vevS*Sin(phi2)*Sin(phiS))))/
              Sqrt(Power(MD33,4) - 2*Power(MD33,2)*Power(MQ33,2) + 
                Power(MQ33,4) + 
                4*Power(Ybtm,2)*
                 (Power(Abtm,2)*Power(v1,2) + 
                   2*Abtm*lambda*
                    (-(v1*v2*vevS*Cos(phi2)*Cos(phiS)) + 
                      v1*v2*vevS*Sin(phi2)*Sin(phiS)) + 
                   Power(lambda,2)*
                    (Power(v2,2)*Power(Cos(phi2),2) + 
                      Power(v2,2)*Power(Sin(phi2),2))*
                    (Power(vevS,2)*Power(Cos(phiS),2) + 
                      Power(vevS,2)*Power(Sin(phiS),2)))),2))/2. + 
        ((4*Power(Ybtm,2) + (4*Power(Ybtm,4)*
                Power(2*Power(Abtm,2)*v1 + 
                  2*Abtm*lambda*
                   (-(v2*vevS*Cos(phi2)*Cos(phiS)) + 
                     v2*vevS*Sin(phi2)*Sin(phiS)),2))/
              Power(Power(MD33,4) - 2*Power(MD33,2)*Power(MQ33,2) + 
                Power(MQ33,4) + 
                4*Power(Ybtm,2)*
                 (Power(Abtm,2)*Power(v1,2) + 
                   2*Abtm*lambda*
                    (-(v1*v2*vevS*Cos(phi2)*Cos(phiS)) + 
                      v1*v2*vevS*Sin(phi2)*Sin(phiS)) + 
                   Power(lambda,2)*
                    (Power(v2,2)*Power(Cos(phi2),2) + 
                      Power(v2,2)*Power(Sin(phi2),2))*
                    (Power(vevS,2)*Power(Cos(phiS),2) + 
                      Power(vevS,2)*Power(Sin(phiS),2))),1.5) - 
             (4*Power(Abtm,2)*Power(Ybtm,2))/
              Sqrt(Power(MD33,4) - 2*Power(MD33,2)*Power(MQ33,2) + 
                Power(MQ33,4) + 
                4*Power(Ybtm,2)*
                 (Power(Abtm,2)*Power(v1,2) + 
                   2*Abtm*lambda*
                    (-(v1*v2*vevS*Cos(phi2)*Cos(phiS)) + 
                      v1*v2*vevS*Sin(phi2)*Sin(phiS)) + 
                   Power(lambda,2)*
                    (Power(v2,2)*Power(Cos(phi2),2) + 
                      Power(v2,2)*Power(Sin(phi2),2))*
                    (Power(vevS,2)*Power(Cos(phiS),2) + 
                      Power(vevS,2)*Power(Sin(phiS),2)))))*
           (Power(MD33,2) + Power(MQ33,2) + 2*Power(v1,2)*Power(Ybtm,2) - 
             Sqrt(Power(MD33,4) - 2*Power(MD33,2)*Power(MQ33,2) + 
               Power(MQ33,4) + 4*Power(Ybtm,2)*
                (Power(Abtm,2)*Power(v1,2) + 
                  2*Abtm*lambda*
                   (-(v1*v2*vevS*Cos(phi2)*Cos(phiS)) + 
                     v1*v2*vevS*Sin(phi2)*Sin(phiS)) + 
                  Power(lambda,2)*
                   (Power(v2,2)*Power(Cos(phi2),2) + 
                     Power(v2,2)*Power(Sin(phi2),2))*
                   (Power(vevS,2)*Power(Cos(phiS),2) + 
                     Power(vevS,2)*Power(Sin(phiS),2))))))/4. + 
        ((-1.5 + Log((Power(MD33,2) + Power(MQ33,2) + 
                 2*Power(v1,2)*Power(Ybtm,2) - 
                 Sqrt(Power(MD33,4) - 2*Power(MD33,2)*Power(MQ33,2) + 
                   Power(MQ33,4) + 
                   4*Power(Ybtm,2)*
                    (Power(Abtm,2)*Power(v1,2) + 
                      2*Abtm*lambda*
                       (-(v1*v2*vevS*Cos(phi2)*Cos(phiS)) + 
                         v1*v2*vevS*Sin(phi2)*Sin(phiS)) + 
                      Power(lambda,2)*
                       (Power(v2,2)*Power(Cos(phi2),2) + 
                         Power(v2,2)*Power(Sin(phi2),2))*
                       (Power(vevS,2)*Power(Cos(phiS),2) + 
                         Power(vevS,2)*Power(Sin(phiS),2)))))/
               (2.*Power(Rsc,2))))*
           (4*Power(Ybtm,2) + (4*Power(Ybtm,4)*
                Power(2*Power(Abtm,2)*v1 + 
                  2*Abtm*lambda*
                   (-(v2*vevS*Cos(phi2)*Cos(phiS)) + 
                     v2*vevS*Sin(phi2)*Sin(phiS)),2))/
              Power(Power(MD33,4) - 2*Power(MD33,2)*Power(MQ33,2) + 
                Power(MQ33,4) + 
                4*Power(Ybtm,2)*
                 (Power(Abtm,2)*Power(v1,2) + 
                   2*Abtm*lambda*
                    (-(v1*v2*vevS*Cos(phi2)*Cos(phiS)) + 
                      v1*v2*vevS*Sin(phi2)*Sin(phiS)) + 
                   Power(lambda,2)*
                    (Power(v2,2)*Power(Cos(phi2),2) + 
                      Power(v2,2)*Power(Sin(phi2),2))*
                    (Power(vevS,2)*Power(Cos(phiS),2) + 
                      Power(vevS,2)*Power(Sin(phiS),2))),1.5) - 
             (4*Power(Abtm,2)*Power(Ybtm,2))/
              Sqrt(Power(MD33,4) - 2*Power(MD33,2)*Power(MQ33,2) + 
                Power(MQ33,4) + 
                4*Power(Ybtm,2)*
                 (Power(Abtm,2)*Power(v1,2) + 
                   2*Abtm*lambda*
                    (-(v1*v2*vevS*Cos(phi2)*Cos(phiS)) + 
                      v1*v2*vevS*Sin(phi2)*Sin(phiS)) + 
                   Power(lambda,2)*
                    (Power(v2,2)*Power(Cos(phi2),2) + 
                      Power(v2,2)*Power(Sin(phi2),2))*
                    (Power(vevS,2)*Power(Cos(phiS),2) + 
                      Power(vevS,2)*Power(Sin(phiS),2)))))*
           (Power(MD33,2) + Power(MQ33,2) + 2*Power(v1,2)*Power(Ybtm,2) - 
             Sqrt(Power(MD33,4) - 2*Power(MD33,2)*Power(MQ33,2) + 
               Power(MQ33,4) + 4*Power(Ybtm,2)*
                (Power(Abtm,2)*Power(v1,2) + 
                  2*Abtm*lambda*
                   (-(v1*v2*vevS*Cos(phi2)*Cos(phiS)) + 
                     v1*v2*vevS*Sin(phi2)*Sin(phiS)) + 
                  Power(lambda,2)*
                   (Power(v2,2)*Power(Cos(phi2),2) + 
                     Power(v2,2)*Power(Sin(phi2),2))*
                   (Power(vevS,2)*Power(Cos(phiS),2) + 
                     Power(vevS,2)*Power(Sin(phiS),2))))))/2. + 
        ((4*Power(Ybtm,2) - (4*Power(Ybtm,4)*
                Power(2*Power(Abtm,2)*v1 + 
                  2*Abtm*lambda*
                   (-(v2*vevS*Cos(phi2)*Cos(phiS)) + 
                     v2*vevS*Sin(phi2)*Sin(phiS)),2))/
              Power(Power(MD33,4) - 2*Power(MD33,2)*Power(MQ33,2) + 
                Power(MQ33,4) + 
                4*Power(Ybtm,2)*
                 (Power(Abtm,2)*Power(v1,2) + 
                   2*Abtm*lambda*
                    (-(v1*v2*vevS*Cos(phi2)*Cos(phiS)) + 
                      v1*v2*vevS*Sin(phi2)*Sin(phiS)) + 
                   Power(lambda,2)*
                    (Power(v2,2)*Power(Cos(phi2),2) + 
                      Power(v2,2)*Power(Sin(phi2),2))*
                    (Power(vevS,2)*Power(Cos(phiS),2) + 
                      Power(vevS,2)*Power(Sin(phiS),2))),1.5) + 
             (4*Power(Abtm,2)*Power(Ybtm,2))/
              Sqrt(Power(MD33,4) - 2*Power(MD33,2)*Power(MQ33,2) + 
                Power(MQ33,4) + 
                4*Power(Ybtm,2)*
                 (Power(Abtm,2)*Power(v1,2) + 
                   2*Abtm*lambda*
                    (-(v1*v2*vevS*Cos(phi2)*Cos(phiS)) + 
                      v1*v2*vevS*Sin(phi2)*Sin(phiS)) + 
                   Power(lambda,2)*
                    (Power(v2,2)*Power(Cos(phi2),2) + 
                      Power(v2,2)*Power(Sin(phi2),2))*
                    (Power(vevS,2)*Power(Cos(phiS),2) + 
                      Power(vevS,2)*Power(Sin(phiS),2)))))*
           (Power(MD33,2) + Power(MQ33,2) + 2*Power(v1,2)*Power(Ybtm,2) + 
             Sqrt(Power(MD33,4) - 2*Power(MD33,2)*Power(MQ33,2) + 
               Power(MQ33,4) + 4*Power(Ybtm,2)*
                (Power(Abtm,2)*Power(v1,2) + 
                  2*Abtm*lambda*
                   (-(v1*v2*vevS*Cos(phi2)*Cos(phiS)) + 
                     v1*v2*vevS*Sin(phi2)*Sin(phiS)) + 
                  Power(lambda,2)*
                   (Power(v2,2)*Power(Cos(phi2),2) + 
                     Power(v2,2)*Power(Sin(phi2),2))*
                   (Power(vevS,2)*Power(Cos(phiS),2) + 
                     Power(vevS,2)*Power(Sin(phiS),2))))))/4. + 
        ((-1.5 + Log((Power(MD33,2) + Power(MQ33,2) + 
                 2*Power(v1,2)*Power(Ybtm,2) + 
                 Sqrt(Power(MD33,4) - 2*Power(MD33,2)*Power(MQ33,2) + 
                   Power(MQ33,4) + 
                   4*Power(Ybtm,2)*
                    (Power(Abtm,2)*Power(v1,2) + 
                      2*Abtm*lambda*
                       (-(v1*v2*vevS*Cos(phi2)*Cos(phiS)) + 
                         v1*v2*vevS*Sin(phi2)*Sin(phiS)) + 
                      Power(lambda,2)*
                       (Power(v2,2)*Power(Cos(phi2),2) + 
                         Power(v2,2)*Power(Sin(phi2),2))*
                       (Power(vevS,2)*Power(Cos(phiS),2) + 
                         Power(vevS,2)*Power(Sin(phiS),2)))))/
               (2.*Power(Rsc,2))))*
           (4*Power(Ybtm,2) - (4*Power(Ybtm,4)*
                Power(2*Power(Abtm,2)*v1 + 
                  2*Abtm*lambda*
                   (-(v2*vevS*Cos(phi2)*Cos(phiS)) + 
                     v2*vevS*Sin(phi2)*Sin(phiS)),2))/
              Power(Power(MD33,4) - 2*Power(MD33,2)*Power(MQ33,2) + 
                Power(MQ33,4) + 
                4*Power(Ybtm,2)*
                 (Power(Abtm,2)*Power(v1,2) + 
                   2*Abtm*lambda*
                    (-(v1*v2*vevS*Cos(phi2)*Cos(phiS)) + 
                      v1*v2*vevS*Sin(phi2)*Sin(phiS)) + 
                   Power(lambda,2)*
                    (Power(v2,2)*Power(Cos(phi2),2) + 
                      Power(v2,2)*Power(Sin(phi2),2))*
                    (Power(vevS,2)*Power(Cos(phiS),2) + 
                      Power(vevS,2)*Power(Sin(phiS),2))),1.5) + 
             (4*Power(Abtm,2)*Power(Ybtm,2))/
              Sqrt(Power(MD33,4) - 2*Power(MD33,2)*Power(MQ33,2) + 
                Power(MQ33,4) + 
                4*Power(Ybtm,2)*
                 (Power(Abtm,2)*Power(v1,2) + 
                   2*Abtm*lambda*
                    (-(v1*v2*vevS*Cos(phi2)*Cos(phiS)) + 
                      v1*v2*vevS*Sin(phi2)*Sin(phiS)) + 
                   Power(lambda,2)*
                    (Power(v2,2)*Power(Cos(phi2),2) + 
                      Power(v2,2)*Power(Sin(phi2),2))*
                    (Power(vevS,2)*Power(Cos(phiS),2) + 
                      Power(vevS,2)*Power(Sin(phiS),2)))))*
           (Power(MD33,2) + Power(MQ33,2) + 2*Power(v1,2)*Power(Ybtm,2) + 
             Sqrt(Power(MD33,4) - 2*Power(MD33,2)*Power(MQ33,2) + 
               Power(MQ33,4) + 4*Power(Ybtm,2)*
                (Power(Abtm,2)*Power(v1,2) + 
                  2*Abtm*lambda*
                   (-(v1*v2*vevS*Cos(phi2)*Cos(phiS)) + 
                     v1*v2*vevS*Sin(phi2)*Sin(phiS)) + 
                  Power(lambda,2)*
                   (Power(v2,2)*Power(Cos(phi2),2) + 
                     Power(v2,2)*Power(Sin(phi2),2))*
                   (Power(vevS,2)*Power(Cos(phiS),2) + 
                     Power(vevS,2)*Power(Sin(phiS),2))))))/2.))/
    (64.*Power(Pi,2))