>	restart:grtw();

 

 

 

 

 

 

	(1)

 

>	l:=1:grload(AdS,`/home/grumil/projects/BHfreeCCTMG/AdS3lccx.mpl`):

 

 

 

 

 

 

 

 

 

Calculated ds for AdS (0.000000 sec.)
	(2)

 

>	hdp:=-h-hp:

 

>	bhdp:=-bh-bhp:

 

>	grload(log1,`/home/grumil/projects/BHfreeCCTMG/Logasymptoticsimple.mpl`):

 

 

 

 

 

 

 

 

 

Calculated ds for log1 (0.000000 sec.)
	(3)

 

>	grdef(`Psi1{a b}:=g<2>{a b}`);

 

Created definition for Psi1(dn,dn)

 

>	grcalc(2=log1,Psi1(dn,dn));

 

 

Calculated Psi1(dn,dn) for log1 (0.000000 sec.)
	(4)

 

>	grmetric(AdS);

 

Default metric is now AdS.

 

>	grdef(`Psi1U{a ^b}:=Psi1{a c}*g{^c ^b}`);

 

Created definition for Psi1U(dn,up)

 

>	grdef(`Psi1UU{^a ^b}:=Psi1U{c ^b}*g{^c ^a}`);

 

Created definition for Psi1UU(up,up)

 

>	grdef(`deChr1{^a b c}:=g{^a ^d}*(Psi1{d b ;c}+Psi1{d c ;b}-Psi1{b c ;d})/2`);

 

Created a definition for Psi1(dn,dn,cdn)

 

Created definition for deChr1(up,dn,dn)

 

>	grload(log2,`/home/grumil/projects/BHfreeCCTMG/Logasymptoticprimesimple.mpl`):

 

 

 

 

 

 

 

 

 

Calculated ds for log2 (0.000000 sec.)
	(5)

 

>	grdef(`Psi2{a b}:=g<3>{a b}`);

 

Created definition for Psi2(dn,dn)

 

>	grcalc(3=log2,Psi2(dn,dn));

 

 

Calculated Psi2(dn,dn) for log2 (0.000000 sec.)
	(6)

 

>	grmetric(AdS);

 

Default metric is now AdS.

 

>	grdef(`Psi2U{a ^b}:=Psi2{a c}*g{^c ^b}`);

 

Created definition for Psi2U(dn,up)

 

>	grdef(`Psi2UU{^a ^b}:=Psi2U{c ^b}*g{^c ^a}`);

 

Created definition for Psi2UU(up,up)

 

>	grdef(`deChr2{^a b c}:=g{^a ^d}*(Psi2{d b ;c}+Psi2{d c ;b}-Psi2{b c ;d})/2`);

 

Created a definition for Psi2(dn,dn,cdn) Created definition for deChr2(up,dn,dn)

 

>	grload(right,`/home/grumil/projects/BHfreeCCTMG/Rasymptoticdoubleprimesimple.mpl`):

 

 

 

 

 

 

 

 

 

Calculated ds for right (0.000000 sec.)
	(7)

 

>	grdef(`PsiR{a b}:=g<4>{a b}`);

 

Created definition for PsiR(dn,dn)

 

>	grcalc(4=right,PsiR(dn,dn));

 

 

Calculated PsiR(dn,dn) for right (0.000000 sec.)
	(8)

 

>	grmetric(AdS);

 

Default metric is now AdS.

 

>	grdef(`PsiRU{a ^b}:=PsiR{a c}*g{^c ^b}`);

 

Created definition for PsiRU(dn,up)

 

>	grdef(`PsiRUU{^a ^b}:=PsiRU{c ^b}*g{^c ^a}`);

 

Created definition for PsiRUU(up,up)

 

>	grdef(`deChrR{^a b c}:=g{^a ^d}*(PsiR{d b ;c}+PsiR{d c ;b}-PsiR{b c ;d})/2`);

 

Created a definition for PsiR(dn,dn,cdn) Created definition for deChrR(up,dn,dn)

 

>	grload(left1,`/home/grumil/projects/BHfreeCCTMG/Lasymptoticsimple.mpl`):

 

 

 

 

 

 

 

 

 

Calculated ds for left1 (0.000000 sec.)
	(9)

 

>	grdef(`PsiL1{a b}:=g<5>{a b}`);

 

Created definition for PsiL1(dn,dn)

 

>	grcalc(5=left1,PsiL1(dn,dn));

 

 

Calculated PsiL1(dn,dn) for left1 (0.000000 sec.)
	(10)

 

>	grmetric(AdS);

 

Default metric is now AdS.

 

