coeffrr = 0 ; 
coeffaa = 0 ; 
coeffra = 0 ; 

dEqXr[0][0] = 0 ;
dEqXr[0][1] = 0 ;
dEqXr[0][2] = 0 ;
dEqXr[0][3] = 0 ;
dEqXr[0][4] = 0 ;
dEqXr[1][0] = 0 ;
dEqXr[1][1] = 0 ;
dEqXr[1][2] = 0 ;
dEqXr[1][3] = 0 ;
dEqXr[1][4] = 0 ;
dEqXr[2][0] = 0 ;
dEqXr[2][1] = 0 ;
dEqXr[2][2] = 0 ;
dEqXr[2][3] = 0 ;
dEqXr[2][4] = 0 ;
dEqXr[3][0] = 0 ;
dEqXr[3][1] = 0 ;
dEqXr[3][2] = 0 ;
dEqXr[3][3] = 0 ;
dEqXr[3][4] = 0 ;
dEqXr[4][0] = 0 ;
dEqXr[4][1] = (Ar*Power(r,4)*Power(-1 + r2,4)*Power(rh,2)*expl(-(A*r2) - T))/8. - (Br*Power(-1 + r2,4)*r2*Power(rh,2)*expl(-(A*r2) - T))/8. - (3*Power(-1 + r2,4)*r2*Power(rh,2)*Sr*expl(-(A*r2) - T))/8. - (Fa*Power(r,3)*Power(-1 + r2,3)*Power(rh,2)*expl(-B - A*r2 - S - T))/16. + (r*Power(-1 + r2,3)*Power(rh,2)*expl(-B - A*r2 - S - T)*((1 + A*Power(r,4) - (-5 + A)*r2)*expl(B + S) + (-1 + r2)*expl(B + A*r2 + S) - 2*r2*expl(A*r2 + T) - 4*r2*expl(B + A*r2 + T)))/4. ;
dEqXr[4][2] = (Br*Power(-1 + r2,4)*Power(rh,2)*expl(-(A*r2) - T))/8. - (Ar*Power(-1 + r2,4)*r2*Power(rh,2)*expl(-(A*r2) - T))/8. + (3*Power(-1 + r2,4)*Power(rh,2)*Sr*expl(-(A*r2) - T))/8. + (Fa*r*Power(-1 + r2,3)*Power(rh,2)*expl(-B - A*r2 - S - T))/16. - (Power(-1 + r2,3)*Power(rh,2)*expl(-B - A*r2 - S - T)*((1 + A*Power(r,4) - (-5 + A)*r2)*expl(B + S) + (-1 + r2)*expl(B + A*r2 + S) - 2*r2*expl(A*r2 + T) - 4*r2*expl(B + A*r2 + T)))/(4.*r) ;
dEqXr[4][3] = 0 ;
dEqXr[4][4] = (3*Br*Power(-1 + r2,4)*Power(rh,2)*expl(-(A*r2) - T))/8. - (3*Ar*Power(-1 + r2,4)*r2*Power(rh,2)*expl(-(A*r2) - T))/8. + (9*Power(-1 + r2,4)*Power(rh,2)*Sr*expl(-(A*r2) - T))/8. + (3*Fa*r*Power(-1 + r2,3)*Power(rh,2)*expl(-B - A*r2 - S - T))/16. - (3*Power(-1 + r2,3)*Power(rh,2)*expl(-B - A*r2 - S - T)*((1 + A*Power(r,4) - (-5 + A)*r2)*expl(B + S) + (-1 + r2)*expl(B + A*r2 + S) - 2*r2*expl(A*r2 + T) - 4*r2*expl(B + A*r2 + T)))/(4.*r) ;

