MX11=(4.d0*mH1**2+3.d0*g2**2*v1**2-g2**2*v2**2+4.d0*lambda**2*
     . v2**2+g1**2*(3.d0*v1**2-v2**2)+4.d0*lambda**2*vevS**2)/4.d0+
     . (3.d0*(-28.d0*v1**2*Ybtm**4+(24.d0*lambda**2*vevS**2*Ytop**4*
     . (lambda*v1*vevS-Atop*v2*DCOS(phi2+phiS))**2)/(MQ33**4-2.d0*
     . MQ33**2*MU33**2+MU33**4+4.d0*Atop**2*v2**2*Ytop**2+4.d0*
     . lambda**2*v1**2*vevS**2*Ytop**2-8.d0*Atop*lambda*v1*v2*vevS*
     . Ytop**2*DCOS(phi2+phiS))-12.d0*v1**2*Ybtm**4*(-3.d0+2.d0*
     . DLOG((v1**2*Ybtm**2)/Rsc**2))+(4.d0*lambda**2*vevS**2*Ytop**4*
     . (lambda*v1*vevS-Atop*v2*DCOS(phi2+phiS))**2*(-3.d0+2.d0*
     . DLOG((MQ33**2+MU33**2+2.d0*v2**2*Ytop**2-DSQRT(MQ33**4-2.d0*
     . MQ33**2*MU33**2+MU33**4+4.d0*Atop**2*v2**2*Ytop**2+4.d0*
     . lambda**2*v1**2*vevS**2*Ytop**2-8.d0*Atop*lambda*v1*v2*vevS*
     . Ytop**2*DCOS(phi2+phiS)))/(2.d0*Rsc**2))))/(MQ33**4-2.d0*MQ33**2
     . *MU33**2+MU33**4+4.d0*Atop**2*v2**2*Ytop**2+4.d0*lambda**2*v1**2
     . *vevS**2*Ytop**2-8.d0*Atop*lambda*v1*v2*vevS*Ytop**2*DCOS(phi2+
     . phiS))+(4.d0*lambda**2*vevS**2*Ytop**4*(lambda*v1*vevS-Atop*v2*
     . DCOS(phi2+phiS))**2*(-3.d0+2.d0*DLOG((MQ33**2+MU33**2+2.d0*v2**2
     . *Ytop**2+DSQRT(MQ33**4-2.d0*MQ33**2*MU33**2+MU33**4+4.d0*Atop**2
     . *v2**2*Ytop**2+4.d0*lambda**2*v1**2*vevS**2*Ytop**2-8.d0*Atop*
     . lambda*v1*v2*vevS*Ytop**2*DCOS(phi2+phiS)))/(2.d0*Rsc**2))))/
     . (MQ33**4-2.d0*MQ33**2*MU33**2+MU33**4+4.d0*Atop**2*v2**2*Ytop**2
     . +4.d0*lambda**2*v1**2*vevS**2*Ytop**2-8.d0*Atop*lambda*v1*v2*
     . vevS*Ytop**2*DCOS(phi2+phiS))+12.d0*Ybtm**4*(v1+(Abtm*(Abtm*v1-
     . lambda*v2*vevS*DCOS(phi2+phiS)))/DSQRT(MD33**4-2.d0*MD33**2*
     . MQ33**2+MQ33**4+4.d0*Abtm**2*v1**2*Ybtm**2+4.d0*lambda**2*v2**2*
     . vevS**2*Ybtm**2-8.d0*Abtm*lambda*v1*v2*vevS*Ybtm**2*DCOS(phi2+
     . phiS)))**2+8.d0*Ybtm**4*(-1.5d0+DLOG((MD33**2+MQ33**2+2.d0*v1**2
     . *Ybtm**2+DSQRT(MD33**4-2.d0*MD33**2*MQ33**2+MQ33**4+4.d0*Abtm**2
     . *v1**2*Ybtm**2+4.d0*lambda**2*v2**2*vevS**2*Ybtm**2-8.d0*Abtm*
     . lambda*v1*v2*vevS*Ybtm**2*DCOS(phi2+phiS)))/(2.d0*Rsc**2)))*(v1+
     . (Abtm*(Abtm*v1-lambda*v2*vevS*DCOS(phi2+phiS)))/DSQRT(MD33**4-
     . 2.d0*MD33**2*MQ33**2+MQ33**4+4.d0*Abtm**2*v1**2*Ybtm**2+4.d0*
     . lambda**2*v2**2*vevS**2*Ybtm**2-8.d0*Abtm*lambda*v1*v2*vevS*
     . Ybtm**2*DCOS(phi2+phiS)))**2+12.d0*Ybtm**4*(v1+(Abtm*(-(Abtm*v1)
     . +lambda*v2*vevS*DCOS(phi2+phiS)))/DSQRT(MD33**4-2.d0*MD33**2*
     . MQ33**2+MQ33**4+4.d0*Abtm**2*v1**2*Ybtm**2+4.d0*lambda**2*v2**2*
     . vevS**2*Ybtm**2-8.d0*Abtm*lambda*v1*v2*vevS*Ybtm**2*DCOS(phi2+
     . phiS)))**2+8.d0*Ybtm**4*(-1.5d0+DLOG((MD33**2+MQ33**2+2.d0*v1**2
     . *Ybtm**2-DSQRT(MD33**4-2.d0*MD33**2*MQ33**2+MQ33**4+4.d0*Abtm**2
     . *v1**2*Ybtm**2+4.d0*lambda**2*v2**2*vevS**2*Ybtm**2-8.d0*Abtm*
     . lambda*v1*v2*vevS*Ybtm**2*DCOS(phi2+phiS)))/(2.d0*Rsc**2)))*(v1+
     . (Abtm*(-(Abtm*v1)+lambda*v2*vevS*DCOS(phi2+phiS)))/DSQRT(MD33**4
     . -2.d0*MD33**2*MQ33**2+MQ33**4+4.d0*Abtm**2*v1**2*Ybtm**2+4.d0*
     . lambda**2*v2**2*vevS**2*Ybtm**2-8.d0*Abtm*lambda*v1*v2*vevS*
     . Ybtm**2*DCOS(phi2+phiS)))**2+Ybtm**2*(1.d0+(4.d0*Abtm**2*Ybtm**2
     . *(Abtm*v1-lambda*v2*vevS*DCOS(phi2+phiS))**2)/(MD33**4-2.d0*
     . MD33**2*MQ33**2+MQ33**4+4.d0*Abtm**2*v1**2*Ybtm**2+4.d0*
     . lambda**2*v2**2*vevS**2*Ybtm**2-8.d0*Abtm*lambda*v1*v2*vevS*
     . Ybtm**2*DCOS(phi2+phiS))**1.5-Abtm**2/DSQRT(MD33**4-2.d0*
     . MD33**2*MQ33**2+MQ33**4+4.d0*Abtm**2*v1**2*Ybtm**2+4.d0*
     . lambda**2*v2**2*vevS**2*Ybtm**2-8.d0*Abtm*lambda*v1*v2*vevS*
     . Ybtm**2*DCOS(phi2+phiS)))*(MD33**2+MQ33**2+2.d0*v1**2*Ybtm**2-
     . DSQRT(MD33**4-2.d0*MD33**2*MQ33**2+MQ33**4+4.d0*Abtm**2*v1**2*
     . Ybtm**2+4.d0*lambda**2*v2**2*vevS**2*Ybtm**2-8.d0*Abtm*lambda*v1
     . *v2*vevS*Ybtm**2*DCOS(phi2+phiS)))+2.d0*Ybtm**2*(-1.5d0+
     . DLOG((MD33**2+MQ33**2+2.d0*v1**2*Ybtm**2-DSQRT(MD33**4-2.d0*
     . MD33**2*MQ33**2+MQ33**4+4.d0*Abtm**2*v1**2*Ybtm**2+4.d0*
     . lambda**2*v2**2*vevS**2*Ybtm**2-8.d0*Abtm*lambda*v1*v2*vevS*
     . Ybtm**2*DCOS(phi2+phiS)))/(2.d0*Rsc**2)))*(1.d0+(4.d0*Abtm**2*
     . Ybtm**2*(Abtm*v1-lambda*v2*vevS*DCOS(phi2+phiS))**2)/(MD33**4-
     . 2.d0*MD33**2*MQ33**2+MQ33**4+4.d0*Abtm**2*v1**2*Ybtm**2+4.d0*
     . lambda**2*v2**2*vevS**2*Ybtm**2-8.d0*Abtm*lambda*v1*v2*vevS*
     . Ybtm**2*DCOS(phi2+phiS))**1.5-Abtm**2/DSQRT(MD33**4-2.d0*
     . MD33**2*MQ33**2+MQ33**4+4.d0*Abtm**2*v1**2*Ybtm**2+4.d0*
     . lambda**2*v2**2*vevS**2*Ybtm**2-8.d0*Abtm*lambda*v1*v2*vevS*
     . Ybtm**2*DCOS(phi2+phiS)))*(MD33**2+MQ33**2+2.d0*v1**2*Ybtm**2-
     . DSQRT(MD33**4-2.d0*MD33**2*MQ33**2+MQ33**4+4.d0*Abtm**2*v1**2*
     . Ybtm**2+4.d0*lambda**2*v2**2*vevS**2*Ybtm**2-8.d0*Abtm*lambda*v1
     . *v2*vevS*Ybtm**2*DCOS(phi2+phiS)))+Ybtm**2*(1.d0-(4.d0*Abtm**2*
     . Ybtm**2*(Abtm*v1-lambda*v2*vevS*DCOS(phi2+phiS))**2)/(MD33**4-
     . 2.d0*MD33**2*MQ33**2+MQ33**4+4.d0*Abtm**2*v1**2*Ybtm**2+4.d0*
     . lambda**2*v2**2*vevS**2*Ybtm**2-8.d0*Abtm*lambda*v1*v2*vevS*
     . Ybtm**2*DCOS(phi2+phiS))**1.5+Abtm**2/DSQRT(MD33**4-2.d0*
     . MD33**2*MQ33**2+MQ33**4+4.d0*Abtm**2*v1**2*Ybtm**2+4.d0*
     . lambda**2*v2**2*vevS**2*Ybtm**2-8.d0*Abtm*lambda*v1*v2*vevS*
     . Ybtm**2*DCOS(phi2+phiS)))*(MD33**2+MQ33**2+2.d0*v1**2*Ybtm**2+
     . DSQRT(MD33**4-2.d0*MD33**2*MQ33**2+MQ33**4+4.d0*Abtm**2*v1**2*
     . Ybtm**2+4.d0*lambda**2*v2**2*vevS**2*Ybtm**2-8.d0*Abtm*lambda*v1
     . *v2*vevS*Ybtm**2*DCOS(phi2+phiS)))+2.d0*Ybtm**2*(-1.5d0+
     . DLOG((MD33**2+MQ33**2+2.d0*v1**2*Ybtm**2+DSQRT(MD33**4-2.d0*
     . MD33**2*MQ33**2+MQ33**4+4.d0*Abtm**2*v1**2*Ybtm**2+4.d0*
     . lambda**2*v2**2*vevS**2*Ybtm**2-8.d0*Abtm*lambda*v1*v2*vevS*
     . Ybtm**2*DCOS(phi2+phiS)))/(2.d0*Rsc**2)))*(1.d0-(4.d0*Abtm**2*
     . Ybtm**2*(Abtm*v1-lambda*v2*vevS*DCOS(phi2+phiS))**2)/(MD33**4-
     . 2.d0*MD33**2*MQ33**2+MQ33**4+4.d0*Abtm**2*v1**2*Ybtm**2+4.d0*
     . lambda**2*v2**2*vevS**2*Ybtm**2-8.d0*Abtm*lambda*v1*v2*vevS*
     . Ybtm**2*DCOS(phi2+phiS))**1.5+Abtm**2/DSQRT(MD33**4-2.d0*
     . MD33**2*MQ33**2+MQ33**4+4.d0*Abtm**2*v1**2*Ybtm**2+4.d0*
     . lambda**2*v2**2*vevS**2*Ybtm**2-8.d0*Abtm*lambda*v1*v2*vevS*
     . Ybtm**2*DCOS(phi2+phiS)))*(MD33**2+MQ33**2+2.d0*v1**2*Ybtm**2+
     . DSQRT(MD33**4-2.d0*MD33**2*MQ33**2+MQ33**4+4.d0*Abtm**2*v1**2*
     . Ybtm**2+4.d0*lambda**2*v2**2*vevS**2*Ybtm**2-8.d0*Abtm*lambda*v1
     . *v2*vevS*Ybtm**2*DCOS(phi2+phiS)))+(4.d0*lambda**2*vevS**2*
     . Ytop**4*(lambda*v1*vevS-Atop*v2*DCOS(phi2+phiS))**2*(MQ33**2+
     . MU33**2+2.d0*v2**2*Ytop**2-DSQRT(MQ33**4-2.d0*MQ33**2*MU33**2+
     . MU33**4+4.d0*Atop**2*v2**2*Ytop**2+4.d0*lambda**2*v1**2*vevS**2*
     . Ytop**2-8.d0*Atop*lambda*v1*v2*vevS*Ytop**2*DCOS(phi2+phiS))))/
     . (MQ33**4-2.d0*MQ33**2*MU33**2+MU33**4+4.d0*Atop**2*v2**2*Ytop**2
     . +4.d0*lambda**2*v1**2*vevS**2*Ytop**2-8.d0*Atop*lambda*v1*v2*
     . vevS*Ytop**2*DCOS(phi2+phiS))**1.5+(4.d0*lambda**2*vevS**2*
     . Ytop**4*(lambda*v1*vevS-Atop*v2*DCOS(phi2+phiS))**2*(-3.d0+2.d0*
     . DLOG((MQ33**2+MU33**2+2.d0*v2**2*Ytop**2-DSQRT(MQ33**4-2.d0*
     . MQ33**2*MU33**2+MU33**4+4.d0*Atop**2*v2**2*Ytop**2+4.d0*
     . lambda**2*v1**2*vevS**2*Ytop**2-8.d0*Atop*lambda*v1*v2*vevS*
     . Ytop**2*DCOS(phi2+phiS)))/(2.d0*Rsc**2)))*(MQ33**2+MU33**2+2.d0*
     . v2**2*Ytop**2-DSQRT(MQ33**4-2.d0*MQ33**2*MU33**2+MU33**4+4.d0*
     . Atop**2*v2**2*Ytop**2+4.d0*lambda**2*v1**2*vevS**2*Ytop**2-8.d0*
     . Atop*lambda*v1*v2*vevS*Ytop**2*DCOS(phi2+phiS))))/(MQ33**4-2.d0*
     . MQ33**2*MU33**2+MU33**4+4.d0*Atop**2*v2**2*Ytop**2+4.d0*
     . lambda**2*v1**2*vevS**2*Ytop**2-8.d0*Atop*lambda*v1*v2*vevS*
     . Ytop**2*DCOS(phi2+phiS))**1.5-(lambda**2*vevS**2*Ytop**2*
     . (MQ33**2+MU33**2+2.d0*v2**2*Ytop**2-DSQRT(MQ33**4-2.d0*MQ33**2*
     . MU33**2+MU33**4+4.d0*Atop**2*v2**2*Ytop**2+4.d0*lambda**2*v1**2*
     . vevS**2*Ytop**2-8.d0*Atop*lambda*v1*v2*vevS*Ytop**2*DCOS(phi2+
     . phiS))))/DSQRT(MQ33**4-2.d0*MQ33**2*MU33**2+MU33**4+4.d0*Atop**2
     . *v2**2*Ytop**2+4.d0*lambda**2*v1**2*vevS**2*Ytop**2-8.d0*Atop*
     . lambda*v1*v2*vevS*Ytop**2*DCOS(phi2+phiS))-(lambda**2*vevS**2*
     . Ytop**2*(-3.d0+2.d0*DLOG((MQ33**2+MU33**2+2.d0*v2**2*Ytop**2-
     . DSQRT(MQ33**4-2.d0*MQ33**2*MU33**2+MU33**4+4.d0*Atop**2*v2**2*
     . Ytop**2+4.d0*lambda**2*v1**2*vevS**2*Ytop**2-8.d0*Atop*lambda*v1
     . *v2*vevS*Ytop**2*DCOS(phi2+phiS)))/(2.d0*Rsc**2)))*(MQ33**2+
     . MU33**2+2.d0*v2**2*Ytop**2-DSQRT(MQ33**4-2.d0*MQ33**2*MU33**2+
     . MU33**4+4.d0*Atop**2*v2**2*Ytop**2+4.d0*lambda**2*v1**2*vevS**2*
     . Ytop**2-8.d0*Atop*lambda*v1*v2*vevS*Ytop**2*DCOS(phi2+phiS))))/
     . DSQRT(MQ33**4-2.d0*MQ33**2*MU33**2+MU33**4+4.d0*Atop**2*v2**2*
     . Ytop**2+4.d0*lambda**2*v1**2*vevS**2*Ytop**2-8.d0*Atop*lambda*v1
     . *v2*vevS*Ytop**2*DCOS(phi2+phiS))-(4.d0*lambda**2*vevS**2*
     . Ytop**4*(lambda*v1*vevS-Atop*v2*DCOS(phi2+phiS))**2*(MQ33**2+
     . MU33**2+2.d0*v2**2*Ytop**2+DSQRT(MQ33**4-2.d0*MQ33**2*MU33**2+
     . MU33**4+4.d0*Atop**2*v2**2*Ytop**2+4.d0*lambda**2*v1**2*vevS**2*
     . Ytop**2-8.d0*Atop*lambda*v1*v2*vevS*Ytop**2*DCOS(phi2+phiS))))/
     . (MQ33**4-2.d0*MQ33**2*MU33**2+MU33**4+4.d0*Atop**2*v2**2*Ytop**2
     . +4.d0*lambda**2*v1**2*vevS**2*Ytop**2-8.d0*Atop*lambda*v1*v2*
     . vevS*Ytop**2*DCOS(phi2+phiS))**1.5-(4.d0*lambda**2*vevS**2*
     . Ytop**4*(lambda*v1*vevS-Atop*v2*DCOS(phi2+phiS))**2*(-3.d0+2.d0*
     . DLOG((MQ33**2+MU33**2+2.d0*v2**2*Ytop**2+DSQRT(MQ33**4-2.d0*
     . MQ33**2*MU33**2+MU33**4+4.d0*Atop**2*v2**2*Ytop**2+4.d0*
     . lambda**2*v1**2*vevS**2*Ytop**2-8.d0*Atop*lambda*v1*v2*vevS*
     . Ytop**2*DCOS(phi2+phiS)))/(2.d0*Rsc**2)))*(MQ33**2+MU33**2+2.d0*
     . v2**2*Ytop**2+DSQRT(MQ33**4-2.d0*MQ33**2*MU33**2+MU33**4+4.d0*
     . Atop**2*v2**2*Ytop**2+4.d0*lambda**2*v1**2*vevS**2*Ytop**2-8.d0*
     . Atop*lambda*v1*v2*vevS*Ytop**2*DCOS(phi2+phiS))))/(MQ33**4-2.d0*
     . MQ33**2*MU33**2+MU33**4+4.d0*Atop**2*v2**2*Ytop**2+4.d0*
     . lambda**2*v1**2*vevS**2*Ytop**2-8.d0*Atop*lambda*v1*v2*vevS*
     . Ytop**2*DCOS(phi2+phiS))**1.5+(lambda**2*vevS**2*Ytop**2*
     . (MQ33**2+MU33**2+2.d0*v2**2*Ytop**2+DSQRT(MQ33**4-2.d0*MQ33**2*
     . MU33**2+MU33**4+4.d0*Atop**2*v2**2*Ytop**2+4.d0*lambda**2*v1**2*
     . vevS**2*Ytop**2-8.d0*Atop*lambda*v1*v2*vevS*Ytop**2*DCOS(phi2+
     . phiS))))/DSQRT(MQ33**4-2.d0*MQ33**2*MU33**2+MU33**4+4.d0*Atop**2
     . *v2**2*Ytop**2+4.d0*lambda**2*v1**2*vevS**2*Ytop**2-8.d0*Atop*
     . lambda*v1*v2*vevS*Ytop**2*DCOS(phi2+phiS))+(lambda**2*vevS**2*
     . Ytop**2*(-3.d0+2.d0*DLOG((MQ33**2+MU33**2+2.d0*v2**2*Ytop**2+
     . DSQRT(MQ33**4-2.d0*MQ33**2*MU33**2+MU33**4+4.d0*Atop**2*v2**2*
     . Ytop**2+4.d0*lambda**2*v1**2*vevS**2*Ytop**2-8.d0*Atop*lambda*v1
     . *v2*vevS*Ytop**2*DCOS(phi2+phiS)))/(2.d0*Rsc**2)))*(MQ33**2+
     . MU33**2+2.d0*v2**2*Ytop**2+DSQRT(MQ33**4-2.d0*MQ33**2*MU33**2+
     . MU33**4+4.d0*Atop**2*v2**2*Ytop**2+4.d0*lambda**2*v1**2*vevS**2*
     . Ytop**2-8.d0*Atop*lambda*v1*v2*vevS*Ytop**2*DCOS(phi2+phiS))))/
     . DSQRT(MQ33**4-2.d0*MQ33**2*MU33**2+MU33**4+4.d0*Atop**2*v2**2*
     . Ytop**2+4.d0*lambda**2*v1**2*vevS**2*Ytop**2-8.d0*Atop*lambda*v1
     . *v2*vevS*Ytop**2*DCOS(phi2+phiS))))/(64.d0*Pi**2)