>	grdef(`PsiL1U{a ^b}:=PsiL1{a c}*g{^c ^b}`);

 

Created definition for PsiL1U(dn,up)

 

>	grdef(`PsiL1UU{^a ^b}:=PsiL1U{c ^b}*g{^c ^a}`);

 

Created definition for PsiL1UU(up,up)

 

>	grload(left2,`/home/grumil/projects/BHfreeCCTMG/Lasymptoticprimesimple.mpl`):

 

 

 

 

 

 

 

 

 

Calculated ds for left2 (0.000000 sec.)
	(11)

 

>	grdef(`PsiL2{a b}:=g<6>{a b}`);

 

Created definition for PsiL2(dn,dn)

 

>	grcalc(6=left2,PsiL2(dn,dn));

 

 

Calculated PsiL2(dn,dn) for left2 (0.000000 sec.)
	(12)

 

>	grmetric(AdS);

 

Default metric is now AdS.

 

>	grdef(`PsiL2U{a ^b}:=PsiL2{a c}*g{^c ^b}`);

 

Created definition for PsiL2U(dn,up)

 

>	grdef(`PsiL2UU{^a ^b}:=PsiL2U{c ^b}*g{^c ^a}`);

 

Created definition for PsiL2UU(up,up)

 

>	grdef(`deRiem1{^a b c d}:=deChr1{^a b d ;c}-deChr1{^a b c ;d}`);

 

Created a definition for deChr1(up,dn,dn,cdn) Created definition for deRiem1(up,dn,dn,dn)

 

>	grdef(`deRiem2{^a b c d}:=deChr2{^a b d ;c}-deChr2{^a b c ;d}`);

 

Created a definition for deChr2(up,dn,dn,cdn) Created definition for deRiem2(up,dn,dn,dn)

 

>	grdef(`deRiemR{^a b c d}:=deChrR{^a b d ;c}-deChrR{^a b c ;d}`);

 

Created a definition for deChrR(up,dn,dn,cdn) Created definition for deRiemR(up,dn,dn,dn)

 

>

grdef(`deRic12{(a b)}:=Psi1U{c ^d}*deRiem2{^c (a b) d}+(2*Psi1U{a ^c ;d}*Psi2{c b ;e}*g{^d ^e}-2*Psi1U{a ^c ;d}*Psi2U{b ^d ;c}+Psi1UU{^c ^d ;a}*Psi2{c d ;b})/4+Psi2U{c ^d}*deRiem1{^c (a b) d}+(2*Psi2U{a ^c ;d}*Psi1{c b ;e}*g{^d ^e}-2*Psi2U{a ^c ;d}*Psi1U{b ^d ;c}+Psi2UU{^c ^d ;a}*Psi1{c d ;b})/4`);grdef(`deRic2{(a b)}:=PsiRU{c ^d}*deRiem2{^c (a b) d}+(2*PsiRU{a ^c ;d}*Psi2{c b ;e}*g{^d ^e}-2*PsiRU{a ^c ;d}*Psi2U{b ^d ;c}+PsiRUU{^c ^d ;a}*Psi2{c d ;b})/4+Psi2U{c ^d}*deRiemR{^c (a b) d}+(2*Psi2U{a ^c ;d}*PsiR{c b ;e}*g{^d ^e}-2*Psi2U{a ^c ;d}*PsiRU{b ^d ;c}+Psi2UU{^c ^d ;a}*PsiR{c d ;b})/4`);grdef(`deRic1{(a b)}:=PsiRU{c ^d}*deRiem1{^c (a b) d}+(2*PsiRU{a ^c ;d}*Psi1{c b ;e}*g{^d ^e}-2*PsiRU{a ^c ;d}*Psi1U{b ^d ;c}+PsiRUU{^c ^d ;a}*Psi1{c d ;b})/4+Psi1U{c ^d}*deRiemR{^c (a b) d}+(2*Psi1U{a ^c ;d}*PsiR{c b ;e}*g{^d ^e}-2*Psi1U{a ^c ;d}*PsiRU{b ^d ;c}+Psi1UU{^c ^d ;a}*PsiR{c d ;b})/4`);

 

Created a definition for Psi1U(dn,up,cdn) Created a definition for Psi2U(dn,up,cdn) Created a definition for Psi1UU(up,up,cdn) Created a definition for Psi2UU(up,up,cdn)

 

Created definition for deRic12(dn,dn) Created a definition for PsiRU(dn,up,cdn) Created a definition for PsiRUU(up,up,cdn) Created definition for deRic2(dn,dn) Created definition for deRic1(dn,dn)

 

>