dEqXa[0][0] = 1 ;
dEqXa[0][1] = 0 ;
dEqXa[0][2] = 0 ;
dEqXa[0][3] = 0 ;
dEqXa[0][4] = 0 ;
dEqXa[1][0] = 0 ;
dEqXa[1][1] = 0 ;
dEqXa[1][2] = 0 ;
dEqXa[1][3] = 0 ;
dEqXa[1][4] = 1 ;
dEqXa[2][0] = 0 ;
dEqXa[2][1] = 1 ;
dEqXa[2][2] = 0 ;
dEqXa[2][3] = 0 ;
dEqXa[2][4] = 0 ;
dEqXa[3][0] = 0 ;
dEqXa[3][1] = 0 ;
dEqXa[3][2] = 0 ;
dEqXa[3][3] = 0 ;
dEqXa[3][4] = 0 ;
dEqXa[4][0] = -(Ba*Power(-1 + r2,2)*Power(rh,2)*expl(-B - S))/4. + (Aa*Power(-1 + r2,2)*r2*Power(rh,2)*expl(-B - S))/4. + (Power(-1 + r2,2)*Power(rh,2)*Sa*expl(-B - S))/4. + (Power(-1 + r2,2)*Power(rh,2)*Ta*expl(-B - S))/2. + ((-1 + r2)*rh*expl(-B - A*r2 - S - T)*(2*expl(B + S) + 2*expl(B + A*r2 + S) + (-8 + B*(rh - r2*rh))*expl(A*r2 + T) + 4*expl(B + A*r2 + T)))/2. ;
dEqXa[4][1] = (Aa*Power(r,4)*Power(-1 + r2,2)*Power(rh,2)*expl(-B - S))/8. - (Ba*Power(-1 + r2,2)*r2*Power(rh,2)*expl(-B - S))/8. + (Power(-1 + r2,2)*r2*Power(rh,2)*Sa*expl(-B - S))/8. + (Power(-1 + r2,2)*r2*Power(rh,2)*Ta*expl(-B - S))/4. + ((-1 + r2)*r2*rh*expl(-B - A*r2 - S - T)*(2*expl(B + S) + 2*expl(B + A*r2 + S) + (-8 + B*(rh - r2*rh))*expl(A*r2 + T) + 4*expl(B + A*r2 + T)))/4. ;
dEqXa[4][2] = (Ba*Power(-1 + r2,2)*Power(rh,2)*expl(-B - S))/8. - (Aa*Power(-1 + r2,2)*r2*Power(rh,2)*expl(-B - S))/8. - (Power(-1 + r2,2)*Power(rh,2)*Sa*expl(-B - S))/8. - (Power(-1 + r2,2)*Power(rh,2)*Ta*expl(-B - S))/4. + ((-1 + r2)*rh*expl(-B - A*r2 - S - T)*(-2*expl(B + S) - 2*expl(B + A*r2 + S) + (8 + B*(-1 + r2)*rh)*expl(A*r2 + T) - 4*expl(B + A*r2 + T)))/4. ;
dEqXa[4][3] = (Fa*Power(-1 + r2,2)*r2*Power(rh,2)*expl(-2*B - A*r2 - 2*S - T))/32. + (Br*r*Power(-1 + r2,3)*Power(rh,2)*expl(-B - A*r2 - S - T))/16. - (Ar*Power(r,3)*Power(-1 + r2,3)*Power(rh,2)*expl(-B - A*r2 - S - T))/16. + (3*r*Power(-1 + r2,3)*Power(rh,2)*Sr*expl(-B - A*r2 - S - T))/16. - (Power(-1 + r2,2)*Power(rh,2)*expl(-2*B - A*r2 - 2*S - T)*((1 + A*Power(r,4) - (-5 + A)*r2)*expl(B + S) + (-1 + r2)*expl(B + A*r2 + S) - 2*r2*expl(A*r2 + T) - 4*r2*expl(B + A*r2 + T)))/8. ;
dEqXa[4][4] = -(Ba*Power(-1 + r2,2)*Power(rh,2)*expl(-B - S))/8. + (Aa*Power(-1 + r2,2)*r2*Power(rh,2)*expl(-B - S))/8. + (Power(-1 + r2,2)*Power(rh,2)*Sa*expl(-B - S))/8. + (Power(-1 + r2,2)*Power(rh,2)*Ta*expl(-B - S))/4. + ((-1 + r2)*rh*expl(-B - A*r2 - S - T)*(2*expl(B + S) + 2*expl(B + A*r2 + S) + (-8 + B*(rh - r2*rh))*expl(A*r2 + T) + 4*expl(B + A*r2 + T)))/4. ;

