(* Results for TTP13-008 with general group factors for the QCD part *)
(* In QCD: *)
(*
 dR=3;
 ca=3;
 cf=4/3;
 tr=1/2;
 *)
(* Number of generations in the SM *)
(*
 NG=3;
 *)
(*
 h=1/(16 Pi^2);
 z3=Zeta[3];
 *)

(* Beta function for lambda *)

betalambdal1 =(

       + ytau^4*h * (
          - 1
          )

       + yb^4*h * (
          - dR
          )

       + g2^4*h * (
          + 9/16
          )

       + g1^2*g2^2*h * (
          + 3/8
          )

       + g1^4*h * (
          + 3/16
          )

       + lambda*ytau^2*h * (
          + 2
          )

       + lambda*yb^2*h * (
          + 2*dR
          )

       + lambda*g2^2*h * (
          - 9/2
          )

       + lambda*g1^2*h * (
          - 3/2
          )

       + lambda^2*h * (
          + 12
          )

       + yt^2*lambda*h * (
          + 2*dR
          )

       + yt^4*h * (
          - dR
		     ));


betalambdal2 =(

       + ytau^6*h^2 * (
          + 5
          )

       + yb^6*h^2 * (
          + 5*dR
          )

       + g2^4*ytau^2*h^2 * (
          - 3/8
          )

       + g2^4*yb^2*h^2 * (
          - 3/8*dR
          )

       + g2^6*h^2 * (
          + 497/32
          - 1/2*NG
          - 1/2*NG*dR
          )

       + g1^2*ytau^4*h^2 * (
          - 2
          )

       + g1^2*yb^4*h^2 * (
          + 2/9*dR
          )

       + g1^2*g2^2*ytau^2*h^2 * (
          + 11/4
          )

       + g1^2*g2^2*yb^2*h^2 * (
          + 3/4*dR
          )

       + g1^2*g2^4*h^2 * (
          - 97/96
          - 1/6*NG
          - 1/6*NG*dR
          )

       + g1^4*ytau^2*h^2 * (
          - 25/8
          )

       + g1^4*yb^2*h^2 * (
          + 5/24*dR
          )

       + g1^4*g2^2*h^2 * (
          - 239/96
          - 1/2*NG
          - 11/54*NG*dR
          )

       + g1^6*h^2 * (
          - 59/96
          - 1/2*NG
          - 11/54*NG*dR
          )

       + lambda*ytau^4*h^2 * (
          - 1/2
          )

       + lambda*yb^4*h^2 * (
          - 1/2*dR
          )

       + lambda*g2^2*ytau^2*h^2 * (
          + 15/4
          )

       + lambda*g2^2*yb^2*h^2 * (
          + 15/4*dR
          )

       + lambda*g2^4*h^2 * (
          - 313/16
          + 5/4*NG
          + 5/4*NG*dR
          )

       + lambda*g1^2*ytau^2*h^2 * (
          + 25/4
          )

       + lambda*g1^2*yb^2*h^2 * (
          + 25/36*dR
          )

       + lambda*g1^2*g2^2*h^2 * (
          + 39/8
          )

       + lambda*g1^4*h^2 * (
          + 229/48
          + 5/4*NG
          + 55/108*NG*dR
          )

       + lambda^2*ytau^2*h^2 * (
          - 24
          )

       + lambda^2*yb^2*h^2 * (
          - 24*dR
          )

       + lambda^2*g2^2*h^2 * (
          + 54
          )

       + lambda^2*g1^2*h^2 * (
          + 18
          )

       + lambda^3*h^2 * (
          - 156
          )

       + yt^2*yb^4*h^2 * (
          - dR
          )

       + yt^2*g2^4*h^2 * (
          - 3/8*dR
          )

       + yt^2*g1^2*g2^2*h^2 * (
          + 7/4*dR
          )

       + yt^2*g1^4*h^2 * (
          - 19/24*dR
          )

       + yt^2*lambda*yb^2*h^2 * (
          - 7*dR
          )

       + yt^2*lambda*g2^2*h^2 * (
          + 15/4*dR
          )

       + yt^2*lambda*g1^2*h^2 * (
          + 85/36*dR
          )

       + yt^2*lambda^2*h^2 * (
          - 24*dR
          )

       + yt^4*yb^2*h^2 * (
          - dR
          )

       + yt^4*g1^2*h^2 * (
          - 4/9*dR
          )

       + yt^4*lambda*h^2 * (
          - 1/2*dR
          )

       + yt^6*h^2 * (
          + 5*dR
          )

       + gs^2*yb^4*h^2 * (
          - 4*cf*dR
          )

       + gs^2*lambda*yb^2*h^2 * (
          + 10*cf*dR
          )

       + gs^2*yt^2*lambda*h^2 * (
          + 10*cf*dR
          )

       + gs^2*yt^4*h^2 * (
          - 4*cf*dR
			    ));