grdef(`De12{(a b)}:=(-(Psi2U{a ^c}*PsiL1{c b})+(PsiL1{c d}*LevC{a ^c ^e}*deChr2{^d b e})-(Psi1U{a ^c}*PsiL2{c b})+(PsiL2{c d}*LevC{a ^c ^e}*deChr1{^d b e}))`);grdef(`De2{(a b)}:=-(PsiRU{a ^c}*PsiL2{c b})+(PsiL2{c d}*LevC{a ^c ^e}*deChrR{^d b e})`);grdef(`De1{(a b)}:=-(PsiRU{a ^c}*PsiL1{c b})+(PsiL1{c d}*LevC{a ^c ^e}*deChrR{^d b e})`);

 

Created definition for LevC(dn,up,up) Created definition for De12(dn,dn) Created definition for De2(dn,dn) Created definition for De1(dn,dn)

 

>

grcalcalter(detg,simplify);grcalcalter(g(up,up),simplify);grcalcalter(Psi1(dn,dn),simplify);grcalcalter(Psi1U(dn,up),simplify);grcalcalter(Psi1UU(up,up),simplify);grcalcalter(Psi2(dn,dn),simplify);grcalcalter(Psi2U(dn,up),simplify);grcalcalter(Psi2UU(up,up),simplify);grcalcalter(PsiL1(dn,dn),simplify);grcalcalter(PsiL1U(dn,up),simplify);grcalcalter(PsiL1UU(up,up),simplify);grcalcalter(PsiL2(dn,dn),simplify);grcalcalter(PsiL2U(dn,up),simplify);grcalcalter(PsiL2UU(up,up),simplify);grcalcalter(Chr(dn,dn,dn),simplify);grcalcalter(Chr(dn,dn,up),simplify);grcalcalter(deChr1(up,dn,dn),simplify);grcalcalter(deChr2(up,dn,dn),simplify);grcalcalter(deChrR(up,dn,dn),simplify);grcalcalter(De2(dn,dn),simplify);grcalcalter(De1(dn,dn),simplify);grcalcalter(De12(dn,dn),simplify);

 

 

Simplification will be applied during calculation.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Applying routine simplify to object detg Calculated detg for AdS (0.000000 sec.)
Simplification will be applied during calculation. Applying routine simplify to object g(up,up) Calculated g(up,up) for AdS (0.004000 sec.)
Simplification will be applied during calculation. Applying routine simplify to object Psi1(dn,dn) Calculated Psi1(dn,dn) for AdS (0.016000 sec.)
Simplification will be applied during calculation.
Applying routine simplify to object Psi1U(dn,up) Calculated Psi1U(dn,up) for AdS (0.028000 sec.)
Simplification will be applied during calculation. Applying routine simplify to object Psi1UU(up,up) Calculated Psi1UU(up,up) for AdS (0.040000 sec.)
Simplification will be applied during calculation. Applying routine simplify to object Psi2(dn,dn) Calculated Psi2(dn,dn) for AdS (0.012000 sec.)
Simplification will be applied during calculation. Applying routine simplify to object Psi2U(dn,up) Calculated Psi2U(dn,up) for AdS (0.024000 sec.)
Simplification will be applied during calculation. Applying routine simplify to object Psi2UU(up,up) Calculated Psi2UU(up,up) for AdS (0.036000 sec.)
Simplification will be applied during calculation. Applying routine simplify to object PsiL1(dn,dn) Calculated PsiL1(dn,dn) for AdS (0.008000 sec.)
Simplification will be applied during calculation. Applying routine simplify to object PsiL1U(dn,up) Calculated PsiL1U(dn,up) for AdS (0.008000 sec.)
Simplification will be applied during calculation. Applying routine simplify to object PsiL1UU(up,up) Calculated PsiL1UU(up,up) for AdS (0.012000 sec.)
Simplification will be applied during calculation. Applying routine simplify to object PsiL2(dn,dn) Calculated PsiL2(dn,dn) for AdS (0.008000 sec.)
Simplification will be applied during calculation. Applying routine simplify to object PsiL2U(dn,up) Calculated PsiL2U(dn,up) for AdS (0.012000 sec.)
Simplification will be applied during calculation. Applying routine simplify to object PsiL2UU(up,up) Calculated PsiL2UU(up,up) for AdS (0.016000 sec.)
Simplification will be applied during calculation. Applying routine simplify to object g(dn,dn,pdn) Calculated g(dn,dn,pdn) for AdS (0.000000 sec.) Applying routine simplify to object Chr(dn,dn,dn) Calculated Chr(dn,dn,dn) for AdS (0.000000 sec.)
Simplification will be applied during calculation. Applying routine simplify to object Chr(dn,dn,up) Calculated Chr(dn,dn,up) for AdS (0.000000 sec.)
Simplification will be applied during calculation. Applying routine simplify to object Psi1(dn,dn,cdn) Calculated Psi1(dn,dn,cdn) for AdS (0.136000 sec.)
Applying routine simplify to object deChr1(up,dn,dn) Calculated deChr1(up,dn,dn) for AdS (0.268000 sec.)
Simplification will be applied during calculation. Applying routine simplify to object Psi2(dn,dn,cdn) Calculated Psi2(dn,dn,cdn) for AdS (0.128000 sec.) Applying routine simplify to object deChr2(up,dn,dn) Calculated deChr2(up,dn,dn) for AdS (0.256000 sec.)
Simplification will be applied during calculation. Applying routine simplify to object PsiR(dn,dn) Calculated PsiR(dn,dn) for AdS (0.008000 sec.) Applying routine simplify to object PsiR(dn,dn,cdn) Calculated PsiR(dn,dn,cdn) for AdS (0.076000 sec.) Applying routine simplify to object deChrR(up,dn,dn) Calculated deChrR(up,dn,dn) for AdS (0.172000 sec.)
Simplification will be applied during calculation. Applying routine simplify to object LevC(up,up,up) Calculated LevC(up,up,up) for AdS (0.000000 sec.) Applying routine simplify to object LevC(dn,up,up) Calculated LevC(dn,up,up) for AdS (0.000000 sec.) Applying routine simplify to object PsiRU(dn,up) Calculated PsiRU(dn,up) for AdS (0.016000 sec.) Applying routine simplify to object De2(dn,dn)
Calculated De2(dn,dn) for AdS (0.796000 sec.)
Simplification will be applied during calculation. Applying routine simplify to object De1(dn,dn)
Calculated De1(dn,dn) for AdS (0.764000 sec.)
Simplification will be applied during calculation. Applying routine simplify to object De12(dn,dn)
Calculated De12(dn,dn) for AdS (1.556000 sec.)
	(13)

 