dEqT[0] = 0 ;
dEqT[1] = 0 ;
dEqT[2] = 0 ;
dEqT[3] = 0 ;
dEqT[4] = -(Power(Br,2)*Power(-1 + r2,4)*Power(rh,2)*expl(-(A*r2) - T))/16. - (Power(Ar,2)*Power(r,4)*Power(-1 + r2,4)*Power(rh,2)*expl(-(A*r2) - T))/16. + (Ar*Br*Power(-1 + r2,4)*r2*Power(rh,2)*expl(-(A*r2) - T))/8. - (3*Br*Power(-1 + r2,4)*Power(rh,2)*Sr*expl(-(A*r2) - T))/8. + (3*Ar*Power(-1 + r2,4)*r2*Power(rh,2)*Sr*expl(-(A*r2) - T))/8. - (9*Power(-1 + r2,4)*Power(rh,2)*Power(Sr,2)*expl(-(A*r2) - T))/16. + (Ba*(-1 + r2)*rh*(1 + expl(A*r2))*expl(-(A*r2) - T))/2. - (Aa*(-1 + r2)*r2*rh*(1 + expl(A*r2))*expl(-(A*r2) - T))/2. - ((-1 + r2)*rh*Sa*(1 + expl(A*r2))*expl(-(A*r2) - T))/2. - (-1 + r2)*rh*Ta*(1 + expl(A*r2))*expl(-(A*r2) - T) + (Br*Power(-1 + r2,3)*Power(rh,2)*(1 + A*Power(r,4) - (-5 + A)*r2 + (-1 + r2)*expl(A*r2))*expl(-(A*r2) - T))/(4.*r) - (Ar*r*Power(-1 + r2,3)*Power(rh,2)*(1 + A*Power(r,4) - (-5 + A)*r2 + (-1 + r2)*expl(A*r2))*expl(-(A*r2) - T))/4. + (3*Power(-1 + r2,3)*Power(rh,2)*Sr*(1 + A*Power(r,4) - (-5 + A)*r2 + (-1 + r2)*expl(A*r2))*expl(-(A*r2) - T))/(4.*r) - (Power(Fa,2)*Power(-1 + r2,2)*r2*Power(rh,2)*expl(-2*B - A*r2 - 2*S - T))/64. - (Br*Fa*r*Power(-1 + r2,3)*Power(rh,2)*expl(-B - A*r2 - S - T))/16. + (Ar*Fa*Power(r,3)*Power(-1 + r2,3)*Power(rh,2)*expl(-B - A*r2 - S - T))/16. - (3*Fa*r*Power(-1 + r2,3)*Power(rh,2)*Sr*expl(-B - A*r2 - S - T))/16. + (Fa*Power(-1 + r2,2)*Power(rh,2)*(1 + A*Power(r,4) - (-5 + A)*r2 + (-1 + r2)*expl(A*r2))*expl(-B - A*r2 - S - T))/8. + (Power(rh - r2*rh,6)*expl(-2*(B + A*r2 + S + T))*(-4*expl(B + S)*Power(4*expl(B + S) + 4*expl(B + A*r2 + S) - 2*(8 + B*(-1 + r2)*rh)*expl(A*r2 + T) + 8*expl(B + A*r2 + T),2) + 32*Power(-1 + r2,2)*Power(rh,2)*(1 + 2*expl(B))*expl(2*(A*r2 + T))*((-1 - A*Power(r,4) + (-5 + A)*r2)*expl(B + S) - (-1 + r2)*expl(B + A*r2 + S) + 2*r2*expl(A*r2 + T) + 4*r2*expl(B + A*r2 + T)) - (8*Power(-1 + r2,2)*Power(rh,2)*expl(A*r2 + T)*Power((-1 - A*Power(r,4) + (-5 + A)*r2)*expl(B + S) - (-1 + r2)*expl(B + A*r2 + S) + 2*r2*expl(A*r2 + T) + 4*r2*expl(B + A*r2 + T),2))/Power(r,2) + 8*(-8 + B*(rh - r2*rh) + 4*expl(B))*(4*expl(B + S) + 4*expl(B + A*r2 + S) - 2*(8 + B*(-1 + r2)*rh)*expl(A*r2 + T) + 8*expl(B + A*r2 + T))*expl(B + A*r2 + S + T)))/(32.*Power(-1 + r2,6)*Power(rh,6)) ;