betalambdal3 =(

       + ytau^8*h^3 * (
          - 143/8
          - 12*z3
          )

       + yb^2*ytau^6*h^3 * (
          - 99/8*dR
          )

       + yb^4*ytau^4*h^3 * (
          - 24*dR
          )

       + yb^6*ytau^2*h^3 * (
          - 99/8*dR
          )

       + yb^8*h^3 * (
          + 13/2*dR
          - 195/8*dR^2
          - 12*z3*dR
          )

       + g2^2*ytau^6*h^3 * (
          + 1137/32
          - 9*z3
          )

       + g2^2*yb^6*h^3 * (
          + 1137/32*dR
          - 9*z3*dR
          )

       + g2^4*ytau^4*h^3 * (
          + 4503/128
          - 273/16*z3
          - 13/16*NG
          - 13/16*NG*dR
          )

       + g2^4*yb^2*ytau^2*h^3 * (
          + 3/8*dR
          )

       + g2^4*yb^4*h^3 * (
          + 4479/128*dR
          + 3/16*dR^2
          - 273/16*z3*dR
          - 13/16*NG*dR
          - 13/16*NG*dR^2
          )

       + g2^6*ytau^2*h^3 * (
          - 5739/256
          + 99/4*z3
          + 9/8*NG
          + 9/8*NG*dR
          )

       + g2^6*yb^2*h^3 * (
          - 5739/256*dR
          + 99/4*z3*dR
          + 9/8*NG*dR
          + 9/8*NG*dR^2
          )

       + g2^8*h^3 * (
          + 982291/3072
          - 2781/128*z3
          - 14749/768*NG
          - 14749/768*NG*dR
          - 45/4*NG*z3
          - 45/4*NG*z3*dR
          - 5/48*NG^2
          - 5/24*NG^2*dR
          - 5/48*NG^2*dR^2
          )

       + g1^2*ytau^6*h^3 * (
          + 135/32
          + 33*z3
          )

       + g1^2*yb^6*h^3 * (
          + 5111/288*dR
          - 25/3*z3*dR
          )

       + g1^2*g2^2*ytau^4*h^3 * (
          - 15/64
          - 381/8*z3
          )

       + g1^2*g2^2*yb^2*ytau^2*h^3 * (
          - 5/12*dR
          )

       + g1^2*g2^2*yb^4*h^3 * (
          - 2591/576*dR
          - 3/8*dR^2
          - 311/24*z3*dR
          )

       + g1^2*g2^4*ytau^2*h^3 * (
          + 1833/256
          - 3/2*z3
          - 1/8*NG
          - 1/8*NG*dR
          )

       + g1^2*g2^4*yb^2*h^3 * (
          + 1393/256*dR
          + 3*z3*dR
          + 5/24*NG*dR
          + 5/24*NG*dR^2
          )

       + g1^2*g2^6*h^3 * (
          - 54053/3456
          - 405/32*z3
          - 11095/3456*NG
          - 7423/3456*NG*dR
          + 3/4*NG*z3
          - 1/4*NG*z3*dR
          - 5/216*NG^2
          - 5/108*NG^2*dR
          - 5/216*NG^2*dR^2
          )

       + g1^4*ytau^4*h^3 * (
          + 5697/128
          + 375/16*z3
          + 39/16*NG
          + 143/144*NG*dR
          )

       + g1^4*yb^2*ytau^2*h^3 * (
          + 41/72*dR
          )

       + g1^4*yb^4*h^3 * (
          + 20681/10368*dR
          - 77/432*dR^2
          - 2035/432*z3*dR
          - 83/48*NG*dR
          - 913/1296*NG*dR^2
          )

       + g1^4*g2^2*ytau^2*h^3 * (
          + 6657/256
          - 15/2*z3
          - 3/8*NG
          - 11/72*NG*dR
          )

       + g1^4*g2^2*yb^2*h^3 * (
          + 4403/768*dR
          + 3/2*z3*dR
          + 5/8*NG*dR
          + 55/216*NG*dR^2
          )

       + g1^4*g2^4*h^3 * (
          - 64693/3456
          + 873/64*z3
          - 5129/1728*NG
          + 16579/15552*NG*dR
          + 21/4*NG*z3
          + 7/12*NG*z3*dR
          - 5/72*NG^2
          - 95/972*NG^2*dR
          - 55/1944*NG^2*dR^2
          )

       + g1^6*ytau^2*h^3 * (
          + 3929/256
          - 15/4*z3
          + 33/8*NG
          + 121/72*NG*dR
          )

       + g1^6*yb^2*h^3 * (
          + 12043/6912*dR
          + 5/12*z3*dR
          + 19/24*NG*dR
          + 209/648*NG*dR^2
          )

       + g1^6*g2^2*h^3 * (
          - 29779/6912
          + 75/32*z3
          - 5375/1152*NG
          - 23629/31104*NG*dR
          + 9/2*NG*z3
          + 41/54*NG*z3*dR
          - 5/24*NG^2
          - 55/324*NG^2*dR
          - 605/17496*NG^2*dR^2
          )

       + g1^8*h^3 * (
          - 6845/9216
          + 99/128*z3
          - 5987/768*NG
          - 29057/20736*NG*dR
          + 27/4*NG*z3
          + 137/108*NG*z3*dR
          - 5/16*NG^2
          - 55/216*NG^2*dR
          - 605/11664*NG^2*dR^2
          )

       + lambda*ytau^6*h^3 * (
          - 1241/8
          - 66*z3
          )

       + lambda*yb^2*ytau^4*h^3 * (
          + 80*dR
          )

       + lambda*yb^4*ytau^2*h^3 * (
          + 80*dR
          )

       + lambda*yb^6*h^3 * (
          - 1881/8*dR
          + 80*dR^2
          - 66*z3*dR
          )

       + lambda*g2^2*ytau^4*h^3 * (
          - 1587/8
          + 171*z3
          )

       + lambda*g2^2*yb^2*ytau^2*h^3 * (
          - 9*dR
          )

       + lambda*g2^2*yb^4*h^3 * (
          - 1551/8*dR
          - 9/2*dR^2
          + 171*z3*dR
          )

       + lambda*g2^4*ytau^2*h^3 * (
          - 1311/64
          - 117/2*z3
          - 21/16*NG
          - 21/16*NG*dR
          )

       + lambda*g2^4*yb^2*h^3 * (
          - 1311/64*dR
          - 117/2*z3*dR
          - 21/16*NG*dR
          - 21/16*NG*dR^2
          )

       + lambda*g2^6*h^3 * (
          - 46489/288
          + 2259/8*z3
          + 3515/144*NG
          + 3515/144*NG*dR
          + 45/2*NG*z3
          + 45/2*NG*z3*dR
          + 35/72*NG^2
          + 35/36*NG^2*dR
          + 35/72*NG^2*dR^2
          )

       + lambda*g1^2*ytau^4*h^3 * (
          + 507/8
          - 117*z3
          )

       + lambda*g1^2*yb^2*ytau^2*h^3 * (
          - 3*dR
          )

       + lambda*g1^2*yb^4*h^3 * (
          - 5413/72*dR
          - 3/2*dR^2
          + 83*z3*dR
          )

       + lambda*g1^2*g2^2*ytau^2*h^3 * (
          - 3771/32
          + 126*z3
          )

       + lambda*g1^2*g2^2*yb^2*h^3 * (
          - 1003/32*dR
          + 4*z3*dR
          )

       + lambda*g1^2*g2^4*h^3 * (
          + 4553/32
          - 249/8*z3
          + 111/16*NG
          + 51/16*NG*dR
          - 9/2*NG*z3
          - 1/2*NG*z3*dR
          )

       + lambda*g1^4*ytau^2*h^3 * (
          - 1783/64
          - 123/2*z3
          - 117/16*NG
          - 143/48*NG*dR
          )

       + lambda*g1^4*yb^2*h^3 * (
          - 127303/5184*dR
          - 47/18*z3*dR
          - 31/48*NG*dR
          - 341/1296*NG*dR^2
          )

       + lambda*g1^4*g2^2*h^3 * (
          + 979/8
          - 3/8*z3
          + 99/8*NG
          + 91/24*NG*dR
          - 9/2*NG*z3
          - 1/2*NG*z3*dR
          )

       + lambda*g1^6*h^3 * (
          + 12679/432
          + 9/8*z3
          + 1469/72*NG
          + 10759/1944*NG*dR
          - 27/2*NG*z3
          - 137/54*NG*z3*dR
          + 35/24*NG^2
          + 385/324*NG^2*dR
          + 4235/17496*NG^2*dR^2
          )

       + lambda^2*ytau^4*h^3 * (
          + 717/2
          + 252*z3
          )

       + lambda^2*yb^2*ytau^2*h^3 * (
          - 72*dR
          )

       + lambda^2*yb^4*h^3 * (
          + 789/2*dR
          - 36*dR^2
          + 252*z3*dR
          )

       + lambda^2*g2^2*ytau^2*h^3 * (
          + 213/4
          - 144*z3
          )

       + lambda^2*g2^2*yb^2*h^3 * (
          + 213/4*dR
          - 144*z3*dR
          )

       + lambda^2*g2^4*h^3 * (
          + 1995/8
          - 513*z3
          - 141/4*NG
          - 141/4*NG*dR
          )

       + lambda^2*g1^2*ytau^2*h^3 * (
          - 541/4
          + 96*z3
          )

       + lambda^2*g1^2*yb^2*h^3 * (
          + 139/4*dR
          - 64*z3*dR
          )

       + lambda^2*g1^2*g2^2*h^3 * (
          - 333
          - 162*z3
          )

       + lambda^2*g1^4*h^3 * (
          - 183
          - 81*z3
          - 141/4*NG
          - 517/36*NG*dR
          )

       + lambda^3*ytau^2*h^3 * (
          + 291
          )

       + lambda^3*yb^2*h^3 * (
          + 291*dR
          )

       + lambda^3*g2^2*h^3 * (
          - 474
          + 72*z3
          )

       + lambda^3*g1^2*h^3 * (
          - 158
          + 24*z3
          )

       + lambda^4*h^3 * (
          + 3588
          + 2016*z3
          )

       + yt^2*ytau^6*h^3 * (
          - 99/8*dR
          )

       + yt^2*yb^2*ytau^4*h^3 * (
          + 4*dR
          )

       + yt^2*yb^4*ytau^2*h^3 * (
          + 15/8*dR
          )

       + yt^2*yb^6*h^3 * (
          - 83/8*dR
          - 13/2*dR^2
          - 12*z3*dR
          )

       + yt^2*g2^2*yb^4*h^3 * (
          + 159/32*dR
          )

       + yt^2*g2^4*ytau^2*h^3 * (
          + 3/8*dR
          )

       + yt^2*g2^4*yb^2*h^3 * (
          - 189/64*dR
          + 3/8*dR^2
          + 39/2*z3*dR
          - NG*dR
          - NG*dR^2
          )

       + yt^2*g2^6*h^3 * (
          - 5739/256*dR
          + 99/4*z3*dR
          + 9/8*NG*dR
          + 9/8*NG*dR^2
          )

       + yt^2*g1^2*yb^4*h^3 * (
          - 2299/288*dR
          + 26/3*z3*dR
          )

       + yt^2*g1^2*g2^2*ytau^2*h^3 * (
          + 29/12*dR
          )

       + yt^2*g1^2*g2^2*yb^2*h^3 * (
          + 1505/288*dR
          - 7/12*dR^2
          + 31/6*z3*dR
          )

       + yt^2*g1^2*g2^4*h^3 * (
          + 3103/768*dR
          + 9/4*z3*dR
          + 1/24*NG*dR
          + 1/24*NG*dR^2
          )

       + yt^2*g1^4*ytau^2*h^3 * (
          + 701/72*dR
          )

       + yt^2*g1^4*yb^2*h^3 * (
          - 1991/576*dR
          - 17/216*dR^2
          - 1/3*z3*dR
          )

       + yt^2*g1^4*g2^2*h^3 * (
          + 23521/2304*dR
          - z3*dR
          + 1/8*NG*dR
          + 11/216*NG*dR^2
          )

       + yt^2*g1^6*h^3 * (
          + 42943/6912*dR
          - 5/6*z3*dR
          + 43/24*NG*dR
          + 473/648*NG*dR^2
          )

       + yt^2*lambda*ytau^4*h^3 * (
          + 80*dR
          )

       + yt^2*lambda*yb^2*ytau^2*h^3 * (
          + 7*dR
          )

       + yt^2*lambda*yb^4*h^3 * (
          + 45/8*dR
          + 87*dR^2
          + 48*z3*dR
          )

       + yt^2*lambda*g2^2*ytau^2*h^3 * (
          - 9*dR
          )

       + yt^2*lambda*g2^2*yb^2*h^3 * (
          - 69/4*dR
          - 9*dR^2
          + 18*z3*dR
          )

       + yt^2*lambda*g2^4*h^3 * (
          - 1311/64*dR
          - 117/2*z3*dR
          - 21/16*NG*dR
          - 21/16*NG*dR^2
          )

       + yt^2*lambda*g1^2*ytau^2*h^3 * (
          - 3*dR
          )

       + yt^2*lambda*g1^2*yb^2*h^3 * (
          - 605/36*dR
          - 3*dR^2
          - 2/3*z3*dR
          )

       + yt^2*lambda*g1^2*g2^2*h^3 * (
          - 6509/96*dR
          + 59*z3*dR
          )

       + yt^2*lambda*g1^4*h^3 * (
          - 112447/5184*dR
          - 449/18*z3*dR
          - 127/48*NG*dR
          - 1397/1296*NG*dR^2
          )

       + yt^2*lambda^2*ytau^2*h^3 * (
          - 72*dR
          )

       + yt^2*lambda^2*yb^2*h^3 * (
          + 255*dR
          - 72*dR^2
          - 288*z3*dR
          )

       + yt^2*lambda^2*g2^2*h^3 * (
          + 213/4*dR
          - 144*z3*dR
          )

       + yt^2*lambda^2*g1^2*h^3 * (
          - 65/4*dR
          - 16*z3*dR
          )

       + yt^2*lambda^3*h^3 * (
          + 291*dR
          )

       + yt^4*ytau^4*h^3 * (
          - 24*dR
          )

       + yt^4*yb^2*ytau^2*h^3 * (
          + 15/8*dR
          )

       + yt^4*yb^4*h^3 * (
          + 147/4*dR
          - 49/4*dR^2
          + 24*z3*dR
          )

       + yt^4*g2^2*yb^2*h^3 * (
          + 159/32*dR
          )

       + yt^4*g2^4*h^3 * (
          + 4479/128*dR
          + 3/16*dR^2
          - 273/16*z3*dR
          - 13/16*NG*dR
          - 13/16*NG*dR^2
          )

       + yt^4*g1^2*yb^2*h^3 * (
          + 1337/288*dR
          - 28/3*z3*dR
          )

       + yt^4*g1^2*g2^2*h^3 * (
          - 1871/576*dR
          + 11/24*dR^2
          - 743/24*z3*dR
          )

       + yt^4*g1^4*h^3 * (
          + 56345/10368*dR
          + 619/432*dR^2
          + 2957/432*z3*dR
          - 23/48*NG*dR
          - 253/1296*NG*dR^2
          )

       + yt^4*lambda*ytau^2*h^3 * (
          + 80*dR
          )

       + yt^4*lambda*yb^2*h^3 * (
          + 45/8*dR
          + 87*dR^2
          + 48*z3*dR
          )

       + yt^4*lambda*g2^2*h^3 * (
          - 1551/8*dR
          - 9/2*dR^2
          + 171*z3*dR
          )

       + yt^4*lambda*g1^2*h^3 * (
          - 2161/72*dR
          - 3/2*dR^2
          + 19*z3*dR
          )

       + yt^4*lambda^2*h^3 * (
          + 789/2*dR
          - 36*dR^2
          + 252*z3*dR
          )

       + yt^6*ytau^2*h^3 * (
          - 99/8*dR
          )

       + yt^6*yb^2*h^3 * (
          - 83/8*dR
          - 13/2*dR^2
          - 12*z3*dR
          )

       + yt^6*g2^2*h^3 * (
          + 1137/32*dR
          - 9*z3*dR
          )

       + yt^6*g1^2*h^3 * (
          + 3467/288*dR
          + 17/3*z3*dR
          )

       + yt^6*lambda*h^3 * (
          - 1881/8*dR
          + 80*dR^2
          - 66*z3*dR
          )

       + yt^8*h^3 * (
          + 13/2*dR
          - 195/8*dR^2
          - 12*z3*dR
          )

       + gs^2*yb^6*h^3 * (
          - 19/2*cf*dR
          + 60*z3*cf*dR
          )

       + gs^2*g2^2*yb^4*h^3 * (
          - 31/8*cf*dR
          + 6*z3*cf*dR
          )

       + gs^2*g2^4*yb^2*h^3 * (
          + 651/32*cf*dR
          - 27/2*z3*cf*dR
          )

       + gs^2*g2^6*h^3 * (
          - 153/32*NG*cf*dR
          + 9/2*NG*z3*cf*dR
          )

       + gs^2*g1^2*yb^4*h^3 * (
          - 641/72*cf*dR
          + 34/3*z3*cf*dR
          )

       + gs^2*g1^2*g2^2*yb^2*h^3 * (
          + 233/16*cf*dR
          - 9*z3*cf*dR
          )

       + gs^2*g1^2*g2^4*h^3 * (
          - 51/32*NG*cf*dR
          + 3/2*NG*z3*cf*dR
          )

       + gs^2*g1^4*yb^2*h^3 * (
          + 683/96*cf*dR
          - 9/2*z3*cf*dR
          )

       + gs^2*g1^4*g2^2*h^3 * (
          - 187/96*NG*cf*dR
          + 11/6*NG*z3*cf*dR
          )

       + gs^2*g1^6*h^3 * (
          - 187/96*NG*cf*dR
          + 11/6*NG*z3*cf*dR
          )

       + gs^2*lambda*yb^4*h^3 * (
          + 895/4*cf*dR
          - 324*z3*cf*dR
          )

       + gs^2*lambda*g2^2*yb^2*h^3 * (
          - 489/8*cf*dR
          + 54*z3*cf*dR
          )

       + gs^2*lambda*g2^4*h^3 * (
          + 135/8*NG*cf*dR
          - 18*NG*z3*cf*dR
          )

       + gs^2*lambda*g1^2*yb^2*h^3 * (
          - 991/72*cf*dR
          + 10*z3*cf*dR
          )

       + gs^2*lambda*g1^4*h^3 * (
          + 55/8*NG*cf*dR
          - 22/3*NG*z3*cf*dR
          )

       + gs^2*lambda^2*yb^2*h^3 * (
          - 306*cf*dR
          + 288*z3*cf*dR
          )

       + gs^2*yt^2*yb^4*h^3 * (
          - 1/2*cf*dR
          - 12*z3*cf*dR
          )

       + gs^2*yt^2*g2^2*yb^2*h^3 * (
          - 2*cf*dR
          + 24*z3*cf*dR
          )

       + gs^2*yt^2*g2^4*h^3 * (
          + 651/32*cf*dR
          - 27/2*z3*cf*dR
          )

       + gs^2*yt^2*g1^2*g2^2*h^3 * (
          + 249/16*cf*dR
          - 9*z3*cf*dR
          )

       + gs^2*yt^2*g1^4*h^3 * (
          + 587/96*cf*dR
          - 9/2*z3*cf*dR
          )

       + gs^2*yt^2*lambda*yb^2*h^3 * (
          + 41/2*cf*dR
          - 24*z3*cf*dR
          )

       + gs^2*yt^2*lambda*g2^2*h^3 * (
          - 489/8*cf*dR
          + 54*z3*cf*dR
          )

       + gs^2*yt^2*lambda*g1^2*h^3 * (
          - 2419/72*cf*dR
          + 34*z3*cf*dR
          )

       + gs^2*yt^2*lambda^2*h^3 * (
          - 306*cf*dR
          + 288*z3*cf*dR
          )

       + gs^2*yt^4*yb^2*h^3 * (
          - 1/2*cf*dR
          - 12*z3*cf*dR
          )

       + gs^2*yt^4*g2^2*h^3 * (
          - 31/8*cf*dR
          + 6*z3*cf*dR
          )

       + gs^2*yt^4*g1^2*h^3 * (
          + 931/72*cf*dR
          - 14/3*z3*cf*dR
          )

       + gs^2*yt^4*lambda*h^3 * (
          + 895/4*cf*dR
          - 324*z3*cf*dR
          )

       + gs^2*yt^6*h^3 * (
          - 19/2*cf*dR
          + 60*z3*cf*dR
          )

       + gs^4*yb^4*h^3 * (
          + 131/2*cf^2*dR
          + 48*tr*cf*dR
          - 109/2*ca*cf*dR
          - 48*z3*cf^2*dR
          + 24*z3*ca*cf*dR
          + 20*NG*tr*cf*dR
          )

       + gs^4*lambda*yb^2*h^3 * (
          - 119/2*cf^2*dR
          + 77*ca*cf*dR
          + 72*z3*cf^2*dR
          - 36*z3*ca*cf*dR
          - 32*NG*tr*cf*dR
          )

       + gs^4*yt^2*yb^2*h^3 * (
          + 96*tr*cf*dR
          )

       + gs^4*yt^2*lambda*h^3 * (
          - 119/2*cf^2*dR
          + 77*ca*cf*dR
          + 72*z3*cf^2*dR
          - 36*z3*ca*cf*dR
          - 32*NG*tr*cf*dR
          )

       + gs^4*yt^4*h^3 * (
          + 131/2*cf^2*dR
          + 48*tr*cf*dR
          - 109/2*ca*cf*dR
          - 48*z3*cf^2*dR
          + 24*z3*ca*cf*dR
          + 20*NG*tr*cf*dR
			    ));