>	grcalcalter(deRiem1(up,dn,dn,dn),simplify);grcalcalter(deRiem2(up,dn,dn,dn),simplify);grcalcalter(deRiemR(up,dn,dn,dn),simplify);grcalcalter(deRic2(dn,dn),simplify);grcalcalter(deRic1(dn,dn),simplify);grcalcalter(deRic12(dn,dn),simplify);

 

 

Simplification will be applied during calculation.

 

Applying routine simplify to object deChr1(up,dn,dn,cdn)

 

Calculated deChr1(up,dn,dn,cdn) for AdS (1.968000 sec.) Applying routine simplify to object deRiem1(up,dn,dn,dn)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Calculated deRiem1(up,dn,dn,dn) for AdS (1.336000 sec.)
Simplification will be applied during calculation. Applying routine simplify to object deChr2(up,dn,dn,cdn)
Calculated deChr2(up,dn,dn,cdn) for AdS (1.984000 sec.) Applying routine simplify to object deRiem2(up,dn,dn,dn)
Calculated deRiem2(up,dn,dn,dn) for AdS (1.336000 sec.)
Simplification will be applied during calculation. Applying routine simplify to object deChrR(up,dn,dn,cdn)
Calculated deChrR(up,dn,dn,cdn) for AdS (1.680000 sec.) Applying routine simplify to object deRiemR(up,dn,dn,dn)
Calculated deRiemR(up,dn,dn,dn) for AdS (1.389000 sec.)
Simplification will be applied during calculation. Applying routine simplify to object PsiRU(dn,up,cdn)
Calculated PsiRU(dn,up,cdn) for AdS (0.212000 sec.) Applying routine simplify to object Psi2U(dn,up,cdn)
Calculated Psi2U(dn,up,cdn) for AdS (0.240000 sec.) Applying routine simplify to object Psi2UU(up,up,cdn)
Calculated Psi2UU(up,up,cdn) for AdS (0.216000 sec.) Applying routine simplify to object PsiRUU(up,up) Calculated PsiRUU(up,up) for AdS (0.020000 sec.) Applying routine simplify to object PsiRUU(up,up,cdn)
Calculated PsiRUU(up,up,cdn) for AdS (0.168000 sec.) Applying routine simplify to object deRic2(dn,dn)
Calculated deRic2(dn,dn) for AdS (40.406000 sec.)
Simplification will be applied during calculation. Applying routine simplify to object Psi1U(dn,up,cdn)
Calculated Psi1U(dn,up,cdn) for AdS (0.276000 sec.) Applying routine simplify to object Psi1UU(up,up,cdn)
Calculated Psi1UU(up,up,cdn) for AdS (0.244000 sec.) Applying routine simplify to object deRic1(dn,dn)
Calculated deRic1(dn,dn) for AdS (41.507000 sec.)
Simplification will be applied during calculation. Applying routine simplify to object deRic12(dn,dn)
Calculated deRic12(dn,dn) for AdS (41.630000 sec.)
	(14)

 