dEqA[0] = 0 ;
dEqA[1] = 0 ;
dEqA[2] = 0 ;
dEqA[3] = 0 ;
dEqA[4] = -(Aa*Power(r,4)*(-1 + r2)*rh*expl(-(A*r2) - T))/2. + (Ba*(-1 + r2)*r2*rh*expl(-(A*r2) - T))/2. + (Ar*Br*Power(r,4)*Power(-1 + r2,4)*Power(rh,2)*expl(-(A*r2) - T))/8. - (Power(Ar,2)*Power(r,6)*Power(-1 + r2,4)*Power(rh,2)*expl(-(A*r2) - T))/16. - (Power(Br,2)*Power(-1 + r2,4)*r2*Power(rh,2)*expl(-(A*r2) - T))/16. + (Br*r*Power(-1 + r2,3)*(2 + A*Power(r,4) - (-4 + A)*r2)*Power(rh,2)*expl(-(A*r2) - T))/4. - (Ar*Power(r,3)*Power(-1 + r2,3)*(2 + A*Power(r,4) - (-4 + A)*r2)*Power(rh,2)*expl(-(A*r2) - T))/4. - ((-1 + r2)*r2*rh*Sa*expl(-(A*r2) - T))/2. + (3*Ar*Power(r,4)*Power(-1 + r2,4)*Power(rh,2)*Sr*expl(-(A*r2) - T))/8. - (3*Br*Power(-1 + r2,4)*r2*Power(rh,2)*Sr*expl(-(A*r2) - T))/8. + (3*r*Power(-1 + r2,3)*(2 + A*Power(r,4) - (-4 + A)*r2)*Power(rh,2)*Sr*expl(-(A*r2) - T))/4. - (9*Power(-1 + r2,4)*r2*Power(rh,2)*Power(Sr,2)*expl(-(A*r2) - T))/16. - (-1 + r2)*r2*rh*Ta*expl(-(A*r2) - T) - (Power(Fa,2)*Power(r,4)*Power(-1 + r2,2)*Power(rh,2)*expl(-2*B - A*r2 - 2*S - T))/64. - (Br*Fa*Power(r,3)*Power(-1 + r2,3)*Power(rh,2)*expl(-B - A*r2 - S - T))/16. + (Ar*Fa*Power(r,5)*Power(-1 + r2,3)*Power(rh,2)*expl(-B - A*r2 - S - T))/16. + (Fa*Power(-1 + r2,2)*r2*(2 + A*Power(r,4) - (-4 + A)*r2)*Power(rh,2)*expl(-B - A*r2 - S - T))/8. - (3*Fa*Power(r,3)*Power(-1 + r2,3)*Power(rh,2)*Sr*expl(-B - A*r2 - S - T))/16. + (Power(rh - r2*rh,6)*(16*r2*expl(B + A*r2 + S)*(2*expl(B + S) + (-8 + B*(rh - r2*rh))*expl(T) + 4*expl(B + T))*(4*expl(B + S) + 4*expl(B + A*r2 + S) - 2*(8 + B*(-1 + r2)*rh)*expl(A*r2 + T) + 8*expl(B + A*r2 + T)) + (32*Power(-1 + r2,2)*Power(rh,2)*expl(A*r2 + T)*((-Power(r,4) + r2)*expl(B + S) + (-Power(r,4) + r2)*expl(B + A*r2 + S) + 2*Power(r,4)*expl(A*r2 + T) + 4*Power(r,4)*expl(B + A*r2 + T))*((-1 - A*Power(r,4) + (-5 + A)*r2)*expl(B + S) - (-1 + r2)*expl(B + A*r2 + S) + 2*r2*expl(A*r2 + T) + 4*r2*expl(B + A*r2 + T)))/Power(r,2) + 16*Power(-1 + r2,2)*Power(rh,2)*expl(A*r2 + T)*Power((-1 - A*Power(r,4) + (-5 + A)*r2)*expl(B + S) - (-1 + r2)*expl(B + A*r2 + S) + 2*r2*expl(A*r2 + T) + 4*r2*expl(B + A*r2 + T),2) - 2*r2*(4*expl(B + S)*Power(4*expl(B + S) + 4*expl(B + A*r2 + S) - 2*(8 + B*(-1 + r2)*rh)*expl(A*r2 + T) + 8*expl(B + A*r2 + T),2) + (16*Power(-1 + r2,2)*Power(rh,2)*expl(A*r2 + T)*Power((-1 - A*Power(r,4) + (-5 + A)*r2)*expl(B + S) - (-1 + r2)*expl(B + A*r2 + S) + 2*r2*expl(A*r2 + T) + 4*r2*expl(B + A*r2 + T),2))/Power(r,2)))*expl(-2*(B + A*r2 + S + T)))/(64.*Power(-1 + r2,6)*Power(rh,6)) ;