(* Beta function for m^2 *)
betam2l1 =(

       + ytau^2*h * (
          + 1
          )

       + yb^2*h * (
          + dR
          )

       + g2^2*h * (
          - 9/4
          )

       + g1^2*h * (
          - 3/4
          )

       + lambda*h * (
          + 6
          )

       + yt^2*h * (
          + dR
		     ));


betam2l2 =(

       + ytau^4*h^2 * (
          - 9/4
          )

       + yb^4*h^2 * (
          - 9/4*dR
          )

       + g2^2*ytau^2*h^2 * (
          + 15/8
          )

       + g2^2*yb^2*h^2 * (
          + 15/8*dR
          )

       + g2^4*h^2 * (
          - 385/32
          + 5/8*NG
          + 5/8*NG*dR
          )

       + g1^2*ytau^2*h^2 * (
          + 25/8
          )

       + g1^2*yb^2*h^2 * (
          + 25/72*dR
          )

       + g1^2*g2^2*h^2 * (
          + 15/16
          )

       + g1^4*h^2 * (
          + 157/96
          + 5/8*NG
          + 55/216*NG*dR
          )

       + lambda*ytau^2*h^2 * (
          - 12
          )

       + lambda*yb^2*h^2 * (
          - 12*dR
          )

       + lambda*g2^2*h^2 * (
          + 36
          )

       + lambda*g1^2*h^2 * (
          + 12
          )

       + lambda^2*h^2 * (
          - 30
          )

       + yt^2*yb^2*h^2 * (
          - 7/2*dR
          )

       + yt^2*g2^2*h^2 * (
          + 15/8*dR
          )

       + yt^2*g1^2*h^2 * (
          + 85/72*dR
          )

       + yt^2*lambda*h^2 * (
          - 12*dR
          )

       + yt^4*h^2 * (
          - 9/4*dR
          )

       + gs^2*yb^2*h^2 * (
          + 5*cf*dR
          )

       + gs^2*yt^2*h^2 * (
          + 5*cf*dR
			    ));