>

 

>	grdef(`Corpart1:=2*c*PsiL1UU{^a ^b}*deRic2{(a b)}`);

 

Created definition for Corpart1

 

>	grdef(`Corpart2:=2*c*Psi1UU{^a ^b}*De2{(a b)}`);

 

Created definition for Corpart2

 

>	grdef(`Corpart3:=2*c*PsiL2UU{^a ^b}*deRic1{(a b)}`);

 

Created definition for Corpart3

 

>	grdef(`Corpart4:=2*c*Psi2UU{^a ^b}*De1{(a b)}`);

 

Created definition for Corpart4

 

>	grdef(`Corpart5:=2*PsiRUU{^a ^b}*deRic12{(a b)}`);

 

Created definition for Corpart5

 

>	grdef(`Corpart6:=2*c*PsiRUU{^a ^b}*De12{(a b)}`);

 

Created definition for Corpart6

 

>	grcalcalter(Corpart1,simplify);

 

 

Simplification will be applied during calculation. Applying routine simplify to object Corpart1

 

Calculated Corpart1 for AdS (61.728000 sec.)
	(15)

 

>	grcalcalter(Corpart2,simplify);

 

 

Simplification will be applied during calculation. Applying routine simplify to object Corpart2

 

Calculated Corpart2 for AdS (8.329000 sec.)
	(16)

 

>	grcalcalter(Corpart3,simplify);

 

 

Simplification will be applied during calculation. Applying routine simplify to object Corpart3

 

Calculated Corpart3 for AdS (61.724000 sec.)
	(17)

 

>	grcalcalter(Corpart4,simplify);

 

 

Simplification will be applied during calculation. Applying routine simplify to object Corpart4

 

Calculated Corpart4 for AdS (8.464000 sec.)
	(18)

 

>	grcalcalter(Corpart5,simplify);

 

 

Simplification will be applied during calculation. Applying routine simplify to object Corpart5

 

Calculated Corpart5 for AdS (521.497000 sec.)
	(19)

 

>	grcalcalter(Corpart6,simplify);

 

 

Simplification will be applied during calculation. Applying routine simplify to object Corpart6

 

Calculated Corpart6 for AdS (10.036000 sec.)
	(20)

 

>	integrandCorpart1:=series(simplify(subs(u=t+phi,(subs(v=t-phi,(grcomponent(Corpart1))*2*Pi/8)))),xx=infinity,11);

 

(21)

 

>	integrandCorpart2:=series(simplify(subs(u=t+phi,(subs(v=t-phi,(grcomponent(Corpart2))*2*Pi/8)))),xx=infinity,11);

 

(22)

 

>	integrandCorpart3:=series(simplify(subs(u=t+phi,(subs(v=t-phi,(grcomponent(Corpart3))*2*Pi/8)))),xx=infinity,11);

 

(23)

 

>	integrandCorpart4:=series(simplify(subs(u=t+phi,(subs(v=t-phi,(grcomponent(Corpart4))*2*Pi/8)))),xx=infinity,11);

 

(24)

 

>	integrandCorpart5:=series(simplify(subs(u=t+phi,(subs(v=t-phi,(grcomponent(Corpart5))*2*Pi/8)))),xx=infinity,11);

 

(25)

 

>	integrandCorpart6:=series(simplify(subs(u=t+phi,(subs(v=t-phi,(grcomponent(Corpart6))*2*Pi/8)))),xx=infinity,11);

 

(26)

 

>	fullintegrand:=collect(collect(simplify(integrandCorpart1+integrandCorpart2+integrandCorpart3+integrandCorpart4+integrandCorpart5+integrandCorpart6),log(xx)),xx);

 

(27)

 

>

 

>