dEqB[0] = (-2*Sqrt(expl(B + S)))/(rh - r2*rh) ;
dEqB[1] = (-2*Sqrt(expl(B + S)))/(rh - r2*rh) ;
dEqB[2] = 0 ;
dEqB[3] = 0 ;
dEqB[4] = -(Power(Ba,2)*Power(-1 + r2,2)*Power(rh,2)*expl(-B - S))/16. - (Power(Aa,2)*Power(r,4)*Power(-1 + r2,2)*Power(rh,2)*expl(-B - S))/16. - (Br*r*Power(-1 + r2,3)*Power(rh,2)*expl(-B - S))/2. + (Ar*Power(r,3)*Power(-1 + r2,3)*Power(rh,2)*expl(-B - S))/2. + (Aa*Ba*Power(-1 + r2,2)*r2*Power(rh,2)*expl(-B - S))/8. - (Ba*(-1 + r2)*rh*(8 + (-1 + B)*(-1 + r2)*rh)*expl(-B - S))/4. + (Aa*(-1 + r2)*r2*rh*(8 + (-1 + B)*(-1 + r2)*rh)*expl(-B - S))/4. + (Ba*Power(-1 + r2,2)*Power(rh,2)*Sa*expl(-B - S))/8. - (Aa*Power(-1 + r2,2)*r2*Power(rh,2)*Sa*expl(-B - S))/8. + ((-1 + r2)*rh*(8 + (-1 + B)*(-1 + r2)*rh)*Sa*expl(-B - S))/4. - (Power(-1 + r2,2)*Power(rh,2)*Power(Sa,2)*expl(-B - S))/16. - (3*r*Power(-1 + r2,3)*Power(rh,2)*Sr*expl(-B - S))/2. + (Ba*Power(-1 + r2,2)*Power(rh,2)*Ta*expl(-B - S))/4. - (Aa*Power(-1 + r2,2)*r2*Power(rh,2)*Ta*expl(-B - S))/4. + ((-1 + r2)*rh*(8 + (-1 + B)*(-1 + r2)*rh)*Ta*expl(-B - S))/2. - (Power(-1 + r2,2)*Power(rh,2)*Sa*Ta*expl(-B - S))/4. - (Power(-1 + r2,2)*Power(rh,2)*Power(Ta,2)*expl(-B - S))/4. - (Power(Fa,2)*Power(-1 + r2,2)*r2*Power(rh,2)*expl(-2*B - A*r2 - 2*S - T))/32. - (Br*Fa*r*Power(-1 + r2,3)*Power(rh,2)*expl(-B - A*r2 - S - T))/16. + (Ar*Fa*Power(r,3)*Power(-1 + r2,3)*Power(rh,2)*expl(-B - A*r2 - S - T))/16. - (3*Fa*r*Power(-1 + r2,3)*Power(rh,2)*Sr*expl(-B - A*r2 - S - T))/16. + (Fa*Power(-1 + r2,2)*Power(rh,2)*expl(-2*B - A*r2 - 2*S - T)*((1 + A*Power(r,4) - (-5 + A)*r2)*expl(B + S) + (-1 + r2)*expl(B + A*r2 + S) - 4*r2*expl(A*r2 + T) - 4*r2*expl(B + A*r2 + T)))/8. + (expl(-2*(B + A*r2 + S + T))*(4*expl(3*(B + S)) + 8*expl(3*B + A*r2 + 3*S) + 4*expl(3*B + 2*A*r2 + 3*S) - 8*Power(-1 + r2,2)*r2*Power(rh,2)*expl(3*(A*r2 + T)) - 4*(-1 + r2)*rh*expl(2*B + A*r2 + 2*S + T) + 16*expl(3*B + A*r2 + 2*S + T) - 4*(-1 + r2)*rh*expl(2*B + 2*A*r2 + 2*S + T) + 16*expl(3*B + 2*A*r2 + 2*S + T) + (-64 - 16*(-1 + B)*(-1 + r2)*rh + Power(-1 + r2,2)*(2*B - Power(B,2) + 4*(1 + A*Power(r,4) - (-5 + A)*r2))*Power(rh,2))*expl(B + 2*A*r2 + S + 2*T) - 8*(-1 + r2)*rh*expl(2*B + 2*A*r2 + S + 2*T) + 16*expl(3*B + 2*A*r2 + S + 2*T) + 4*Power(-1 + r2,3)*Power(rh,2)*expl(B + 3*A*r2 + S + 2*T) - 16*Power(-1 + r2,2)*r2*Power(rh,2)*expl(B + 3*A*r2 + 3*T)))/4. ;