betam2l3 =(

       + ytau^6*h^3 * (
          - 233/16
          + 15*z3
          )

       + yb^2*ytau^4*h^3 * (
          + 24*dR
          )

       + yb^4*ytau^2*h^3 * (
          + 24*dR
          )

       + yb^6*h^3 * (
          - 617/16*dR
          + 24*dR^2
          + 15*z3*dR
          )

       + g2^2*ytau^4*h^3 * (
          - 987/16
          + 54*z3
          )

       + g2^2*yb^2*ytau^2*h^3 * (
          - 9/2*dR
          )

       + g2^2*yb^4*h^3 * (
          - 951/16*dR
          - 9/4*dR^2
          + 54*z3*dR
          )

       + g2^4*ytau^2*h^3 * (
          - 255/128
          - 81/4*z3
          - 21/32*NG
          - 21/32*NG*dR
          )

       + g2^4*yb^2*h^3 * (
          - 255/128*dR
          - 81/4*z3*dR
          - 21/32*NG*dR
          - 21/32*NG*dR^2
          )

       + g2^6*h^3 * (
          - 39415/576
          + 711/16*z3
          + 2867/288*NG
          + 2867/288*NG*dR
          + 45/4*NG*z3
          + 45/4*NG*z3*dR
          + 35/144*NG^2
          + 35/72*NG^2*dR
          + 35/144*NG^2*dR^2
          )

       + g1^2*ytau^4*h^3 * (
          + 291/16
          - 36*z3
          )

       + g1^2*yb^2*ytau^2*h^3 * (
          - 3/2*dR
          )

       + g1^2*yb^4*h^3 * (
          - 959/48*dR
          - 3/4*dR^2
          + 24*z3*dR
          )

       + g1^2*g2^2*ytau^2*h^3 * (
          - 2331/64
          + 45*z3
          )

       + g1^2*g2^2*yb^2*h^3 * (
          - 865/192*dR
          )

       + g1^2*g2^4*h^3 * (
          + 2691/64
          - 405/16*z3
          + 87/32*NG
          + 27/32*NG*dR
          - 9/4*NG*z3
          - 1/4*NG*z3*dR
          )

       + g1^4*ytau^2*h^3 * (
          - 3607/128
          - 15/4*z3
          - 117/32*NG
          - 143/96*NG*dR
          )

       + g1^4*yb^2*h^3 * (
          - 79207/10368*dR
          - 35/36*z3*dR
          - 31/96*NG*dR
          - 341/2592*NG*dR^2
          )

       + g1^4*g2^2*h^3 * (
          + 1053/32
          - 207/16*z3
          + 63/16*NG
          + 47/48*NG*dR
          - 9/4*NG*z3
          - 1/4*NG*z3*dR
          )

       + g1^6*h^3 * (
          + 839/108
          - 51/16*z3
          + 1145/144*NG
          + 7195/3888*NG*dR
          - 27/4*NG*z3
          - 137/108*NG*z3*dR
          + 35/48*NG^2
          + 385/648*NG^2*dR
          + 4235/34992*NG^2*dR^2
          )

       + lambda*ytau^4*h^3 * (
          + 261/4
          + 72*z3
          )

       + lambda*yb^2*ytau^2*h^3 * (
          - 36*dR
          )

       + lambda*yb^4*h^3 * (
          + 333/4*dR
          - 18*dR^2
          + 72*z3*dR
          )

       + lambda*g2^2*ytau^2*h^3 * (
          + 189/8
          - 108*z3
          )

       + lambda*g2^2*yb^2*h^3 * (
          + 189/8*dR
          - 108*z3*dR
          )

       + lambda*g2^4*h^3 * (
          + 11511/32
          - 162*z3
          - 153/8*NG
          - 153/8*NG*dR
          )

       + lambda*g1^2*ytau^2*h^3 * (
          - 549/8
          + 36*z3
          )

       + lambda*g1^2*yb^2*h^3 * (
          + 131/8*dR
          - 44*z3*dR
          )

       + lambda*g1^2*g2^2*h^3 * (
          - 1701/16
          + 36*z3
          )

       + lambda*g1^4*h^3 * (
          - 1077/32
          - 18*z3
          - 153/8*NG
          - 187/24*NG*dR
          )

       + lambda^2*ytau^2*h^3 * (
          + 99/2
          )

       + lambda^2*yb^2*h^3 * (
          + 99/2*dR
          )

       + lambda^2*g2^2*h^3 * (
          - 63
          - 108*z3
          )

       + lambda^2*g1^2*h^3 * (
          - 21
          - 36*z3
          )

       + lambda^3*h^3 * (
          + 1026
          )

       + yt^2*ytau^4*h^3 * (
          + 24*dR
          )

       + yt^2*yb^2*ytau^2*h^3 * (
          + 7/2*dR
          )

       + yt^2*yb^4*h^3 * (
          + 29/16*dR
          + 55/2*dR^2
          + 12*z3*dR
          )

       + yt^2*g2^2*ytau^2*h^3 * (
          - 9/2*dR
          )

       + yt^2*g2^2*yb^2*h^3 * (
          + 27/8*dR
          - 9/2*dR^2
          - 9*z3*dR
          )

       + yt^2*g2^4*h^3 * (
          - 255/128*dR
          - 81/4*z3*dR
          - 21/32*NG*dR
          - 21/32*NG*dR^2
          )

       + yt^2*g1^2*ytau^2*h^3 * (
          - 3/2*dR
          )

       + yt^2*g1^2*yb^2*h^3 * (
          - 605/72*dR
          - 3/2*dR^2
          - 1/3*z3*dR
          )

       + yt^2*g1^2*g2^2*h^3 * (
          - 3277/192*dR
          + 39/2*z3*dR
          )

       + yt^2*g1^4*h^3 * (
          - 123103/10368*dR
          - 149/36*z3*dR
          - 127/96*NG*dR
          - 1397/2592*NG*dR^2
          )

       + yt^2*lambda*ytau^2*h^3 * (
          - 36*dR
          )

       + yt^2*lambda*yb^2*h^3 * (
          + 111/2*dR
          - 36*dR^2
          - 72*z3*dR
          )

       + yt^2*lambda*g2^2*h^3 * (
          + 189/8*dR
          - 108*z3*dR
          )

       + yt^2*lambda*g1^2*h^3 * (
          - 73/8*dR
          - 20*z3*dR
          )

       + yt^2*lambda^2*h^3 * (
          + 99/2*dR
          )

       + yt^4*ytau^2*h^3 * (
          + 24*dR
          )

       + yt^4*yb^2*h^3 * (
          + 29/16*dR
          + 55/2*dR^2
          + 12*z3*dR
          )

       + yt^4*g2^2*h^3 * (
          - 951/16*dR
          - 9/4*dR^2
          + 54*z3*dR
          )

       + yt^4*g1^2*h^3 * (
          - 323/48*dR
          - 3/4*dR^2
          + 4*z3*dR
          )

       + yt^4*lambda*h^3 * (
          + 333/4*dR
          - 18*dR^2
          + 72*z3*dR
          )

       + yt^6*h^3 * (
          - 617/16*dR
          + 24*dR^2
          + 15*z3*dR
          )

       + gs^2*yb^4*h^3 * (
          + 447/8*cf*dR
          - 90*z3*cf*dR
          )

       + gs^2*g2^2*yb^2*h^3 * (
          - 489/16*cf*dR
          + 27*z3*cf*dR
          )

       + gs^2*g2^4*h^3 * (
          + 135/16*NG*cf*dR
          - 9*NG*z3*cf*dR
          )

       + gs^2*g1^2*yb^2*h^3 * (
          - 991/144*cf*dR
          + 5*z3*cf*dR
          )

       + gs^2*g1^4*h^3 * (
          + 55/16*NG*cf*dR
          - 11/3*NG*z3*cf*dR
          )

       + gs^2*lambda*yb^2*h^3 * (
          - 153*cf*dR
          + 144*z3*cf*dR
          )

       + gs^2*yt^2*yb^2*h^3 * (
          + 41/4*cf*dR
          - 12*z3*cf*dR
          )

       + gs^2*yt^2*g2^2*h^3 * (
          - 489/16*cf*dR
          + 27*z3*cf*dR
          )

       + gs^2*yt^2*g1^2*h^3 * (
          - 2419/144*cf*dR
          + 17*z3*cf*dR
          )

       + gs^2*yt^2*lambda*h^3 * (
          - 153*cf*dR
          + 144*z3*cf*dR
          )

       + gs^2*yt^4*h^3 * (
          + 447/8*cf*dR
          - 90*z3*cf*dR
          )

       + gs^4*yb^2*h^3 * (
          - 119/4*cf^2*dR
          + 77/2*ca*cf*dR
          + 36*z3*cf^2*dR
          - 18*z3*ca*cf*dR
          - 16*NG*tr*cf*dR
          )

       + gs^4*yt^2*h^3 * (
          - 119/4*cf^2*dR
          + 77/2*ca*cf*dR
          + 36*z3*cf^2*dR
          - 18*z3*ca*cf*dR
          - 16*NG*tr*cf*dR
			    ));