LO:=simplify((128*c^2*Pi*ln(xx)^2+(2*(-64*c^2*Pi*bh*bhp*hp^2+(256*I)*c^2*Pi*bhp*bh*t+64*c^2*Pi*bh*bhp*hp*h+96*c^2*Pi*bh*bhp+32*h*c^2*Pi*bh*bhp+32*c^2*Pi*bh*bhp*hp-64*c^2*Pi*bh*bhp*h^2))*ln(xx)/(bh*bhp)+(2*(32*h*c^2*Pi*bh*bhp-256*c^2*bh*bhp*Pi*t^2+(64*I)*c^2*bh*bhp*Pi*h*t-128*c^2*Pi*bh*bhp*hp*h+144*c^2*Pi*bh*bhp*h^2+144*c^2*Pi*bh*bhp*hp^2+(192*I)*c^2*Pi*bhp*bh*t+(64*I)*hp*c^2*Pi*bhp*bh*t+32*c^2*Pi*bh*bhp+(128*I)*c^2*bh*bhp*Pi*h*hp*t-(128*I)*c^2*Pi*bh*hp^2*bhp*t-(128*I)*c^2*bh*bhp*Pi*h^2*t+32*c^2*Pi*bh*bhp*hp))/(bh*bhp))/xx);

 

(28)

 

>