dEqF[0] = 0 ;
dEqF[1] = 0 ;
dEqF[2] = 0 ;
dEqF[3] = 1 ;
dEqF[4] = 0 ;

dEqS[0] = (-2*Sqrt(expl(B + S)))/(rh - r2*rh) ;
dEqS[1] = (-2*Sqrt(expl(B + S)))/(rh - r2*rh) ;
dEqS[2] = 0 ;
dEqS[3] = 0 ;
dEqS[4] = -(Power(Ba,2)*Power(-1 + r2,2)*Power(rh,2)*expl(-B - S))/16. - (Power(Aa,2)*Power(r,4)*Power(-1 + r2,2)*Power(rh,2)*expl(-B - S))/16. + (Aa*Ba*Power(-1 + r2,2)*r2*Power(rh,2)*expl(-B - S))/8. + (Ba*Power(-1 + r2,2)*Power(rh,2)*Sa*expl(-B - S))/8. - (Aa*Power(-1 + r2,2)*r2*Power(rh,2)*Sa*expl(-B - S))/8. - (Power(-1 + r2,2)*Power(rh,2)*Power(Sa,2)*expl(-B - S))/16. + (Ba*Power(-1 + r2,2)*Power(rh,2)*Ta*expl(-B - S))/4. - (Aa*Power(-1 + r2,2)*r2*Power(rh,2)*Ta*expl(-B - S))/4. - (Power(-1 + r2,2)*Power(rh,2)*Sa*Ta*expl(-B - S))/4. - (Power(-1 + r2,2)*Power(rh,2)*Power(Ta,2)*expl(-B - S))/4. + (Aa*(-1 + r2)*r2*rh*(8 + B*(-1 + r2)*rh - 4*expl(B))*expl(-B - S))/4. + ((-1 + r2)*rh*Sa*(8 + B*(-1 + r2)*rh - 4*expl(B))*expl(-B - S))/4. + ((-1 + r2)*rh*Ta*(8 + B*(-1 + r2)*rh - 4*expl(B))*expl(-B - S))/2. - (Br*r*Power(-1 + r2,3)*Power(rh,2)*(1 + 2*expl(B))*expl(-B - S))/2. + (Ar*Power(r,3)*Power(-1 + r2,3)*Power(rh,2)*(1 + 2*expl(B))*expl(-B - S))/2. - (3*r*Power(-1 + r2,3)*Power(rh,2)*Sr*(1 + 2*expl(B))*expl(-B - S))/2. + (Ba*(-1 + r2)*rh*(-8 + B*(rh - r2*rh) + 4*expl(B))*expl(-B - S))/4. - (Power(Fa,2)*Power(-1 + r2,2)*r2*Power(rh,2)*expl(-2*B - A*r2 - 2*S - T))/32. - (Br*Fa*r*Power(-1 + r2,3)*Power(rh,2)*expl(-B - A*r2 - S - T))/16. + (Ar*Fa*Power(r,3)*Power(-1 + r2,3)*Power(rh,2)*expl(-B - A*r2 - S - T))/16. - (3*Fa*r*Power(-1 + r2,3)*Power(rh,2)*Sr*expl(-B - A*r2 - S - T))/16. + (Fa*Power(-1 + r2,2)*Power(rh,2)*expl(-2*B - A*r2 - 2*S - T)*((1 + A*Power(r,4) - (-5 + A)*r2)*expl(B + S) + (-1 + r2)*expl(B + A*r2 + S) - 4*r2*expl(A*r2 + T) - 8*r2*expl(B + A*r2 + T)))/8. + (expl(-2*(B + A*r2 + S + T))*(4*expl(3*(B + S)) + 8*expl(3*B + A*r2 + 3*S) + 4*expl(3*B + 2*A*r2 + 3*S) - 8*Power(-1 + r2,2)*r2*Power(rh,2)*expl(3*(A*r2 + T)) - (16*B*(-1 + r2)*rh + Power(B,2)*Power(-1 + r2,2)*Power(rh,2) - 4*(-16 + Power(-1 + r2,2)*(1 + A*Power(r,4) - (-5 + A)*r2)*Power(rh,2)))*expl(B + 2*A*r2 + S + 2*T) + 8*(8 + B*(-1 + r2)*rh + Power(-1 + r2,2)*(1 + A*Power(r,4) - (-5 + A)*r2)*Power(rh,2))*expl(2*B + 2*A*r2 + S + 2*T) - 16*expl(3*B + 2*A*r2 + S + 2*T) + 4*Power(-1 + r2,3)*Power(rh,2)*expl(B + 3*A*r2 + S + 2*T) + 8*Power(-1 + r2,3)*Power(rh,2)*expl(2*B + 3*A*r2 + S + 2*T) - 32*Power(-1 + r2,2)*r2*Power(rh,2)*expl(B + 3*A*r2 + 3*T) - 32*Power(-1 + r2,2)*r2*Power(rh,2)*expl(2*B + 3*A*r2 + 3*T)))/4. ;