NLO:=collect(collect(simplify((((2*(128*c^2*Pi*hp^3*bh^2*bhp+96*c^2*Pi*bh*bhp^2*hp+96*c^2*Pi*h^2*bh*bhp^2*hp+96*c^2*Pi*bh^2*bhp*h+128*c^2*Pi*bh*h^3*bhp^2+32*c^2*Pi*hp^3*bh*bhp^2-32*c^2*Pi*bh^2*bhp*hp*h^2-64*c^2*Pi*bh^2*bhp*hp-64*c^2*Pi*bh*bhp^2*h-32*c^2*Pi*bh*bhp^2*h*hp^2+32*c^2*Pi*h^3*bh^2*bhp+96*c^2*Pi*hp^2*h*bh^2*bhp))*ln(xx)^2/(bh*bhp)+(2*((384*I)*c^2*Pi*bh*bhp^2*hp*t+(512*I)*c^2*Pi*h^3*bhp^2*bh*t+(512*I)*c^2*Pi*hp^3*bh^2*bhp*t-16*c^2*Pi*bh^2*bhp*h+80*c^2*Pi*bh^2*bhp*hp+96*c^2*Pi*hp*bhp*h*bh^2+(128*I)*c^2*Pi*bh*hp^3*bhp^2*t-(256*I)*c^2*Pi*bh^2*bhp*hp*t+(128*I)*c^2*Pi*bhp*h^3*bh^2*t+272*c^2*Pi*bh^2*bhp*hp*h^2-544*c^2*Pi*h^2*bh*bhp^2*hp-544*c^2*Pi*hp^2*h*bh^2*bhp+48*c^2*Pi*hp^3*bh+96*c^2*Pi*bh^3*hp^3+96*h^3*c^2*Pi*bhp^3+64*c^2*Pi*bh*h^3*bhp-208*c^2*Pi*h^3*bh^2*bhp-208*c^2*Pi*hp^3*bh*bhp^2-(128*I)*h^2*c^2*Pi*bh^2*bhp*hp*t+(384*I)*h*c^2*Pi*bh^2*bhp*hp^2*t+(384*I)*h^2*c^2*Pi*bh*bhp^2*hp*t-(256*I)*h*c^2*Pi*bh*bhp^2*t+48*c^2*Pi*h^3*bhp+64*c^2*Pi*bhp*hp^3*bh-64*c^2*Pi*bh*h^2*hp*bhp-64*c^2*Pi*bhp*hp^2*h*bh+272*c^2*Pi*bh*bhp^2*h*hp^2-96*c^2*Pi*h^2*bh*bhp^2+96*c^2*Pi*hp*bhp^2*bh*h-96*c^2*Pi*hp^2*bhp^2*bh-96*c^2*Pi*hp^2*bhp*bh^2-400*c^2*Pi*hp^3*bh^2*bhp-400*c^2*Pi*bh*h^3*bhp^2-16*c^2*Pi*bh*bhp^2*hp+80*c^2*Pi*bh*bhp^2*h-96*c^2*Pi*h^2*bh^2*bhp-(128*I)*h*c^2*Pi*bh*bhp^2*hp^2*t+(384*I)*h*c^2*Pi*bh^2*bhp*t))*ln(xx)/(bh*bhp)+(2*(-(192*I)*c^2*Pi*bh^2*bhp*h^2*t-(128*I)*h*c^2*Pi*bh*hp^2*bhp*t+(128*I)*h^3*c^2*Pi*bhp*bh*t+(160*I)*c^2*Pi*bh^2*bhp*hp*t+(128*I)*c^2*Pi*bh^3*hp^3*t-32*c^2*Pi*h^3*bhp^2-32*c^2*Pi*hp^3*bh^2+32*c^2*Pi*bh^2*bhp*h+32*c^2*Pi*bh^2*bhp*hp-96*c^2*Pi*hp*bhp*h*bh^2-(1088*I)*h*c^2*Pi*bh^2*bhp*hp^2*t-432*c^2*Pi*bh^2*bhp*hp*h^2+672*c^2*Pi*h^2*bh*bhp^2*hp+672*c^2*Pi*hp^2*h*bh^2*bhp-(672*I)*c^2*Pi*h^3*bhp^2*bh*t+48*c^2*Pi*hp^3*bh-64*c^2*Pi*bh^2*bhp-160*c^2*Pi*bh^3*hp^3-160*h^3*c^2*Pi*bhp^3+(128*I)*c^2*Pi*bh*hp^3*bhp*t+32*c^2*Pi*hp^3*bh*ln(-hp*bhp)*bhp^2-64*c^2*Pi*bh*h^3*bhp-(672*I)*c^2*Pi*hp^3*bh^2*bhp*t-32*c^2*Pi*bh^3*hp^3*ln(-hp*bhp)+48*c^2*Pi*h^3*bh^2*bhp+48*c^2*Pi*hp^3*bh*bhp^2-512*c^2*Pi*hp^3*bh^2*bhp*t^2+256*c^2*Pi*hp*bhp*bh^2*t^2-384*c^2*Pi*hp*bhp^2*bh*t^2+64*c^2*Pi*hp^3*ln(-hp*bhp)*bh^2+128*h*c^2*Pi*hp^2*bhp^2*bh*t^2-(32*I)*h*c^2*Pi*bh^2*bhp*t+128*h^2*c^2*Pi*hp*bhp*bh^2*t^2-384*h^2*c^2*Pi*hp*bhp^2*bh*t^2-384*h*c^2*Pi*bh^2*bhp*t^2-384*h*c^2*Pi*hp^2*bhp*bh^2*t^2+64*h^3*c^2*Pi*ln(-h*bh)*bhp^2*bh-128*c^2*Pi*hp^3*bhp^2*bh*t^2+(192*I)*c^2*Pi*bh^2*bhp*hp*h*t+(160*I)*h*c^2*Pi*bh*bhp^2*t-(128*I)*h^2*hp*c^2*Pi*bhp*bh*t+(128*I)*c^2*Pi*hp^3*bh^2*t-(352*I)*c^2*Pi*bhp*h^3*bh^2*t-(352*I)*c^2*Pi*bh*hp^3*bhp^2*t+(96*I)*c^2*Pi*bh*hp^3*t-(192*I)*c^2*Pi*bh^2*bhp*hp^2*t-(1088*I)*h^2*c^2*Pi*bh*bhp^2*hp*t-(32*I)*c^2*Pi*bh*bhp^2*hp*t+(544*I)*h^2*c^2*Pi*bh^2*bhp*hp*t+256*h*c^2*Pi*bh*bhp^2*t^2-144*h*c^2*Pi*bh*bhp+48*c^2*Pi*h^3*bhp-512*h^3*c^2*Pi*bh*bhp^2*t^2-128*h^3*c^2*Pi*bh^2*bhp*t^2+64*h^3*c^2*Pi*ln(-h*bh)*bhp^2+32*h^3*c^2*Pi*ln(-h*bh)*bhp*bh^2-32*h^3*c^2*Pi*ln(-h*bh)*bhp^3-64*c^2*Pi*bhp*hp^3*bh-144*c^2*Pi*bh*bhp*hp+64*c^2*Pi*hp^3*ln(-hp*bhp)*bh^2*bhp+224*c^2*Pi*bh*h^2*hp*bhp+224*c^2*Pi*bhp*hp^2*h*bh+(192*I)*c^2*Pi*bh*bhp^2*h*hp*t-(192*I)*c^2*Pi*bh*bhp^2*h^2*t-64*c^2*Pi*bh*bhp^2-432*c^2*Pi*bh*bhp^2*h*hp^2+(544*I)*h*c^2*Pi*bh*bhp^2*hp^2*t+(128*I)*h^3*c^2*Pi*bhp^3*t-32*c^2*Pi*h^2*bh*bhp^2-96*c^2*Pi*hp*bhp^2*bh*h+160*c^2*Pi*hp^2*bhp^2*bh-128*c^2*Pi*bh*bhp*hp*h-32*c^2*Pi*hp^2*bhp*bh^2+80*c^2*Pi*hp^3*bh^2*bhp+64*c^2*Pi*bh*bhp*hp^2+80*c^2*Pi*bh*h^3*bhp^2+32*c^2*Pi*bh*bhp^2*hp+32*c^2*Pi*bh*bhp^2*h+160*c^2*Pi*h^2*bh^2*bhp+64*c^2*Pi*bh*bhp*h^2+(96*I)*c^2*Pi*bhp*h^3*t+(128*I)*c^2*Pi*h^3*bhp^2*t-(192*I)*c^2*Pi*bh*bhp^2*hp^2*t))/(bh*bhp))/xx^2)),t),log(xx));

 

(29)

 

>	series(subs(h=a*h,subs(bh=a*bh,subs(hp=a*hp,subs(bhp=a*bhp,NLO)))),a=infinity,3);

 

(30)

 

>	NLOlogxsquared:=simplify(32*c^2*Pi*(-2*bh*bhp^2*h*hp^2+6*hp^2*h*bh^2*bhp-2*bh^2*bhp*hp*h^2+8*bh*h^3*bhp^2+2*h^3*bh^2*bhp+6*h^2*bh*bhp^2*hp+8*hp^3*bh^2*bhp+2*hp^3*bh*bhp^2)/bh/bhp);

 

(31)

 

>	NLOtsquared:=simplify(32*c^2*Pi*(-24*h^2*bh*bhp^2*hp-24*hp^2*h*bh^2*bhp-8*hp^3*bh*bhp^2-32*hp^3*bh^2*bhp+8*bh*bhp^2*h*hp^2+8*bh^2*bhp*hp*h^2-32*bh*h^3*bhp^2-8*h^3*bh^2*bhp)/bh/bhp);

 

(32)

 

>	NLOlogxt:=simplify(32*c^2*Pi*(-(8*I)*h*bh*bhp^2*hp^2+(8*I)*bhp*h^3*bh^2-(8*I)*h^2*bh^2*bhp*hp+(8*I)*bh*hp^3*bhp^2+(32*I)*h^3*bhp^2*bh+(24*I)*h*bh^2*bhp*hp^2+(24*I)*h^2*bh*bhp^2*hp+(32*I)*hp^3*bh^2*bhp)/bh/bhp);

 

(33)

 

>	NLOlogx:=simplify(32*c^2*Pi*(6*bh^3*hp^3-25*hp^3*bh^2*bhp+17*bh*bhp^2*h*hp^2+6*h^3*bhp^3-13*hp^3*bh*bhp^2-13*h^3*bh^2*bhp-25*bh*h^3*bhp^2+17*bh^2*bhp*hp*h^2-34*h^2*bh*bhp^2*hp-34*hp^2*h*bh^2*bhp)/bh/bhp);

 

(34)

 

>	noncontactNLOlogx:=simplify(32*c^2*Pi*(6*h^3*bhp^3+6*bh^3*hp^3)/bh/bhp);

 

(35)

 

>	NLOt:=simplify(32*c^2*Pi*((8*I)*h^3*bhp^3+(34*I)*h^2*bh^2*bhp*hp-(22*I)*bhp*h^3*bh^2+(8*I)*bh^3*hp^3-(68*I)*h^2*bh*bhp^2*hp+(34*I)*h*bh*bhp^2*hp^2-(68*I)*h*bh^2*bhp*hp^2-(42*I)*hp^3*bh^2*bhp-(42*I)*h^3*bhp^2*bh-(22*I)*bh*hp^3*bhp^2)/bh/bhp);

 

(36)

 

>	noncontactNLOt:=simplify(64*I*c^2*Pi*(4*h^3*bhp^3+4*bh^3*hp^3)/bh/bhp);

 

(37)

 

>

NLOrest:=collect(collect(simplify(32*c^2*Pi*(-10*bh^3*hp^3+4*hp^3*ln(-a^2*hp*bhp)*bh^2*bhp+2*hp^3*bh*ln(-a^2*hp*bhp)*bhp^2-10*h^3*bhp^3+3*hp^3*bh*bhp^2+3*h^3*bh^2*bhp-27*bh*bhp^2*h*hp^2+5*bh*h^3*bhp^2-27*bh^2*bhp*hp*h^2+42*h^2*bh*bhp^2*hp+5*hp^3*bh^2*bhp+42*hp^2*h*bh^2*bhp-2*bh^3*hp^3*ln(-a^2*hp*bhp)-2*h^3*ln(-a^2*h*bh)*bhp^3+2*h^3*ln(-a^2*h*bh)*bhp*bh^2+4*h^3*ln(-a^2*h*bh)*bhp^2*bh)/bh/bhp),log(-a^2*hp*bhp)),log(-a^2*h*bh));

 

(38)

 

>	noncontactNLOrest:=simplify(32*Pi*c^2*((-2*h^3*bhp^3)*log(-a^2*h*bh)+(-2*bh^3*hp^3)*log(-a^2*hp*bhp)+(-10*bh^3*hp^3-10*h^3*bhp^3))/bh/bhp);

 

(39)

 

>

 

(40)

 

>

 

>

 

>