* * epsilon-expansion of hypergometric functions 3F2 with half-integer value of parameter. * * The functions F32, and any two of the (F32A,F32B,F32C) are basis for differential reduction * of 3F2[I1+a1*ep,I2+a2*ep,I3+a3*ep; 1/2+I4+f1*ep,I5+b1*ep, z] * **************************************************************************************** * Modifications: * 15.04.2009: presenation for F32 is improved; * 25.03.2009: functions "num" and "den" are added; * **************************************************************************************** * * F32(a1?,a2?,a3?,f1?,c1?,z?,ep?) = Hypergeometric3F2[a1*ep,a2*ep, a3*ep; 1/2+f1*ep, 1+c1*ep;z] * * F32A(a1?,a2?,a3?,f1?,c1?,z?,ep?) = Hypergeometric3F2[1+a1*ep,1+a2*ep, 1+a3*ep; 3/2+f1*ep, 2+c1*ep;z] * * F32B(a1?,a2?,a3?,f1?,e1?,z?,ep?) = Hypergeometri3F2[1+a1*ep,1+a2*ep, 1+a3*ep; 3/2+f1*ep, 1+e1*ep;z] * * F32C(a1?,a2?,d1?,f1?,c1?,z?,ep?) = Hypergeometri3F2[1+a1*ep,1+a2*ep, 2+d1*ep; 3/2+f1*ep, 2+c1*ep;z] * ***************************************************************************************** * * "z" is argument of hypergeometric function * "ep" is a small parametr * * "y" is auxiliary variable defined by as * y = (sqrt(z-1)-sqrt(z) )/(sqrt(z-1) + sqrt(z) ) = [1 - \sqrt(z/(z-1))]/[1 + \sqrt(z/(z-1))] * and inverse relation is * z = -(1-y)^2/4y; * [1-y] = 1-y; * [1+y] = 1+y; * den(a,b,ep,j)=1/(a+b*ep)^j * num(a,b,ep,j)=(a+b*ep)^j * "zt" is Zeta-function, "ln" is logarithm * "Li" and "S" are classical and Nielsen polylgarithms. * "H" is a harmonic polylogarithm H(-1,0,0,1,z) = \int_0^z dt \frac{Li(3,t)}{1+t}, * See Eq.(2.47) in Davydychev,Kalmykov, Nucl.Phys. B699 (2004) 3. **************************************************************************************** *s a1,a2,a3,b1,c1,d1,f1,e1,ep,z,y,[1-y],[1+y],[1+d1*ep],[1+c1*ep],[1+2*f1*ep]; *c Li,S,zt,ln,H; *c F32, F32A, F32B, F32C, F43A; *c den,num; * * Example of calling: * g TEST1 = F32(1,2,3,2,5,z,ep); * g TEST2 = F32A(1,2,3,2,5,z,ep); * g TEST3 = F32B(3,2,3,2,5,z,ep); * g TEST4 = F32C(14,3,3,2,5,z,ep); * ***************************************************************************************** * F32 = Hypergeometric3F2[a1*ep,a2*ep, a3*ep; 1/2+f1*ep, 1+c1*ep;z] **************************************************************************************** id F32(a1?,a2?,a3?,f1?,c1?,z?,ep?) = 1 + a1*a2*a3*ep^3*den(1,c1,ep,1)*2*den(1,2*f1,ep,1)*( -[1-y]^2/(4*y) )*F43A(0,a1,a2,a3,f1,0,c1,z,ep); .sort ***************************************************************************************** * F32A = Hypergeometric3F2[1+a1*ep,1+a2*ep, 1+a3*ep; 3/2+f1*ep, 2+c1*ep;z] ***************************************************************************************** id F32A(a1?,a2?,a3?,f1?,c1?,z?,ep?) = (-2*y/[1-y]^2)*num(1,c1,ep,1)*num(1,2*f1,ep,1)*( - 1/2*ln(y)*ln(y) + ep * ( ln(y)*ln(y)*ln(1 - y)*c1 - 1/6*ln(y)*ln(y)*ln(y)*a1 - 1/6*ln(y )*ln(y)*ln(y)*a2 - 1/6*ln(y)*ln(y)*ln(y)*a3 + ln(y)*zt(2)*a1 + ln(y)* zt(2)*a2 + ln(y)*zt(2)*a3 - ln(y)*zt(2)*c1 - 3*ln(y)*zt(2)*f1 - 2*ln( y)*Li(2, - y)*a1 - 2*ln(y)*Li(2, - y)*a2 - 2*ln(y)*Li(2, - y)*a3 + 2* ln(y)*Li(2, - y)*c1 + 2*ln(y)*Li(2, - y)*f1 + 2*ln(y)*Li(2,y)*c1 - 2* ln(y)*Li(2,y)*f1 + 3*zt(3)*a1 + 3*zt(3)*a2 + 3*zt(3)*a3 - zt(3)*c1 - 7*zt(3)*f1 + 4*Li(3, - y)*a1 + 4*Li(3, - y)*a2 + 4*Li(3, - y)*a3 - 4* Li(3, - y)*c1 - 4*Li(3, - y)*f1 - 2*Li(3,y)*c1 + 4*Li(3,y)*f1 ) + ep^2 * ( - 8*S(1,2, - y)*ln(y)*a1*a2 - 8*S(1,2, - y)*ln(y)*a1*a3 + 4 *S(1,2, - y)*ln(y)*a1*c1 + 12*S(1,2, - y)*ln(y)*a1*f1 - 4*S(1,2, - y) *ln(y)*a1^2 - 8*S(1,2, - y)*ln(y)*a2*a3 + 4*S(1,2, - y)*ln(y)*a2*c1 + 12*S(1,2, - y)*ln(y)*a2*f1 - 4*S(1,2, - y)*ln(y)*a2^2 + 4*S(1,2, - y)*ln(y)*a3*c1 + 12*S(1,2, - y)*ln(y)*a3*f1 - 4*S(1,2, - y)*ln(y)* a3^2 - 8*S(1,2, - y)*ln(y)*c1*f1 - 8*S(1,2, - y)*ln(y)*f1^2 + 2*S(1,2 ,y^2)*ln(y)*a1*c1 - 2*S(1,2,y^2)*ln(y)*a1*f1 + 2*S(1,2,y^2)*ln(y)*a2* c1 - 2*S(1,2,y^2)*ln(y)*a2*f1 + 2*S(1,2,y^2)*ln(y)*a3*c1 - 2*S(1,2, y^2)*ln(y)*a3*f1 - 2*S(1,2,y^2)*ln(y)*c1^2 + 2*S(1,2,y^2)*ln(y)*f1^2 - 4*S(1,2,y)*ln(y)*a1*c1 + 4*S(1,2,y)*ln(y)*a1*f1 - 4*S(1,2,y)*ln(y) *a2*c1 + 4*S(1,2,y)*ln(y)*a2*f1 - 4*S(1,2,y)*ln(y)*a3*c1 + 4*S(1,2,y) *ln(y)*a3*f1 + 8*S(1,2,y)*ln(y)*c1*f1 - 8*S(1,2,y)*ln(y)*f1^2 + 4*S(2 ,2, - y)*a1*c1 + 4*S(2,2, - y)*a2*c1 + 4*S(2,2, - y)*a3*c1 - 4*S(2,2, - y)*c1*f1 - 4*S(2,2, - y)*c1^2 - S(2,2,y^2)*a1*c1 - S(2,2,y^2)*a2* c1 - S(2,2,y^2)*a3*c1 + S(2,2,y^2)*c1*f1 + S(2,2,y^2)*c1^2 + 4*S(2,2, y)*a1*c1 + 4*S(2,2,y)*a2*c1 + 4*S(2,2,y)*a3*c1 - 12*S(2,2,y)*c1*f1 - 4*H( - y)*a1*c1 - 4*H( - y)*a2*c1 - 4*H( - y)*a3*c1 + 4*H( - y)*c1*f1 + 4*H( - y)*c1^2 - 6*ln(1 - y)*zt(3)*a1*c1 - 6*ln(1 - y)*zt(3)*a2*c1 - 6*ln(1 - y)*zt(3)*a3*c1 + 14*ln(1 - y)*zt(3)*c1*f1 + 2*ln(1 - y)* zt(3)*c1^2 - 4*ln(1 - y)*Li(3, - y)*a1*c1 - 4*ln(1 - y)*Li(3, - y)*a2 *c1 - 4*ln(1 - y)*Li(3, - y)*a3*c1 + 4*ln(1 - y)*Li(3, - y)*c1*f1 + 4 *ln(1 - y)*Li(3, - y)*c1^2 - 8*ln(1 - y)*Li(3,y)*c1*f1 + 4*ln(1 - y)* Li(3,y)*c1^2 - 2*ln(y)*ln(1 - y)*zt(2)*a1*c1 - 2*ln(y)*ln(1 - y)*zt(2 )*a2*c1 - 2*ln(y)*ln(1 - y)*zt(2)*a3*c1 + 6*ln(y)*ln(1 - y)*zt(2)*c1* f1 + 2*ln(y)*ln(1 - y)*zt(2)*c1^2 + 4*ln(y)*ln(1 - y)*Li(2, - y)*a1* c1 + 4*ln(y)*ln(1 - y)*Li(2, - y)*a2*c1 + 4*ln(y)*ln(1 - y)*Li(2, - y )*a3*c1 - 4*ln(y)*ln(1 - y)*Li(2, - y)*c1*f1 - 4*ln(y)*ln(1 - y)*Li(2 , - y)*c1^2 + 4*ln(y)*ln(1 - y)*Li(2,y)*c1*f1 - 4*ln(y)*ln(1 - y)*Li( 2,y)*c1^2 - ln(y)*ln(y)*ln(1 - y)*ln(1 - y)*c1^2 + 1/3*ln(y)*ln(y)* ln(y)*ln(1 - y)*a1*c1 + 1/3*ln(y)*ln(y)*ln(y)*ln(1 - y)*a2*c1 + 1/3* ln(y)*ln(y)*ln(y)*ln(1 - y)*a3*c1 - 1/24*ln(y)*ln(y)*ln(y)*ln(y)*a1* a2 - 1/24*ln(y)*ln(y)*ln(y)*ln(y)*a1*a3 - 1/24*ln(y)*ln(y)*ln(y)*ln(y )*a1^2 - 1/24*ln(y)*ln(y)*ln(y)*ln(y)*a2*a3 - 1/24*ln(y)*ln(y)*ln(y)* ln(y)*a2^2 - 1/24*ln(y)*ln(y)*ln(y)*ln(y)*a3^2 + ln(y)*ln(y)*zt(2)*a1 *a2 + ln(y)*ln(y)*zt(2)*a1*a3 - 1/2*ln(y)*ln(y)*zt(2)*a1*c1 - 3/2*ln( y)*ln(y)*zt(2)*a1*f1 + 1/2*ln(y)*ln(y)*zt(2)*a1^2 + ln(y)*ln(y)*zt(2) *a2*a3 - 1/2*ln(y)*ln(y)*zt(2)*a2*c1 - 3/2*ln(y)*ln(y)*zt(2)*a2*f1 + 1/2*ln(y)*ln(y)*zt(2)*a2^2 - 1/2*ln(y)*ln(y)*zt(2)*a3*c1 - 3/2*ln(y)* ln(y)*zt(2)*a3*f1 + 1/2*ln(y)*ln(y)*zt(2)*a3^2 - 2*ln(y)*ln(y)*Li(2, - y)*a1*a2 - 2*ln(y)*ln(y)*Li(2, - y)*a1*a3 + 2*ln(y)*ln(y)*Li(2, - y)*a1*c1 + ln(y)*ln(y)*Li(2, - y)*a1*f1 - ln(y)*ln(y)*Li(2, - y)*a1^2 - 2*ln(y)*ln(y)*Li(2, - y)*a2*a3 + 2*ln(y)*ln(y)*Li(2, - y)*a2*c1 + ln(y)*ln(y)*Li(2, - y)*a2*f1 - ln(y)*ln(y)*Li(2, - y)*a2^2 + 2*ln(y)* ln(y)*Li(2, - y)*a3*c1 + ln(y)*ln(y)*Li(2, - y)*a3*f1 - ln(y)*ln(y)* Li(2, - y)*a3^2 - ln(y)*ln(y)*Li(2, - y)*c1*f1 - ln(y)*ln(y)*Li(2, - y)*c1^2 + ln(y)*ln(y)*Li(2,y)*a1*c1 - ln(y)*ln(y)*Li(2,y)*a1*f1 + ln( y)*ln(y)*Li(2,y)*a2*c1 - ln(y)*ln(y)*Li(2,y)*a2*f1 + ln(y)*ln(y)*Li(2 ,y)*a3*c1 - ln(y)*ln(y)*Li(2,y)*a3*f1 + ln(y)*ln(y)*Li(2,y)*c1*f1 - ln(y)*ln(y)*Li(2,y)*c1^2 + 4*ln(y)*zt(3)*a1*a2 + 4*ln(y)*zt(3)*a1*a3 - ln(y)*zt(3)*a1*c1 - 7*ln(y)*zt(3)*a1*f1 + 2*ln(y)*zt(3)*a1^2 + 4* ln(y)*zt(3)*a2*a3 - ln(y)*zt(3)*a2*c1 - 7*ln(y)*zt(3)*a2*f1 + 2*ln(y) *zt(3)*a2^2 - ln(y)*zt(3)*a3*c1 - 7*ln(y)*zt(3)*a3*f1 + 2*ln(y)*zt(3) *a3^2 + ln(y)*zt(3)*c1^2 + 7*ln(y)*zt(3)*f1^2 + 4*ln(y)*Li(3, - y)*a1 *a2 + 4*ln(y)*Li(3, - y)*a1*a3 - 6*ln(y)*Li(3, - y)*a1*c1 - 2*ln(y)* Li(3, - y)*a1*f1 + 2*ln(y)*Li(3, - y)*a1^2 + 4*ln(y)*Li(3, - y)*a2*a3 - 6*ln(y)*Li(3, - y)*a2*c1 - 2*ln(y)*Li(3, - y)*a2*f1 + 2*ln(y)*Li(3 , - y)*a2^2 - 6*ln(y)*Li(3, - y)*a3*c1 - 2*ln(y)*Li(3, - y)*a3*f1 + 2 *ln(y)*Li(3, - y)*a3^2 + 4*ln(y)*Li(3, - y)*c1*f1 + 4*ln(y)*Li(3, - y )*c1^2 - 2*ln(y)*Li(3,y)*a1*c1 + 2*ln(y)*Li(3,y)*a1*f1 - 2*ln(y)*Li(3 ,y)*a2*c1 + 2*ln(y)*Li(3,y)*a2*f1 - 2*ln(y)*Li(3,y)*a3*c1 + 2*ln(y)* Li(3,y)*a3*f1 - 4*ln(y)*Li(3,y)*c1*f1 + 4*ln(y)*Li(3,y)*c1^2 + 4*zt(2 )*Li(2, - y)*a1*a2 + 4*zt(2)*Li(2, - y)*a1*a3 - 4*zt(2)*Li(2, - y)*a1 *c1 - 8*zt(2)*Li(2, - y)*a1*f1 + 2*zt(2)*Li(2, - y)*a1^2 + 4*zt(2)* Li(2, - y)*a2*a3 - 4*zt(2)*Li(2, - y)*a2*c1 - 8*zt(2)*Li(2, - y)*a2* f1 + 2*zt(2)*Li(2, - y)*a2^2 - 4*zt(2)*Li(2, - y)*a3*c1 - 8*zt(2)*Li( 2, - y)*a3*f1 + 2*zt(2)*Li(2, - y)*a3^2 + 8*zt(2)*Li(2, - y)*c1*f1 + 2*zt(2)*Li(2, - y)*c1^2 + 6*zt(2)*Li(2, - y)*f1^2 - 2*zt(2)*Li(2,y)* a1*c1 + 2*zt(2)*Li(2,y)*a1*f1 - 2*zt(2)*Li(2,y)*a2*c1 + 2*zt(2)*Li(2, y)*a2*f1 - 2*zt(2)*Li(2,y)*a3*c1 + 2*zt(2)*Li(2,y)*a3*f1 + 4*zt(2)* Li(2,y)*c1*f1 + 2*zt(2)*Li(2,y)*c1^2 - 6*zt(2)*Li(2,y)*f1^2 + 5/2*zt( 4)*a1*a2 + 5/2*zt(4)*a1*a3 + 3/2*zt(4)*a1*c1 - 15/2*zt(4)*a1*f1 + 5/4 *zt(4)*a1^2 + 5/2*zt(4)*a2*a3 + 3/2*zt(4)*a2*c1 - 15/2*zt(4)*a2*f1 + 5/4*zt(4)*a2^2 + 3/2*zt(4)*a3*c1 - 15/2*zt(4)*a3*f1 + 5/4*zt(4)*a3^2 - 15/2*zt(4)*c1*f1 + 1/4*zt(4)*c1^2 + 45/4*zt(4)*f1^2 + 4*Li(2, - y) *Li(2, - y)*a1*a2 + 4*Li(2, - y)*Li(2, - y)*a1*a3 - 4*Li(2, - y)*Li(2 , - y)*a1*c1 - 4*Li(2, - y)*Li(2, - y)*a1*f1 + 2*Li(2, - y)*Li(2, - y )*a1^2 + 4*Li(2, - y)*Li(2, - y)*a2*a3 - 4*Li(2, - y)*Li(2, - y)*a2* c1 - 4*Li(2, - y)*Li(2, - y)*a2*f1 + 2*Li(2, - y)*Li(2, - y)*a2^2 - 4 *Li(2, - y)*Li(2, - y)*a3*c1 - 4*Li(2, - y)*Li(2, - y)*a3*f1 + 2*Li(2 , - y)*Li(2, - y)*a3^2 + 4*Li(2, - y)*Li(2, - y)*c1*f1 + 2*Li(2, - y) *Li(2, - y)*c1^2 + 2*Li(2, - y)*Li(2, - y)*f1^2 + 4*Li(2,y)*Li(2, - y )*a1*f1 + 4*Li(2,y)*Li(2, - y)*a2*f1 + 4*Li(2,y)*Li(2, - y)*a3*f1 - 4 *Li(2,y)*Li(2, - y)*c1*f1 - 4*Li(2,y)*Li(2, - y)*f1^2 - 2*Li(2,y)*Li( 2,y)*c1*f1 + 2*Li(2,y)*Li(2,y)*f1^2 + 6*Li(4, - y)*a1*c1 + 6*Li(4, - y)*a2*c1 + 6*Li(4, - y)*a3*c1 - 6*Li(4, - y)*c1*f1 - 6*Li(4, - y)* c1^2 + 2*Li(4,y)*a1*c1 + 2*Li(4,y)*a2*c1 + 2*Li(4,y)*a3*c1 + 6*Li(4,y )*c1*f1 - 6*Li(4,y)*c1^2 ) + ep^3*F32A(a1,a2,a3,f1,c1,z,ep^3) ); .sort **************************************************************************************** * F32B = Hypergeometri3F2[1+a1*ep,1+a2*ep, 1+a3*ep; 3/2+f1*ep, 1+e1*ep;z] ***************************************************************************************** id F32B(a1?,a2?,a3?,f1?,e1?,z?,ep?) = (-2*y/[1-y]^2)*num(1,2*f1,ep,1)*[1-y]/[1+y]* ( + ln(y) + ep * ( - 2*ln(y)*ln(1 - y)*f1 + 2*ln(y)*ln(1 + y)*f1 + 2*ln(y )*ln(1 + y)*e1 - 2*ln(y)*ln(1 + y)*a3 - 2*ln(y)*ln(1 + y)*a2 - 2*ln(y)*ln(1 + y)*a1 - 1/2*ln(y)^2*e1 + 1/2*ln(y)^2*a3 + 1/ 2*ln(y)^2*a2 + 1/2*ln(y)^2*a1 + 3*zt(2)*f1 + zt(2)*e1 - zt(2)* a3 - zt(2)*a2 - zt(2)*a1 + 2*Li(2, - y)*f1 + 2*Li(2, - y)*e1 - 2*Li(2, - y)*a3 - 2*Li(2, - y)*a2 - 2*Li(2, - y)*a1 - 2*Li( 2,y)*f1 ) + ep^2 * ( 8*S(1,2, - y)*f1^2 + 12*S(1,2, - y)*e1*f1 + 4*S(1,2, - y)*e1^2 - 12*S(1,2, - y)*a3*f1 - 8*S(1,2, - y)*a3*e1 + 4*S( 1,2, - y)*a3^2 - 12*S(1,2, - y)*a2*f1 - 8*S(1,2, - y)*a2*e1 + 8*S(1,2, - y)*a2*a3 + 4*S(1,2, - y)*a2^2 - 12*S(1,2, - y)*a1* f1 - 8*S(1,2, - y)*a1*e1 + 8*S(1,2, - y)*a1*a3 + 8*S(1,2, - y) *a1*a2 + 4*S(1,2, - y)*a1^2 - 2*S(1,2,y^2)*f1^2 - 2*S(1,2,y^2) *e1*f1 + 2*S(1,2,y^2)*a3*f1 + 2*S(1,2,y^2)*a2*f1 + 2*S(1,2,y^2 )*a1*f1 + 8*S(1,2,y)*f1^2 + 4*S(1,2,y)*e1*f1 - 4*S(1,2,y)*a3* f1 - 4*S(1,2,y)*a2*f1 - 4*S(1,2,y)*a1*f1 - 6*ln(1 - y)*zt(2)* f1^2 - 2*ln(1 - y)*zt(2)*e1*f1 + 2*ln(1 - y)*zt(2)*a3*f1 + 2* ln(1 - y)*zt(2)*a2*f1 + 2*ln(1 - y)*zt(2)*a1*f1 - 4*ln(1 - y)* Li(2, - y)*f1^2 - 4*ln(1 - y)*Li(2, - y)*e1*f1 + 4*ln(1 - y)* Li(2, - y)*a3*f1 + 4*ln(1 - y)*Li(2, - y)*a2*f1 + 4*ln(1 - y)* Li(2, - y)*a1*f1 + 4*ln(1 - y)*Li(2,y)*f1^2 + 6*ln(1 + y)*zt(2 )*f1^2 + 8*ln(1 + y)*zt(2)*e1*f1 + 2*ln(1 + y)*zt(2)*e1^2 - 8* ln(1 + y)*zt(2)*a3*f1 - 4*ln(1 + y)*zt(2)*a3*e1 + 2*ln(1 + y)* zt(2)*a3^2 - 8*ln(1 + y)*zt(2)*a2*f1 - 4*ln(1 + y)*zt(2)*a2*e1 + 4*ln(1 + y)*zt(2)*a2*a3 + 2*ln(1 + y)*zt(2)*a2^2 - 8*ln(1 + y)*zt(2)*a1*f1 - 4*ln(1 + y)*zt(2)*a1*e1 + 4*ln(1 + y)*zt(2 )*a1*a3 + 4*ln(1 + y)*zt(2)*a1*a2 + 2*ln(1 + y)*zt(2)*a1^2 + 4 *ln(1 + y)*Li(2, - y)*f1^2 + 8*ln(1 + y)*Li(2, - y)*e1*f1 + 4* ln(1 + y)*Li(2, - y)*e1^2 - 8*ln(1 + y)*Li(2, - y)*a3*f1 - 8* ln(1 + y)*Li(2, - y)*a3*e1 + 4*ln(1 + y)*Li(2, - y)*a3^2 - 8* ln(1 + y)*Li(2, - y)*a2*f1 - 8*ln(1 + y)*Li(2, - y)*a2*e1 + 8* ln(1 + y)*Li(2, - y)*a2*a3 + 4*ln(1 + y)*Li(2, - y)*a2^2 - 8* ln(1 + y)*Li(2, - y)*a1*f1 - 8*ln(1 + y)*Li(2, - y)*a1*e1 + 8* ln(1 + y)*Li(2, - y)*a1*a3 + 8*ln(1 + y)*Li(2, - y)*a1*a2 + 4* ln(1 + y)*Li(2, - y)*a1^2 - 4*ln(1 + y)*Li(2,y)*f1^2 - 4*ln(1 + y)*Li(2,y)*e1*f1 + 4*ln(1 + y)*Li(2,y)*a3*f1 + 4*ln(1 + y)* Li(2,y)*a2*f1 + 4*ln(1 + y)*Li(2,y)*a1*f1 + 2*ln(y)*ln(1 - y)^ 2*f1^2 - 4*ln(y)*ln(1 - y)*ln(1 + y)*f1^2 - 4*ln(y)*ln(1 - y)* ln(1 + y)*e1*f1 + 4*ln(y)*ln(1 - y)*ln(1 + y)*a3*f1 + 4*ln(y)* ln(1 - y)*ln(1 + y)*a2*f1 + 4*ln(y)*ln(1 - y)*ln(1 + y)*a1*f1 + 2*ln(y)*ln(1 + y)^2*f1^2 + 4*ln(y)*ln(1 + y)^2*e1*f1 + 2* ln(y)*ln(1 + y)^2*e1^2 - 4*ln(y)*ln(1 + y)^2*a3*f1 - 4*ln(y)* ln(1 + y)^2*a3*e1 + 2*ln(y)*ln(1 + y)^2*a3^2 - 4*ln(y)*ln(1 + y)^2*a2*f1 - 4*ln(y)*ln(1 + y)^2*a2*e1 + 4*ln(y)*ln(1 + y)^2* a2*a3 + 2*ln(y)*ln(1 + y)^2*a2^2 - 4*ln(y)*ln(1 + y)^2*a1*f1 - 4*ln(y)*ln(1 + y)^2*a1*e1 + 4*ln(y)*ln(1 + y)^2*a1*a3 + 4* ln(y)*ln(1 + y)^2*a1*a2 + 2*ln(y)*ln(1 + y)^2*a1^2 + ln(y)^2* ln(1 - y)*e1*f1 - ln(y)^2*ln(1 - y)*a3*f1 - ln(y)^2*ln(1 - y)* a2*f1 - ln(y)^2*ln(1 - y)*a1*f1 - ln(y)^2*ln(1 + y)*e1*f1 - ln(y)^2*ln(1 + y)*e1^2 + ln(y)^2*ln(1 + y)*a3*f1 + 2*ln(y)^2* ln(1 + y)*a3*e1 - ln(y)^2*ln(1 + y)*a3^2 + ln(y)^2*ln(1 + y)* a2*f1 + 2*ln(y)^2*ln(1 + y)*a2*e1 - 2*ln(y)^2*ln(1 + y)*a2*a3 - ln(y)^2*ln(1 + y)*a2^2 + ln(y)^2*ln(1 + y)*a1*f1 + 2*ln(y)^ 2*ln(1 + y)*a1*e1 - 2*ln(y)^2*ln(1 + y)*a1*a3 - 2*ln(y)^2*ln(1 + y)*a1*a2 - ln(y)^2*ln(1 + y)*a1^2 - 1/6*ln(y)^3*a3*e1 + 1/6 *ln(y)^3*a3^2 - 1/6*ln(y)^3*a2*e1 + 1/6*ln(y)^3*a2*a3 + 1/6* ln(y)^3*a2^2 - 1/6*ln(y)^3*a1*e1 + 1/6*ln(y)^3*a1*a3 + 1/6*ln( y)^3*a1*a2 + 1/6*ln(y)^3*a1^2 - 3*ln(y)*zt(2)*e1*f1 - ln(y)* zt(2)*e1^2 + 3*ln(y)*zt(2)*a3*f1 + 2*ln(y)*zt(2)*a3*e1 - ln(y) *zt(2)*a3^2 + 3*ln(y)*zt(2)*a2*f1 + 2*ln(y)*zt(2)*a2*e1 - 2* ln(y)*zt(2)*a2*a3 - ln(y)*zt(2)*a2^2 + 3*ln(y)*zt(2)*a1*f1 + 2 *ln(y)*zt(2)*a1*e1 - 2*ln(y)*zt(2)*a1*a3 - 2*ln(y)*zt(2)*a1*a2 - ln(y)*zt(2)*a1^2 - 7*zt(3)*f1^2 - 7*zt(3)*e1*f1 - 2*zt(3)* e1^2 + 7*zt(3)*a3*f1 + 4*zt(3)*a3*e1 - 2*zt(3)*a3^2 + 7*zt(3)* a2*f1 + 4*zt(3)*a2*e1 - 4*zt(3)*a2*a3 - 2*zt(3)*a2^2 + 7*zt(3) *a1*f1 + 4*zt(3)*a1*e1 - 4*zt(3)*a1*a3 - 4*zt(3)*a1*a2 - 2*zt( 3)*a1^2 - 2*Li(3, - y)*e1*f1 - 2*Li(3, - y)*e1^2 + 2*Li(3, - y )*a3*f1 + 4*Li(3, - y)*a3*e1 - 2*Li(3, - y)*a3^2 + 2*Li(3, - y )*a2*f1 + 4*Li(3, - y)*a2*e1 - 4*Li(3, - y)*a2*a3 - 2*Li(3, - y)*a2^2 + 2*Li(3, - y)*a1*f1 + 4*Li(3, - y)*a1*e1 - 4*Li(3, - y)*a1*a3 - 4*Li(3, - y)*a1*a2 - 2*Li(3, - y)*a1^2 + 2*Li(3,y)* e1*f1 - 2*Li(3,y)*a3*f1 - 2*Li(3,y)*a2*f1 - 2*Li(3,y)*a1*f1 ) + ep^3 * ( 16*S(1,2, - y)*ln(1 + y)*e1*f1^2 + 24*S(1,2, - y)*ln( 1 + y)*e1^2*f1 + 8*S(1,2, - y)*ln(1 + y)*e1^3 - 16*S(1,2, - y) *ln(1 + y)*a3*f1^2 - 48*S(1,2, - y)*ln(1 + y)*a3*e1*f1 - 24*S( 1,2, - y)*ln(1 + y)*a3*e1^2 + 24*S(1,2, - y)*ln(1 + y)*a3^2*f1 + 24*S(1,2, - y)*ln(1 + y)*a3^2*e1 - 8*S(1,2, - y)*ln(1 + y)* a3^3 - 16*S(1,2, - y)*ln(1 + y)*a2*f1^2 - 48*S(1,2, - y)*ln(1 + y)*a2*e1*f1 - 24*S(1,2, - y)*ln(1 + y)*a2*e1^2 + 48*S(1,2, - y)*ln(1 + y)*a2*a3*f1 + 48*S(1,2, - y)*ln(1 + y)*a2*a3*e1 - 24*S(1,2, - y)*ln(1 + y)*a2*a3^2 + 24*S(1,2, - y)*ln(1 + y) *a2^2*f1 + 24*S(1,2, - y)*ln(1 + y)*a2^2*e1 - 24*S(1,2, - y)* ln(1 + y)*a2^2*a3 - 8*S(1,2, - y)*ln(1 + y)*a2^3 - 16*S(1,2, - y)*ln(1 + y)*a1*f1^2 - 48*S(1,2, - y)*ln(1 + y)*a1*e1*f1 - 24*S(1,2, - y)*ln(1 + y)*a1*e1^2 + 48*S(1,2, - y)*ln(1 + y)*a1 *a3*f1 + 48*S(1,2, - y)*ln(1 + y)*a1*a3*e1 - 24*S(1,2, - y)* ln(1 + y)*a1*a3^2 + 48*S(1,2, - y)*ln(1 + y)*a1*a2*f1 + 48*S(1 ,2, - y)*ln(1 + y)*a1*a2*e1 - 48*S(1,2, - y)*ln(1 + y)*a1*a2* a3 - 24*S(1,2, - y)*ln(1 + y)*a1*a2^2 + 24*S(1,2, - y)*ln(1 + y)*a1^2*f1 + 24*S(1,2, - y)*ln(1 + y)*a1^2*e1 - 24*S(1,2, - y) *ln(1 + y)*a1^2*a3 - 24*S(1,2, - y)*ln(1 + y)*a1^2*a2 - 8*S(1, 2, - y)*ln(1 + y)*a1^3 - 2*S(1,2,y^2)*ln(1 - y)*e1*f1^2 - 2*S( 1,2,y^2)*ln(1 - y)*e1^2*f1 + 2*S(1,2,y^2)*ln(1 - y)*a3*f1^2 + 4*S(1,2,y^2)*ln(1 - y)*a3*e1*f1 - 2*S(1,2,y^2)*ln(1 - y)*a3^2* f1 + 2*S(1,2,y^2)*ln(1 - y)*a2*f1^2 + 4*S(1,2,y^2)*ln(1 - y)* a2*e1*f1 - 4*S(1,2,y^2)*ln(1 - y)*a2*a3*f1 - 2*S(1,2,y^2)*ln(1 - y)*a2^2*f1 + 2*S(1,2,y^2)*ln(1 - y)*a1*f1^2 + 4*S(1,2,y^2)* ln(1 - y)*a1*e1*f1 - 4*S(1,2,y^2)*ln(1 - y)*a1*a3*f1 - 4*S(1,2 ,y^2)*ln(1 - y)*a1*a2*f1 - 2*S(1,2,y^2)*ln(1 - y)*a1^2*f1 - 2* S(1,2,y^2)*ln(1 + y)*e1*f1^2 - 2*S(1,2,y^2)*ln(1 + y)*e1^2*f1 + 2*S(1,2,y^2)*ln(1 + y)*a3*f1^2 + 4*S(1,2,y^2)*ln(1 + y)*a3* e1*f1 - 2*S(1,2,y^2)*ln(1 + y)*a3^2*f1 + 2*S(1,2,y^2)*ln(1 + y )*a2*f1^2 + 4*S(1,2,y^2)*ln(1 + y)*a2*e1*f1 - 4*S(1,2,y^2)*ln( 1 + y)*a2*a3*f1 - 2*S(1,2,y^2)*ln(1 + y)*a2^2*f1 + 2*S(1,2,y^2 )*ln(1 + y)*a1*f1^2 + 4*S(1,2,y^2)*ln(1 + y)*a1*e1*f1 - 4*S(1, 2,y^2)*ln(1 + y)*a1*a3*f1 - 4*S(1,2,y^2)*ln(1 + y)*a1*a2*f1 - 2*S(1,2,y^2)*ln(1 + y)*a1^2*f1 + 16*S(1,2,y)*ln(1 - y)*e1*f1^2 + 8*S(1,2,y)*ln(1 - y)*e1^2*f1 - 16*S(1,2,y)*ln(1 - y)*a3* f1^2 - 16*S(1,2,y)*ln(1 - y)*a3*e1*f1 + 8*S(1,2,y)*ln(1 - y)* a3^2*f1 - 16*S(1,2,y)*ln(1 - y)*a2*f1^2 - 16*S(1,2,y)*ln(1 - y )*a2*e1*f1 + 16*S(1,2,y)*ln(1 - y)*a2*a3*f1 + 8*S(1,2,y)*ln(1 - y)*a2^2*f1 - 16*S(1,2,y)*ln(1 - y)*a1*f1^2 - 16*S(1,2,y)* ln(1 - y)*a1*e1*f1 + 16*S(1,2,y)*ln(1 - y)*a1*a3*f1 + 16*S(1,2 ,y)*ln(1 - y)*a1*a2*f1 + 8*S(1,2,y)*ln(1 - y)*a1^2*f1 + 16*S(1 ,3, - y)*e1*f1^2 + 24*S(1,3, - y)*e1^2*f1 + 8*S(1,3, - y)*e1^3 - 16*S(1,3, - y)*a3*f1^2 - 48*S(1,3, - y)*a3*e1*f1 - 24*S(1,3 , - y)*a3*e1^2 + 24*S(1,3, - y)*a3^2*f1 + 24*S(1,3, - y)*a3^2* e1 - 8*S(1,3, - y)*a3^3 - 16*S(1,3, - y)*a2*f1^2 - 48*S(1,3, - y)*a2*e1*f1 - 24*S(1,3, - y)*a2*e1^2 + 48*S(1,3, - y)*a2*a3 *f1 + 48*S(1,3, - y)*a2*a3*e1 - 24*S(1,3, - y)*a2*a3^2 + 24*S( 1,3, - y)*a2^2*f1 + 24*S(1,3, - y)*a2^2*e1 - 24*S(1,3, - y)* a2^2*a3 - 8*S(1,3, - y)*a2^3 - 16*S(1,3, - y)*a1*f1^2 - 48*S(1 ,3, - y)*a1*e1*f1 - 24*S(1,3, - y)*a1*e1^2 + 48*S(1,3, - y)*a1 *a3*f1 + 48*S(1,3, - y)*a1*a3*e1 - 24*S(1,3, - y)*a1*a3^2 + 48 *S(1,3, - y)*a1*a2*f1 + 48*S(1,3, - y)*a1*a2*e1 - 48*S(1,3, - y)*a1*a2*a3 - 24*S(1,3, - y)*a1*a2^2 + 24*S(1,3, - y)*a1^2*f1 + 24*S(1,3, - y)*a1^2*e1 - 24*S(1,3, - y)*a1^2*a3 - 24*S(1,3, - y)*a1^2*a2 - 8*S(1,3, - y)*a1^3 - 2*S(1,3,y^2)*e1*f1^2 - 2* S(1,3,y^2)*e1^2*f1 + 2*S(1,3,y^2)*a3*f1^2 + 4*S(1,3,y^2)*a3*e1 *f1 - 2*S(1,3,y^2)*a3^2*f1 + 2*S(1,3,y^2)*a2*f1^2 + 4*S(1,3, y^2)*a2*e1*f1 - 4*S(1,3,y^2)*a2*a3*f1 - 2*S(1,3,y^2)*a2^2*f1 + 2*S(1,3,y^2)*a1*f1^2 + 4*S(1,3,y^2)*a1*e1*f1 - 4*S(1,3,y^2) *a1*a3*f1 - 4*S(1,3,y^2)*a1*a2*f1 - 2*S(1,3,y^2)*a1^2*f1 + 16* S(1,3,y)*e1*f1^2 + 8*S(1,3,y)*e1^2*f1 - 16*S(1,3,y)*a3*f1^2 - 16*S(1,3,y)*a3*e1*f1 + 8*S(1,3,y)*a3^2*f1 - 16*S(1,3,y)*a2* f1^2 - 16*S(1,3,y)*a2*e1*f1 + 16*S(1,3,y)*a2*a3*f1 + 8*S(1,3,y )*a2^2*f1 - 16*S(1,3,y)*a1*f1^2 - 16*S(1,3,y)*a1*e1*f1 + 16*S( 1,3,y)*a1*a3*f1 + 16*S(1,3,y)*a1*a2*f1 + 8*S(1,3,y)*a1^2*f1 - 8*S(2,2, - y)*e1*f1^2 - 12*S(2,2, - y)*e1^2*f1 - 4*S(2,2, - y) *e1^3 + 8*S(2,2, - y)*a3*f1^2 + 24*S(2,2, - y)*a3*e1*f1 + 12* S(2,2, - y)*a3*e1^2 - 12*S(2,2, - y)*a3^2*f1 - 12*S(2,2, - y)* a3^2*e1 + 4*S(2,2, - y)*a3^3 + 8*S(2,2, - y)*a2*f1^2 + 24*S(2, 2, - y)*a2*e1*f1 + 12*S(2,2, - y)*a2*e1^2 - 24*S(2,2, - y)*a2* a3*f1 - 24*S(2,2, - y)*a2*a3*e1 + 12*S(2,2, - y)*a2*a3^2 - 12* S(2,2, - y)*a2^2*f1 - 12*S(2,2, - y)*a2^2*e1 + 12*S(2,2, - y)* a2^2*a3 + 4*S(2,2, - y)*a2^3 + 8*S(2,2, - y)*a1*f1^2 + 24*S(2, 2, - y)*a1*e1*f1 + 12*S(2,2, - y)*a1*e1^2 - 24*S(2,2, - y)*a1* a3*f1 - 24*S(2,2, - y)*a1*a3*e1 + 12*S(2,2, - y)*a1*a3^2 - 24* S(2,2, - y)*a1*a2*f1 - 24*S(2,2, - y)*a1*a2*e1 + 24*S(2,2, - y )*a1*a2*a3 + 12*S(2,2, - y)*a1*a2^2 - 12*S(2,2, - y)*a1^2*f1 - 12*S(2,2, - y)*a1^2*e1 + 12*S(2,2, - y)*a1^2*a3 + 12*S(2,2, - y)*a1^2*a2 + 4*S(2,2, - y)*a1^3 + S(2,2,y^2)*e1*f1^2 + S(2, 2,y^2)*e1^2*f1 - S(2,2,y^2)*a3*f1^2 - 2*S(2,2,y^2)*a3*e1*f1 + S(2,2,y^2)*a3^2*f1 - S(2,2,y^2)*a2*f1^2 - 2*S(2,2,y^2)*a2*e1* f1 + 2*S(2,2,y^2)*a2*a3*f1 + S(2,2,y^2)*a2^2*f1 - S(2,2,y^2)* a1*f1^2 - 2*S(2,2,y^2)*a1*e1*f1 + 2*S(2,2,y^2)*a1*a3*f1 + 2*S( 2,2,y^2)*a1*a2*f1 + S(2,2,y^2)*a1^2*f1 - 8*S(2,2,y)*e1*f1^2 - 4*S(2,2,y)*e1^2*f1 + 8*S(2,2,y)*a3*f1^2 + 8*S(2,2,y)*a3*e1*f1 - 4*S(2,2,y)*a3^2*f1 + 8*S(2,2,y)*a2*f1^2 + 8*S(2,2,y)*a2*e1* f1 - 8*S(2,2,y)*a2*a3*f1 - 4*S(2,2,y)*a2^2*f1 + 8*S(2,2,y)*a1* f1^2 + 8*S(2,2,y)*a1*e1*f1 - 8*S(2,2,y)*a1*a3*f1 - 8*S(2,2,y)* a1*a2*f1 - 4*S(2,2,y)*a1^2*f1 - 7*ln(1 - y)^2*zt(2)*e1*f1^2 - 3*ln(1 - y)^2*zt(2)*e1^2*f1 + 7*ln(1 - y)^2*zt(2)*a3*f1^2 + 6* ln(1 - y)^2*zt(2)*a3*e1*f1 - 3*ln(1 - y)^2*zt(2)*a3^2*f1 + 7* ln(1 - y)^2*zt(2)*a2*f1^2 + 6*ln(1 - y)^2*zt(2)*a2*e1*f1 - 6* ln(1 - y)^2*zt(2)*a2*a3*f1 - 3*ln(1 - y)^2*zt(2)*a2^2*f1 + 7* ln(1 - y)^2*zt(2)*a1*f1^2 + 6*ln(1 - y)^2*zt(2)*a1*e1*f1 - 6* ln(1 - y)^2*zt(2)*a1*a3*f1 - 6*ln(1 - y)^2*zt(2)*a1*a2*f1 - 3* ln(1 - y)^2*zt(2)*a1^2*f1 - 2*ln(1 - y)^2*Li(2, - y)*e1*f1^2 - 2*ln(1 - y)^2*Li(2, - y)*e1^2*f1 + 2*ln(1 - y)^2*Li(2, - y) *a3*f1^2 + 4*ln(1 - y)^2*Li(2, - y)*a3*e1*f1 - 2*ln(1 - y)^2* Li(2, - y)*a3^2*f1 + 2*ln(1 - y)^2*Li(2, - y)*a2*f1^2 + 4*ln(1 - y)^2*Li(2, - y)*a2*e1*f1 - 4*ln(1 - y)^2*Li(2, - y)*a2*a3* f1 - 2*ln(1 - y)^2*Li(2, - y)*a2^2*f1 + 2*ln(1 - y)^2*Li(2, - y)*a1*f1^2 + 4*ln(1 - y)^2*Li(2, - y)*a1*e1*f1 - 4*ln(1 - y)^2 *Li(2, - y)*a1*a3*f1 - 4*ln(1 - y)^2*Li(2, - y)*a1*a2*f1 - 2* ln(1 - y)^2*Li(2, - y)*a1^2*f1 + 6*ln(1 - y)^2*Li(2,y)*e1*f1^2 + 2*ln(1 - y)^2*Li(2,y)*e1^2*f1 - 6*ln(1 - y)^2*Li(2,y)*a3* f1^2 - 4*ln(1 - y)^2*Li(2,y)*a3*e1*f1 + 2*ln(1 - y)^2*Li(2,y)* a3^2*f1 - 6*ln(1 - y)^2*Li(2,y)*a2*f1^2 - 4*ln(1 - y)^2*Li(2,y )*a2*e1*f1 + 4*ln(1 - y)^2*Li(2,y)*a2*a3*f1 + 2*ln(1 - y)^2* Li(2,y)*a2^2*f1 - 6*ln(1 - y)^2*Li(2,y)*a1*f1^2 - 4*ln(1 - y)^ 2*Li(2,y)*a1*e1*f1 + 4*ln(1 - y)^2*Li(2,y)*a1*a3*f1 + 4*ln(1 - y)^2*Li(2,y)*a1*a2*f1 + 2*ln(1 - y)^2*Li(2,y)*a1^2*f1 + 2* ln(1 - y)*ln(1 + y)*zt(2)*e1*f1^2 + 2*ln(1 - y)*ln(1 + y)*zt(2 )*e1^2*f1 - 2*ln(1 - y)*ln(1 + y)*zt(2)*a3*f1^2 - 4*ln(1 - y)* ln(1 + y)*zt(2)*a3*e1*f1 + 2*ln(1 - y)*ln(1 + y)*zt(2)*a3^2*f1 - 2*ln(1 - y)*ln(1 + y)*zt(2)*a2*f1^2 - 4*ln(1 - y)*ln(1 + y) *zt(2)*a2*e1*f1 + 4*ln(1 - y)*ln(1 + y)*zt(2)*a2*a3*f1 + 2*ln( 1 - y)*ln(1 + y)*zt(2)*a2^2*f1 - 2*ln(1 - y)*ln(1 + y)*zt(2)* a1*f1^2 - 4*ln(1 - y)*ln(1 + y)*zt(2)*a1*e1*f1 + 4*ln(1 - y)* ln(1 + y)*zt(2)*a1*a3*f1 + 4*ln(1 - y)*ln(1 + y)*zt(2)*a1*a2* f1 + 2*ln(1 - y)*ln(1 + y)*zt(2)*a1^2*f1 - 4*ln(1 - y)*ln(1 + y)*Li(2, - y)*e1*f1^2 - 4*ln(1 - y)*ln(1 + y)*Li(2, - y)*e1^2* f1 + 4*ln(1 - y)*ln(1 + y)*Li(2, - y)*a3*f1^2 + 8*ln(1 - y)* ln(1 + y)*Li(2, - y)*a3*e1*f1 - 4*ln(1 - y)*ln(1 + y)*Li(2, - y)*a3^2*f1 + 4*ln(1 - y)*ln(1 + y)*Li(2, - y)*a2*f1^2 + 8*ln(1 - y)*ln(1 + y)*Li(2, - y)*a2*e1*f1 - 8*ln(1 - y)*ln(1 + y)* Li(2, - y)*a2*a3*f1 - 4*ln(1 - y)*ln(1 + y)*Li(2, - y)*a2^2*f1 + 4*ln(1 - y)*ln(1 + y)*Li(2, - y)*a1*f1^2 + 8*ln(1 - y)*ln(1 + y)*Li(2, - y)*a1*e1*f1 - 8*ln(1 - y)*ln(1 + y)*Li(2, - y)* a1*a3*f1 - 8*ln(1 - y)*ln(1 + y)*Li(2, - y)*a1*a2*f1 - 4*ln(1 - y)*ln(1 + y)*Li(2, - y)*a1^2*f1 - 4*ln(1 - y)*ln(1 + y)*Li( 2,y)*e1*f1^2 - 4*ln(1 - y)*ln(1 + y)*Li(2,y)*e1^2*f1 + 4*ln(1 - y)*ln(1 + y)*Li(2,y)*a3*f1^2 + 8*ln(1 - y)*ln(1 + y)*Li(2,y )*a3*e1*f1 - 4*ln(1 - y)*ln(1 + y)*Li(2,y)*a3^2*f1 + 4*ln(1 - y)*ln(1 + y)*Li(2,y)*a2*f1^2 + 8*ln(1 - y)*ln(1 + y)*Li(2,y)* a2*e1*f1 - 8*ln(1 - y)*ln(1 + y)*Li(2,y)*a2*a3*f1 - 4*ln(1 - y )*ln(1 + y)*Li(2,y)*a2^2*f1 + 4*ln(1 - y)*ln(1 + y)*Li(2,y)*a1 *f1^2 + 8*ln(1 - y)*ln(1 + y)*Li(2,y)*a1*e1*f1 - 8*ln(1 - y)* ln(1 + y)*Li(2,y)*a1*a3*f1 - 8*ln(1 - y)*ln(1 + y)*Li(2,y)*a1* a2*f1 - 4*ln(1 - y)*ln(1 + y)*Li(2,y)*a1^2*f1 - 7*ln(1 - y)* zt(3)*e1*f1^2 - 3*ln(1 - y)*zt(3)*e1^2*f1 + 7*ln(1 - y)*zt(3)* a3*f1^2 + 6*ln(1 - y)*zt(3)*a3*e1*f1 - 3*ln(1 - y)*zt(3)*a3^2* f1 + 7*ln(1 - y)*zt(3)*a2*f1^2 + 6*ln(1 - y)*zt(3)*a2*e1*f1 - 6*ln(1 - y)*zt(3)*a2*a3*f1 - 3*ln(1 - y)*zt(3)*a2^2*f1 + 7*ln( 1 - y)*zt(3)*a1*f1^2 + 6*ln(1 - y)*zt(3)*a1*e1*f1 - 6*ln(1 - y )*zt(3)*a1*a3*f1 - 6*ln(1 - y)*zt(3)*a1*a2*f1 - 3*ln(1 - y)* zt(3)*a1^2*f1 + 4*ln(1 - y)*Li(3, - y)*e1*f1^2 + 4*ln(1 - y)* Li(3, - y)*e1^2*f1 - 4*ln(1 - y)*Li(3, - y)*a3*f1^2 - 8*ln(1 - y)*Li(3, - y)*a3*e1*f1 + 4*ln(1 - y)*Li(3, - y)*a3^2*f1 - 4 *ln(1 - y)*Li(3, - y)*a2*f1^2 - 8*ln(1 - y)*Li(3, - y)*a2*e1* f1 + 8*ln(1 - y)*Li(3, - y)*a2*a3*f1 + 4*ln(1 - y)*Li(3, - y)* a2^2*f1 - 4*ln(1 - y)*Li(3, - y)*a1*f1^2 - 8*ln(1 - y)*Li(3, - y)*a1*e1*f1 + 8*ln(1 - y)*Li(3, - y)*a1*a3*f1 + 8*ln(1 - y) *Li(3, - y)*a1*a2*f1 + 4*ln(1 - y)*Li(3, - y)*a1^2*f1 - 4*ln(1 - y)*Li(3,y)*e1*f1^2 + 4*ln(1 - y)*Li(3,y)*a3*f1^2 + 4*ln(1 - y)*Li(3,y)*a2*f1^2 + 4*ln(1 - y)*Li(3,y)*a1*f1^2 + 5*ln(1 + y)^2*zt(2)*e1*f1^2 + 7*ln(1 + y)^2*zt(2)*e1^2*f1 + 2*ln(1 + y)^2*zt(2)*e1^3 - 5*ln(1 + y)^2*zt(2)*a3*f1^2 - 14*ln(1 + y )^2*zt(2)*a3*e1*f1 - 6*ln(1 + y)^2*zt(2)*a3*e1^2 + 7*ln(1 + y) ^2*zt(2)*a3^2*f1 + 6*ln(1 + y)^2*zt(2)*a3^2*e1 - 2*ln(1 + y)^2 *zt(2)*a3^3 - 5*ln(1 + y)^2*zt(2)*a2*f1^2 - 14*ln(1 + y)^2*zt( 2)*a2*e1*f1 - 6*ln(1 + y)^2*zt(2)*a2*e1^2 + 14*ln(1 + y)^2*zt( 2)*a2*a3*f1 + 12*ln(1 + y)^2*zt(2)*a2*a3*e1 - 6*ln(1 + y)^2* zt(2)*a2*a3^2 + 7*ln(1 + y)^2*zt(2)*a2^2*f1 + 6*ln(1 + y)^2* zt(2)*a2^2*e1 - 6*ln(1 + y)^2*zt(2)*a2^2*a3 - 2*ln(1 + y)^2* zt(2)*a2^3 - 5*ln(1 + y)^2*zt(2)*a1*f1^2 - 14*ln(1 + y)^2*zt(2 )*a1*e1*f1 - 6*ln(1 + y)^2*zt(2)*a1*e1^2 + 14*ln(1 + y)^2*zt(2 )*a1*a3*f1 + 12*ln(1 + y)^2*zt(2)*a1*a3*e1 - 6*ln(1 + y)^2*zt( 2)*a1*a3^2 + 14*ln(1 + y)^2*zt(2)*a1*a2*f1 + 12*ln(1 + y)^2* zt(2)*a1*a2*e1 - 12*ln(1 + y)^2*zt(2)*a1*a2*a3 - 6*ln(1 + y)^2 *zt(2)*a1*a2^2 + 7*ln(1 + y)^2*zt(2)*a1^2*f1 + 6*ln(1 + y)^2* zt(2)*a1^2*e1 - 6*ln(1 + y)^2*zt(2)*a1^2*a3 - 6*ln(1 + y)^2* zt(2)*a1^2*a2 - 2*ln(1 + y)^2*zt(2)*a1^3 + 6*ln(1 + y)^2*Li(2, - y)*e1*f1^2 + 10*ln(1 + y)^2*Li(2, - y)*e1^2*f1 + 4*ln(1 + y )^2*Li(2, - y)*e1^3 - 6*ln(1 + y)^2*Li(2, - y)*a3*f1^2 - 20* ln(1 + y)^2*Li(2, - y)*a3*e1*f1 - 12*ln(1 + y)^2*Li(2, - y)*a3 *e1^2 + 10*ln(1 + y)^2*Li(2, - y)*a3^2*f1 + 12*ln(1 + y)^2*Li( 2, - y)*a3^2*e1 - 4*ln(1 + y)^2*Li(2, - y)*a3^3 - 6*ln(1 + y)^ 2*Li(2, - y)*a2*f1^2 - 20*ln(1 + y)^2*Li(2, - y)*a2*e1*f1 - 12 *ln(1 + y)^2*Li(2, - y)*a2*e1^2 + 20*ln(1 + y)^2*Li(2, - y)*a2 *a3*f1 + 24*ln(1 + y)^2*Li(2, - y)*a2*a3*e1 - 12*ln(1 + y)^2* Li(2, - y)*a2*a3^2 + 10*ln(1 + y)^2*Li(2, - y)*a2^2*f1 + 12* ln(1 + y)^2*Li(2, - y)*a2^2*e1 - 12*ln(1 + y)^2*Li(2, - y)* a2^2*a3 - 4*ln(1 + y)^2*Li(2, - y)*a2^3 - 6*ln(1 + y)^2*Li(2, - y)*a1*f1^2 - 20*ln(1 + y)^2*Li(2, - y)*a1*e1*f1 - 12*ln(1 + y)^2*Li(2, - y)*a1*e1^2 + 20*ln(1 + y)^2*Li(2, - y)*a1*a3* f1 + 24*ln(1 + y)^2*Li(2, - y)*a1*a3*e1 - 12*ln(1 + y)^2*Li(2, - y)*a1*a3^2 + 20*ln(1 + y)^2*Li(2, - y)*a1*a2*f1 + 24*ln(1 + y)^2*Li(2, - y)*a1*a2*e1 - 24*ln(1 + y)^2*Li(2, - y)*a1*a2* a3 - 12*ln(1 + y)^2*Li(2, - y)*a1*a2^2 + 10*ln(1 + y)^2*Li(2, - y)*a1^2*f1 + 12*ln(1 + y)^2*Li(2, - y)*a1^2*e1 - 12*ln(1 + y)^2*Li(2, - y)*a1^2*a3 - 12*ln(1 + y)^2*Li(2, - y)*a1^2*a2 - 4*ln(1 + y)^2*Li(2, - y)*a1^3 - 2*ln(1 + y)^2*Li(2,y)*e1*f1^2 - 2*ln(1 + y)^2*Li(2,y)*e1^2*f1 + 2*ln(1 + y)^2*Li(2,y)*a3* f1^2 + 4*ln(1 + y)^2*Li(2,y)*a3*e1*f1 - 2*ln(1 + y)^2*Li(2,y)* a3^2*f1 + 2*ln(1 + y)^2*Li(2,y)*a2*f1^2 + 4*ln(1 + y)^2*Li(2,y )*a2*e1*f1 - 4*ln(1 + y)^2*Li(2,y)*a2*a3*f1 - 2*ln(1 + y)^2* Li(2,y)*a2^2*f1 + 2*ln(1 + y)^2*Li(2,y)*a1*f1^2 + 4*ln(1 + y)^ 2*Li(2,y)*a1*e1*f1 - 4*ln(1 + y)^2*Li(2,y)*a1*a3*f1 - 4*ln(1 + y)^2*Li(2,y)*a1*a2*f1 - 2*ln(1 + y)^2*Li(2,y)*a1^2*f1 - 7* ln(1 + y)*zt(3)*e1*f1^2 - 11*ln(1 + y)*zt(3)*e1^2*f1 - 4*ln(1 + y)*zt(3)*e1^3 + 7*ln(1 + y)*zt(3)*a3*f1^2 + 22*ln(1 + y)* zt(3)*a3*e1*f1 + 12*ln(1 + y)*zt(3)*a3*e1^2 - 11*ln(1 + y)*zt( 3)*a3^2*f1 - 12*ln(1 + y)*zt(3)*a3^2*e1 + 4*ln(1 + y)*zt(3)* a3^3 + 7*ln(1 + y)*zt(3)*a2*f1^2 + 22*ln(1 + y)*zt(3)*a2*e1*f1 + 12*ln(1 + y)*zt(3)*a2*e1^2 - 22*ln(1 + y)*zt(3)*a2*a3*f1 - 24*ln(1 + y)*zt(3)*a2*a3*e1 + 12*ln(1 + y)*zt(3)*a2*a3^2 - 11* ln(1 + y)*zt(3)*a2^2*f1 - 12*ln(1 + y)*zt(3)*a2^2*e1 + 12*ln(1 + y)*zt(3)*a2^2*a3 + 4*ln(1 + y)*zt(3)*a2^3 + 7*ln(1 + y)*zt( 3)*a1*f1^2 + 22*ln(1 + y)*zt(3)*a1*e1*f1 + 12*ln(1 + y)*zt(3)* a1*e1^2 - 22*ln(1 + y)*zt(3)*a1*a3*f1 - 24*ln(1 + y)*zt(3)*a1* a3*e1 + 12*ln(1 + y)*zt(3)*a1*a3^2 - 22*ln(1 + y)*zt(3)*a1*a2* f1 - 24*ln(1 + y)*zt(3)*a1*a2*e1 + 24*ln(1 + y)*zt(3)*a1*a2*a3 + 12*ln(1 + y)*zt(3)*a1*a2^2 - 11*ln(1 + y)*zt(3)*a1^2*f1 - 12*ln(1 + y)*zt(3)*a1^2*e1 + 12*ln(1 + y)*zt(3)*a1^2*a3 + 12* ln(1 + y)*zt(3)*a1^2*a2 + 4*ln(1 + y)*zt(3)*a1^3 - 4*ln(1 + y) *Li(3, - y)*e1*f1^2 - 8*ln(1 + y)*Li(3, - y)*e1^2*f1 - 4*ln(1 + y)*Li(3, - y)*e1^3 + 4*ln(1 + y)*Li(3, - y)*a3*f1^2 + 16* ln(1 + y)*Li(3, - y)*a3*e1*f1 + 12*ln(1 + y)*Li(3, - y)*a3* e1^2 - 8*ln(1 + y)*Li(3, - y)*a3^2*f1 - 12*ln(1 + y)*Li(3, - y )*a3^2*e1 + 4*ln(1 + y)*Li(3, - y)*a3^3 + 4*ln(1 + y)*Li(3, - y)*a2*f1^2 + 16*ln(1 + y)*Li(3, - y)*a2*e1*f1 + 12*ln(1 + y)* Li(3, - y)*a2*e1^2 - 16*ln(1 + y)*Li(3, - y)*a2*a3*f1 - 24*ln( 1 + y)*Li(3, - y)*a2*a3*e1 + 12*ln(1 + y)*Li(3, - y)*a2*a3^2 - 8*ln(1 + y)*Li(3, - y)*a2^2*f1 - 12*ln(1 + y)*Li(3, - y)* a2^2*e1 + 12*ln(1 + y)*Li(3, - y)*a2^2*a3 + 4*ln(1 + y)*Li(3, - y)*a2^3 + 4*ln(1 + y)*Li(3, - y)*a1*f1^2 + 16*ln(1 + y)*Li( 3, - y)*a1*e1*f1 + 12*ln(1 + y)*Li(3, - y)*a1*e1^2 - 16*ln(1 + y)*Li(3, - y)*a1*a3*f1 - 24*ln(1 + y)*Li(3, - y)*a1*a3*e1 + 12*ln(1 + y)*Li(3, - y)*a1*a3^2 - 16*ln(1 + y)*Li(3, - y)* a1*a2*f1 - 24*ln(1 + y)*Li(3, - y)*a1*a2*e1 + 24*ln(1 + y)*Li( 3, - y)*a1*a2*a3 + 12*ln(1 + y)*Li(3, - y)*a1*a2^2 - 8*ln(1 + y)*Li(3, - y)*a1^2*f1 - 12*ln(1 + y)*Li(3, - y)*a1^2*e1 + 12* ln(1 + y)*Li(3, - y)*a1^2*a3 + 12*ln(1 + y)*Li(3, - y)*a1^2*a2 + 4*ln(1 + y)*Li(3, - y)*a1^3 + 4*ln(1 + y)*Li(3,y)*e1*f1^2 + 4*ln(1 + y)*Li(3,y)*e1^2*f1 - 4*ln(1 + y)*Li(3,y)*a3*f1^2 - 8*ln(1 + y)*Li(3,y)*a3*e1*f1 + 4*ln(1 + y)*Li(3,y)*a3^2*f1 - 4*ln(1 + y)*Li(3,y)*a2*f1^2 - 8*ln(1 + y)*Li(3,y)*a2*e1*f1 + 8*ln(1 + y)*Li(3,y)*a2*a3*f1 + 4*ln(1 + y)*Li(3,y)*a2^2*f1 - 4*ln(1 + y)*Li(3,y)*a1*f1^2 - 8*ln(1 + y)*Li(3,y)*a1*e1*f1 + 8*ln(1 + y)*Li(3,y)*a1*a3*f1 + 8*ln(1 + y)*Li(3,y)*a1*a2*f1 + 4*ln(1 + y)*Li(3,y)*a1^2*f1 + 2*ln(y)*ln(1 - y)^3*e1*f1^2 + 2/3*ln(y)*ln(1 - y)^3*e1^2*f1 - 2*ln(y)*ln(1 - y)^3*a3*f1^2 - 4/3*ln(y)*ln(1 - y)^3*a3*e1*f1 + 2/3*ln(y)*ln(1 - y)^3*a3^2 *f1 - 2*ln(y)*ln(1 - y)^3*a2*f1^2 - 4/3*ln(y)*ln(1 - y)^3*a2* e1*f1 + 4/3*ln(y)*ln(1 - y)^3*a2*a3*f1 + 2/3*ln(y)*ln(1 - y)^3 *a2^2*f1 - 2*ln(y)*ln(1 - y)^3*a1*f1^2 - 4/3*ln(y)*ln(1 - y)^3 *a1*e1*f1 + 4/3*ln(y)*ln(1 - y)^3*a1*a3*f1 + 4/3*ln(y)*ln(1 - y)^3*a1*a2*f1 + 2/3*ln(y)*ln(1 - y)^3*a1^2*f1 - 2*ln(y)*ln(1 - y)^2*ln(1 + y)*e1*f1^2 - 2*ln(y)*ln(1 - y)^2*ln(1 + y)*e1^2 *f1 + 2*ln(y)*ln(1 - y)^2*ln(1 + y)*a3*f1^2 + 4*ln(y)*ln(1 - y )^2*ln(1 + y)*a3*e1*f1 - 2*ln(y)*ln(1 - y)^2*ln(1 + y)*a3^2*f1 + 2*ln(y)*ln(1 - y)^2*ln(1 + y)*a2*f1^2 + 4*ln(y)*ln(1 - y)^2 *ln(1 + y)*a2*e1*f1 - 4*ln(y)*ln(1 - y)^2*ln(1 + y)*a2*a3*f1 - 2*ln(y)*ln(1 - y)^2*ln(1 + y)*a2^2*f1 + 2*ln(y)*ln(1 - y)^2 *ln(1 + y)*a1*f1^2 + 4*ln(y)*ln(1 - y)^2*ln(1 + y)*a1*e1*f1 - 4*ln(y)*ln(1 - y)^2*ln(1 + y)*a1*a3*f1 - 4*ln(y)*ln(1 - y)^2* ln(1 + y)*a1*a2*f1 - 2*ln(y)*ln(1 - y)^2*ln(1 + y)*a1^2*f1 - 2 *ln(y)*ln(1 - y)*ln(1 + y)^2*e1*f1^2 - 2*ln(y)*ln(1 - y)*ln(1 + y)^2*e1^2*f1 + 2*ln(y)*ln(1 - y)*ln(1 + y)^2*a3*f1^2 + 4* ln(y)*ln(1 - y)*ln(1 + y)^2*a3*e1*f1 - 2*ln(y)*ln(1 - y)*ln(1 + y)^2*a3^2*f1 + 2*ln(y)*ln(1 - y)*ln(1 + y)^2*a2*f1^2 + 4* ln(y)*ln(1 - y)*ln(1 + y)^2*a2*e1*f1 - 4*ln(y)*ln(1 - y)*ln(1 + y)^2*a2*a3*f1 - 2*ln(y)*ln(1 - y)*ln(1 + y)^2*a2^2*f1 + 2* ln(y)*ln(1 - y)*ln(1 + y)^2*a1*f1^2 + 4*ln(y)*ln(1 - y)*ln(1 + y)^2*a1*e1*f1 - 4*ln(y)*ln(1 - y)*ln(1 + y)^2*a1*a3*f1 - 4* ln(y)*ln(1 - y)*ln(1 + y)^2*a1*a2*f1 - 2*ln(y)*ln(1 - y)*ln(1 + y)^2*a1^2*f1 + 6*ln(y)*ln(1 - y)*zt(2)*e1*f1^2 + 2*ln(y)* ln(1 - y)*zt(2)*e1^2*f1 - 6*ln(y)*ln(1 - y)*zt(2)*a3*f1^2 - 4* ln(y)*ln(1 - y)*zt(2)*a3*e1*f1 + 2*ln(y)*ln(1 - y)*zt(2)*a3^2* f1 - 6*ln(y)*ln(1 - y)*zt(2)*a2*f1^2 - 4*ln(y)*ln(1 - y)*zt(2) *a2*e1*f1 + 4*ln(y)*ln(1 - y)*zt(2)*a2*a3*f1 + 2*ln(y)*ln(1 - y)*zt(2)*a2^2*f1 - 6*ln(y)*ln(1 - y)*zt(2)*a1*f1^2 - 4*ln(y)* ln(1 - y)*zt(2)*a1*e1*f1 + 4*ln(y)*ln(1 - y)*zt(2)*a1*a3*f1 + 4*ln(y)*ln(1 - y)*zt(2)*a1*a2*f1 + 2*ln(y)*ln(1 - y)*zt(2)* a1^2*f1 + 2*ln(y)*ln(1 + y)^3*e1*f1^2 + 10/3*ln(y)*ln(1 + y)^3 *e1^2*f1 + 4/3*ln(y)*ln(1 + y)^3*e1^3 - 2*ln(y)*ln(1 + y)^3*a3 *f1^2 - 20/3*ln(y)*ln(1 + y)^3*a3*e1*f1 - 4*ln(y)*ln(1 + y)^3* a3*e1^2 + 10/3*ln(y)*ln(1 + y)^3*a3^2*f1 + 4*ln(y)*ln(1 + y)^3 *a3^2*e1 - 4/3*ln(y)*ln(1 + y)^3*a3^3 - 2*ln(y)*ln(1 + y)^3*a2 *f1^2 - 20/3*ln(y)*ln(1 + y)^3*a2*e1*f1 - 4*ln(y)*ln(1 + y)^3* a2*e1^2 + 20/3*ln(y)*ln(1 + y)^3*a2*a3*f1 + 8*ln(y)*ln(1 + y)^ 3*a2*a3*e1 - 4*ln(y)*ln(1 + y)^3*a2*a3^2 + 10/3*ln(y)*ln(1 + y )^3*a2^2*f1 + 4*ln(y)*ln(1 + y)^3*a2^2*e1 - 4*ln(y)*ln(1 + y)^ 3*a2^2*a3 - 4/3*ln(y)*ln(1 + y)^3*a2^3 - 2*ln(y)*ln(1 + y)^3* a1*f1^2 - 20/3*ln(y)*ln(1 + y)^3*a1*e1*f1 - 4*ln(y)*ln(1 + y)^ 3*a1*e1^2 + 20/3*ln(y)*ln(1 + y)^3*a1*a3*f1 + 8*ln(y)*ln(1 + y )^3*a1*a3*e1 - 4*ln(y)*ln(1 + y)^3*a1*a3^2 + 20/3*ln(y)*ln(1 + y)^3*a1*a2*f1 + 8*ln(y)*ln(1 + y)^3*a1*a2*e1 - 8*ln(y)*ln(1 + y)^3*a1*a2*a3 - 4*ln(y)*ln(1 + y)^3*a1*a2^2 + 10/3*ln(y)* ln(1 + y)^3*a1^2*f1 + 4*ln(y)*ln(1 + y)^3*a1^2*e1 - 4*ln(y)* ln(1 + y)^3*a1^2*a3 - 4*ln(y)*ln(1 + y)^3*a1^2*a2 - 4/3*ln(y)* ln(1 + y)^3*a1^3 - 6*ln(y)*ln(1 + y)*zt(2)*e1*f1^2 - 8*ln(y)* ln(1 + y)*zt(2)*e1^2*f1 - 2*ln(y)*ln(1 + y)*zt(2)*e1^3 + 6*ln( y)*ln(1 + y)*zt(2)*a3*f1^2 + 16*ln(y)*ln(1 + y)*zt(2)*a3*e1*f1 + 6*ln(y)*ln(1 + y)*zt(2)*a3*e1^2 - 8*ln(y)*ln(1 + y)*zt(2)* a3^2*f1 - 6*ln(y)*ln(1 + y)*zt(2)*a3^2*e1 + 2*ln(y)*ln(1 + y)* zt(2)*a3^3 + 6*ln(y)*ln(1 + y)*zt(2)*a2*f1^2 + 16*ln(y)*ln(1 + y)*zt(2)*a2*e1*f1 + 6*ln(y)*ln(1 + y)*zt(2)*a2*e1^2 - 16* ln(y)*ln(1 + y)*zt(2)*a2*a3*f1 - 12*ln(y)*ln(1 + y)*zt(2)*a2* a3*e1 + 6*ln(y)*ln(1 + y)*zt(2)*a2*a3^2 - 8*ln(y)*ln(1 + y)* zt(2)*a2^2*f1 - 6*ln(y)*ln(1 + y)*zt(2)*a2^2*e1 + 6*ln(y)*ln(1 + y)*zt(2)*a2^2*a3 + 2*ln(y)*ln(1 + y)*zt(2)*a2^3 + 6*ln(y)* ln(1 + y)*zt(2)*a1*f1^2 + 16*ln(y)*ln(1 + y)*zt(2)*a1*e1*f1 + 6*ln(y)*ln(1 + y)*zt(2)*a1*e1^2 - 16*ln(y)*ln(1 + y)*zt(2)*a1* a3*f1 - 12*ln(y)*ln(1 + y)*zt(2)*a1*a3*e1 + 6*ln(y)*ln(1 + y)* zt(2)*a1*a3^2 - 16*ln(y)*ln(1 + y)*zt(2)*a1*a2*f1 - 12*ln(y)* ln(1 + y)*zt(2)*a1*a2*e1 + 12*ln(y)*ln(1 + y)*zt(2)*a1*a2*a3 + 6*ln(y)*ln(1 + y)*zt(2)*a1*a2^2 - 8*ln(y)*ln(1 + y)*zt(2)* a1^2*f1 - 6*ln(y)*ln(1 + y)*zt(2)*a1^2*e1 + 6*ln(y)*ln(1 + y)* zt(2)*a1^2*a3 + 6*ln(y)*ln(1 + y)*zt(2)*a1^2*a2 + 2*ln(y)*ln(1 + y)*zt(2)*a1^3 - ln(y)^2*ln(1 - y)^2*e1*f1^2 + ln(y)^2*ln(1 - y)^2*a3*f1^2 + ln(y)^2*ln(1 - y)^2*a2*f1^2 + ln(y)^2*ln(1 - y)^2*a1*f1^2 + 2*ln(y)^2*ln(1 - y)*ln(1 + y)*e1*f1^2 + 2* ln(y)^2*ln(1 - y)*ln(1 + y)*e1^2*f1 - 2*ln(y)^2*ln(1 - y)*ln(1 + y)*a3*f1^2 - 4*ln(y)^2*ln(1 - y)*ln(1 + y)*a3*e1*f1 + 2*ln( y)^2*ln(1 - y)*ln(1 + y)*a3^2*f1 - 2*ln(y)^2*ln(1 - y)*ln(1 + y)*a2*f1^2 - 4*ln(y)^2*ln(1 - y)*ln(1 + y)*a2*e1*f1 + 4*ln(y)^ 2*ln(1 - y)*ln(1 + y)*a2*a3*f1 + 2*ln(y)^2*ln(1 - y)*ln(1 + y) *a2^2*f1 - 2*ln(y)^2*ln(1 - y)*ln(1 + y)*a1*f1^2 - 4*ln(y)^2* ln(1 - y)*ln(1 + y)*a1*e1*f1 + 4*ln(y)^2*ln(1 - y)*ln(1 + y)* a1*a3*f1 + 4*ln(y)^2*ln(1 - y)*ln(1 + y)*a1*a2*f1 + 2*ln(y)^2* ln(1 - y)*ln(1 + y)*a1^2*f1 - ln(y)^2*ln(1 + y)^2*e1*f1^2 - 2* ln(y)^2*ln(1 + y)^2*e1^2*f1 - ln(y)^2*ln(1 + y)^2*e1^3 + ln(y) ^2*ln(1 + y)^2*a3*f1^2 + 4*ln(y)^2*ln(1 + y)^2*a3*e1*f1 + 3* ln(y)^2*ln(1 + y)^2*a3*e1^2 - 2*ln(y)^2*ln(1 + y)^2*a3^2*f1 - 3*ln(y)^2*ln(1 + y)^2*a3^2*e1 + ln(y)^2*ln(1 + y)^2*a3^3 + ln( y)^2*ln(1 + y)^2*a2*f1^2 + 4*ln(y)^2*ln(1 + y)^2*a2*e1*f1 + 3* ln(y)^2*ln(1 + y)^2*a2*e1^2 - 4*ln(y)^2*ln(1 + y)^2*a2*a3*f1 - 6*ln(y)^2*ln(1 + y)^2*a2*a3*e1 + 3*ln(y)^2*ln(1 + y)^2*a2* a3^2 - 2*ln(y)^2*ln(1 + y)^2*a2^2*f1 - 3*ln(y)^2*ln(1 + y)^2* a2^2*e1 + 3*ln(y)^2*ln(1 + y)^2*a2^2*a3 + ln(y)^2*ln(1 + y)^2* a2^3 + ln(y)^2*ln(1 + y)^2*a1*f1^2 + 4*ln(y)^2*ln(1 + y)^2*a1* e1*f1 + 3*ln(y)^2*ln(1 + y)^2*a1*e1^2 - 4*ln(y)^2*ln(1 + y)^2* a1*a3*f1 - 6*ln(y)^2*ln(1 + y)^2*a1*a3*e1 + 3*ln(y)^2*ln(1 + y )^2*a1*a3^2 - 4*ln(y)^2*ln(1 + y)^2*a1*a2*f1 - 6*ln(y)^2*ln(1 + y)^2*a1*a2*e1 + 6*ln(y)^2*ln(1 + y)^2*a1*a2*a3 + 3*ln(y)^2* ln(1 + y)^2*a1*a2^2 - 2*ln(y)^2*ln(1 + y)^2*a1^2*f1 - 3*ln(y)^ 2*ln(1 + y)^2*a1^2*e1 + 3*ln(y)^2*ln(1 + y)^2*a1^2*a3 + 3*ln(y )^2*ln(1 + y)^2*a1^2*a2 + ln(y)^2*ln(1 + y)^2*a1^3 + 1/3*ln(y) ^3*ln(1 - y)*a3*e1*f1 - 1/3*ln(y)^3*ln(1 - y)*a3^2*f1 + 1/3* ln(y)^3*ln(1 - y)*a2*e1*f1 - 1/3*ln(y)^3*ln(1 - y)*a2*a3*f1 - 1/3*ln(y)^3*ln(1 - y)*a2^2*f1 + 1/3*ln(y)^3*ln(1 - y)*a1*e1*f1 - 1/3*ln(y)^3*ln(1 - y)*a1*a3*f1 - 1/3*ln(y)^3*ln(1 - y)*a1* a2*f1 - 1/3*ln(y)^3*ln(1 - y)*a1^2*f1 - 1/3*ln(y)^3*ln(1 + y)* a3*e1*f1 - 1/3*ln(y)^3*ln(1 + y)*a3*e1^2 + 1/3*ln(y)^3*ln(1 + y)*a3^2*f1 + 2/3*ln(y)^3*ln(1 + y)*a3^2*e1 - 1/3*ln(y)^3*ln(1 + y)*a3^3 - 1/3*ln(y)^3*ln(1 + y)*a2*e1*f1 - 1/3*ln(y)^3*ln(1 + y)*a2*e1^2 + 1/3*ln(y)^3*ln(1 + y)*a2*a3*f1 + ln(y)^3*ln(1 + y)*a2*a3*e1 - 2/3*ln(y)^3*ln(1 + y)*a2*a3^2 + 1/3*ln(y)^3* ln(1 + y)*a2^2*f1 + 2/3*ln(y)^3*ln(1 + y)*a2^2*e1 - 2/3*ln(y)^ 3*ln(1 + y)*a2^2*a3 - 1/3*ln(y)^3*ln(1 + y)*a2^3 - 1/3*ln(y)^3 *ln(1 + y)*a1*e1*f1 - 1/3*ln(y)^3*ln(1 + y)*a1*e1^2 + 1/3*ln(y )^3*ln(1 + y)*a1*a3*f1 + ln(y)^3*ln(1 + y)*a1*a3*e1 - 2/3*ln(y )^3*ln(1 + y)*a1*a3^2 + 1/3*ln(y)^3*ln(1 + y)*a1*a2*f1 + ln(y) ^3*ln(1 + y)*a1*a2*e1 - ln(y)^3*ln(1 + y)*a1*a2*a3 - 2/3*ln(y) ^3*ln(1 + y)*a1*a2^2 + 1/3*ln(y)^3*ln(1 + y)*a1^2*f1 + 2/3*ln( y)^3*ln(1 + y)*a1^2*e1 - 2/3*ln(y)^3*ln(1 + y)*a1^2*a3 - 2/3* ln(y)^3*ln(1 + y)*a1^2*a2 - 1/3*ln(y)^3*ln(1 + y)*a1^3 - 1/24* ln(y)^4*a3^2*e1 + 1/24*ln(y)^4*a3^3 - 1/24*ln(y)^4*a2*a3*e1 + 1/24*ln(y)^4*a2*a3^2 - 1/24*ln(y)^4*a2^2*e1 + 1/24*ln(y)^4* a2^2*a3 + 1/24*ln(y)^4*a2^3 - 1/24*ln(y)^4*a1*a3*e1 + 1/24*ln( y)^4*a1*a3^2 - 1/24*ln(y)^4*a1*a2*e1 + 1/24*ln(y)^4*a1*a2*a3 + 1/24*ln(y)^4*a1*a2^2 - 1/24*ln(y)^4*a1^2*e1 + 1/24*ln(y)^4* a1^2*a3 + 1/24*ln(y)^4*a1^2*a2 + 1/24*ln(y)^4*a1^3 - 3/2*ln(y) ^2*zt(2)*a3*e1*f1 - 1/2*ln(y)^2*zt(2)*a3*e1^2 + 3/2*ln(y)^2* zt(2)*a3^2*f1 + ln(y)^2*zt(2)*a3^2*e1 - 1/2*ln(y)^2*zt(2)*a3^3 - 3/2*ln(y)^2*zt(2)*a2*e1*f1 - 1/2*ln(y)^2*zt(2)*a2*e1^2 + 3/ 2*ln(y)^2*zt(2)*a2*a3*f1 + 3/2*ln(y)^2*zt(2)*a2*a3*e1 - ln(y)^ 2*zt(2)*a2*a3^2 + 3/2*ln(y)^2*zt(2)*a2^2*f1 + ln(y)^2*zt(2)* a2^2*e1 - ln(y)^2*zt(2)*a2^2*a3 - 1/2*ln(y)^2*zt(2)*a2^3 - 3/2 *ln(y)^2*zt(2)*a1*e1*f1 - 1/2*ln(y)^2*zt(2)*a1*e1^2 + 3/2*ln(y )^2*zt(2)*a1*a3*f1 + 3/2*ln(y)^2*zt(2)*a1*a3*e1 - ln(y)^2*zt(2 )*a1*a3^2 + 3/2*ln(y)^2*zt(2)*a1*a2*f1 + 3/2*ln(y)^2*zt(2)*a1* a2*e1 - 3/2*ln(y)^2*zt(2)*a1*a2*a3 - ln(y)^2*zt(2)*a1*a2^2 + 3/ 2*ln(y)^2*zt(2)*a1^2*f1 + ln(y)^2*zt(2)*a1^2*e1 - ln(y)^2*zt(2 )*a1^2*a3 - ln(y)^2*zt(2)*a1^2*a2 - 1/2*ln(y)^2*zt(2)*a1^3 - ln(y)^2*Li(2, - y)*e1^2*f1 - ln(y)^2*Li(2, - y)*e1^3 + ln(y)^2 *Li(2, - y)*a3*e1*f1 + 2*ln(y)^2*Li(2, - y)*a3*e1^2 - ln(y)^2* Li(2, - y)*a3^2*e1 + ln(y)^2*Li(2, - y)*a2*e1*f1 + 2*ln(y)^2* Li(2, - y)*a2*e1^2 - ln(y)^2*Li(2, - y)*a2*a3*f1 - 3*ln(y)^2* Li(2, - y)*a2*a3*e1 + ln(y)^2*Li(2, - y)*a2*a3^2 - ln(y)^2*Li( 2, - y)*a2^2*e1 + ln(y)^2*Li(2, - y)*a2^2*a3 + ln(y)^2*Li(2, - y)*a1*e1*f1 + 2*ln(y)^2*Li(2, - y)*a1*e1^2 - ln(y)^2*Li(2, - y)*a1*a3*f1 - 3*ln(y)^2*Li(2, - y)*a1*a3*e1 + ln(y)^2*Li(2, - y)*a1*a3^2 - ln(y)^2*Li(2, - y)*a1*a2*f1 - 3*ln(y)^2*Li(2, - y)*a1*a2*e1 + 3*ln(y)^2*Li(2, - y)*a1*a2*a3 + ln(y)^2*Li(2, - y)*a1*a2^2 - ln(y)^2*Li(2, - y)*a1^2*e1 + ln(y)^2*Li(2, - y )*a1^2*a3 + ln(y)^2*Li(2, - y)*a1^2*a2 + ln(y)^2*Li(2,y)*e1^2* f1 - ln(y)^2*Li(2,y)*e1^3 - ln(y)^2*Li(2,y)*a3*e1*f1 + ln(y)^2 *Li(2,y)*a3*e1^2 - ln(y)^2*Li(2,y)*a2*e1*f1 + ln(y)^2*Li(2,y)* a2*e1^2 + ln(y)^2*Li(2,y)*a2*a3*f1 - ln(y)^2*Li(2,y)*a2*a3*e1 - ln(y)^2*Li(2,y)*a1*e1*f1 + ln(y)^2*Li(2,y)*a1*e1^2 + ln(y)^ 2*Li(2,y)*a1*a3*f1 - ln(y)^2*Li(2,y)*a1*a3*e1 + ln(y)^2*Li(2,y )*a1*a2*f1 - ln(y)^2*Li(2,y)*a1*a2*e1 + ln(y)^2*Li(2,y)*a1*a2* a3 + 7*ln(y)*zt(3)*e1*f1^2 + ln(y)*zt(3)*e1^3 - 7*ln(y)*zt(3)* a3*f1^2 - 7*ln(y)*zt(3)*a3*e1*f1 - 2*ln(y)*zt(3)*a3*e1^2 + 7* ln(y)*zt(3)*a3^2*f1 + 3*ln(y)*zt(3)*a3^2*e1 - 2*ln(y)*zt(3)* a3^3 - 7*ln(y)*zt(3)*a2*f1^2 - 7*ln(y)*zt(3)*a2*e1*f1 - 2*ln(y )*zt(3)*a2*e1^2 + 7*ln(y)*zt(3)*a2*a3*f1 + 5*ln(y)*zt(3)*a2*a3 *e1 - 3*ln(y)*zt(3)*a2*a3^2 + 7*ln(y)*zt(3)*a2^2*f1 + 3*ln(y)* zt(3)*a2^2*e1 - 3*ln(y)*zt(3)*a2^2*a3 - 2*ln(y)*zt(3)*a2^3 - 7 *ln(y)*zt(3)*a1*f1^2 - 7*ln(y)*zt(3)*a1*e1*f1 - 2*ln(y)*zt(3)* a1*e1^2 + 7*ln(y)*zt(3)*a1*a3*f1 + 5*ln(y)*zt(3)*a1*a3*e1 - 3* ln(y)*zt(3)*a1*a3^2 + 7*ln(y)*zt(3)*a1*a2*f1 + 5*ln(y)*zt(3)* a1*a2*e1 - 5*ln(y)*zt(3)*a1*a2*a3 - 3*ln(y)*zt(3)*a1*a2^2 + 7* ln(y)*zt(3)*a1^2*f1 + 3*ln(y)*zt(3)*a1^2*e1 - 3*ln(y)*zt(3)* a1^2*a3 - 3*ln(y)*zt(3)*a1^2*a2 - 2*ln(y)*zt(3)*a1^3 + 4*ln(y) *Li(3, - y)*e1^2*f1 + 4*ln(y)*Li(3, - y)*e1^3 - 4*ln(y)*Li(3, - y)*a3*e1*f1 - 8*ln(y)*Li(3, - y)*a3*e1^2 + 4*ln(y)*Li(3, - y)*a3^2*e1 - 4*ln(y)*Li(3, - y)*a2*e1*f1 - 8*ln(y)*Li(3, - y)* a2*e1^2 + 4*ln(y)*Li(3, - y)*a2*a3*f1 + 12*ln(y)*Li(3, - y)*a2 *a3*e1 - 4*ln(y)*Li(3, - y)*a2*a3^2 + 4*ln(y)*Li(3, - y)*a2^2* e1 - 4*ln(y)*Li(3, - y)*a2^2*a3 - 4*ln(y)*Li(3, - y)*a1*e1*f1 - 8*ln(y)*Li(3, - y)*a1*e1^2 + 4*ln(y)*Li(3, - y)*a1*a3*f1 + 12*ln(y)*Li(3, - y)*a1*a3*e1 - 4*ln(y)*Li(3, - y)*a1*a3^2 + 4* ln(y)*Li(3, - y)*a1*a2*f1 + 12*ln(y)*Li(3, - y)*a1*a2*e1 - 12* ln(y)*Li(3, - y)*a1*a2*a3 - 4*ln(y)*Li(3, - y)*a1*a2^2 + 4*ln( y)*Li(3, - y)*a1^2*e1 - 4*ln(y)*Li(3, - y)*a1^2*a3 - 4*ln(y)* Li(3, - y)*a1^2*a2 - 4*ln(y)*Li(3,y)*e1^2*f1 + 4*ln(y)*Li(3,y) *e1^3 + 4*ln(y)*Li(3,y)*a3*e1*f1 - 4*ln(y)*Li(3,y)*a3*e1^2 + 4 *ln(y)*Li(3,y)*a2*e1*f1 - 4*ln(y)*Li(3,y)*a2*e1^2 - 4*ln(y)* Li(3,y)*a2*a3*f1 + 4*ln(y)*Li(3,y)*a2*a3*e1 + 4*ln(y)*Li(3,y)* a1*e1*f1 - 4*ln(y)*Li(3,y)*a1*e1^2 - 4*ln(y)*Li(3,y)*a1*a3*f1 + 4*ln(y)*Li(3,y)*a1*a3*e1 - 4*ln(y)*Li(3,y)*a1*a2*f1 + 4*ln( y)*Li(3,y)*a1*a2*e1 - 4*ln(y)*Li(3,y)*a1*a2*a3 - 45/4*zt(4)*e1 *f1^2 - 15*zt(4)*e1^2*f1 + 5/4*zt(4)*e1^3 + 45/4*zt(4)*a3*f1^2 + 15*zt(4)*a3*e1*f1 + 5/4*zt(4)*a3*e1^2 - 1/4*zt(4)*a3^2*e1 - 9/4*zt(4)*a3^3 + 45/4*zt(4)*a2*f1^2 + 15*zt(4)*a2*e1*f1 + 5/ 4*zt(4)*a2*e1^2 - 15*zt(4)*a2*a3*f1 - 3/2*zt(4)*a2*a3*e1 + 1/4 *zt(4)*a2*a3^2 - 1/4*zt(4)*a2^2*e1 + 1/4*zt(4)*a2^2*a3 - 9/4* zt(4)*a2^3 + 45/4*zt(4)*a1*f1^2 + 15*zt(4)*a1*e1*f1 + 5/4*zt(4 )*a1*e1^2 - 15*zt(4)*a1*a3*f1 - 3/2*zt(4)*a1*a3*e1 + 1/4*zt(4) *a1*a3^2 - 15*zt(4)*a1*a2*f1 - 3/2*zt(4)*a1*a2*e1 + 3/2*zt(4)* a1*a2*a3 + 1/4*zt(4)*a1*a2^2 - 1/4*zt(4)*a1^2*e1 + 1/4*zt(4)* a1^2*a3 + 1/4*zt(4)*a1^2*a2 - 9/4*zt(4)*a1^3 - 8*Li(4, - [1-y] *[1+y]^-1)*f1^3 - 12*Li(4, - [1-y]*[1+y]^-1)*e1*f1^2 - 4*Li(4, - [1-y]*[1+y]^-1)*e1^2*f1 + 12*Li(4, - [1-y]*[1+y]^-1)*a3* f1^2 + 8*Li(4, - [1-y]*[1+y]^-1)*a3*e1*f1 - 4*Li(4, - [1-y]* [1+y]^-1)*a3^2*f1 + 12*Li(4, - [1-y]*[1+y]^-1)*a2*f1^2 + 8*Li( 4, - [1-y]*[1+y]^-1)*a2*e1*f1 - 8*Li(4, - [1-y]*[1+y]^-1)*a2* a3*f1 - 4*Li(4, - [1-y]*[1+y]^-1)*a2^2*f1 + 12*Li(4, - [1-y]* [1+y]^-1)*a1*f1^2 + 8*Li(4, - [1-y]*[1+y]^-1)*a1*e1*f1 - 8*Li( 4, - [1-y]*[1+y]^-1)*a1*a3*f1 - 8*Li(4, - [1-y]*[1+y]^-1)*a1* a2*f1 - 4*Li(4, - [1-y]*[1+y]^-1)*a1^2*f1 + 8*Li(4,[1-y]* [1+y]^-1)*f1^3 + 12*Li(4,[1-y]*[1+y]^-1)*e1*f1^2 + 4*Li(4, [1-y]*[1+y]^-1)*e1^2*f1 - 12*Li(4,[1-y]*[1+y]^-1)*a3*f1^2 - 8* Li(4,[1-y]*[1+y]^-1)*a3*e1*f1 + 4*Li(4,[1-y]*[1+y]^-1)*a3^2*f1 - 12*Li(4,[1-y]*[1+y]^-1)*a2*f1^2 - 8*Li(4,[1-y]*[1+y]^-1)*a2 *e1*f1 + 8*Li(4,[1-y]*[1+y]^-1)*a2*a3*f1 + 4*Li(4,[1-y]* [1+y]^-1)*a2^2*f1 - 12*Li(4,[1-y]*[1+y]^-1)*a1*f1^2 - 8*Li(4, [1-y]*[1+y]^-1)*a1*e1*f1 + 8*Li(4,[1-y]*[1+y]^-1)*a1*a3*f1 + 8 *Li(4,[1-y]*[1+y]^-1)*a1*a2*f1 + 4*Li(4,[1-y]*[1+y]^-1)*a1^2* f1 - 6*Li(4, - y)*e1^2*f1 - 6*Li(4, - y)*e1^3 + 4*Li(4, - y)* a3*e1*f1 + 10*Li(4, - y)*a3*e1^2 + 2*Li(4, - y)*a3^2*f1 - 2* Li(4, - y)*a3^2*e1 - 2*Li(4, - y)*a3^3 + 4*Li(4, - y)*a2*e1*f1 + 10*Li(4, - y)*a2*e1^2 - 4*Li(4, - y)*a2*a3*f1 - 12*Li(4, - y)*a2*a3*e1 + 2*Li(4, - y)*a2*a3^2 + 2*Li(4, - y)*a2^2*f1 - 2* Li(4, - y)*a2^2*e1 + 2*Li(4, - y)*a2^2*a3 - 2*Li(4, - y)*a2^3 + 4*Li(4, - y)*a1*e1*f1 + 10*Li(4, - y)*a1*e1^2 - 4*Li(4, - y )*a1*a3*f1 - 12*Li(4, - y)*a1*a3*e1 + 2*Li(4, - y)*a1*a3^2 - 4 *Li(4, - y)*a1*a2*f1 - 12*Li(4, - y)*a1*a2*e1 + 12*Li(4, - y)* a1*a2*a3 + 2*Li(4, - y)*a1*a2^2 + 2*Li(4, - y)*a1^2*f1 - 2*Li( 4, - y)*a1^2*e1 + 2*Li(4, - y)*a1^2*a3 + 2*Li(4, - y)*a1^2*a2 - 2*Li(4, - y)*a1^3 + 6*Li(4,y)*e1^2*f1 - 6*Li(4,y)*e1^3 - 4* Li(4,y)*a3*e1*f1 + 6*Li(4,y)*a3*e1^2 - 2*Li(4,y)*a3^2*f1 - 4* Li(4,y)*a2*e1*f1 + 6*Li(4,y)*a2*e1^2 + 4*Li(4,y)*a2*a3*f1 - 6* Li(4,y)*a2*a3*e1 - 2*Li(4,y)*a2^2*f1 - 4*Li(4,y)*a1*e1*f1 + 6* Li(4,y)*a1*e1^2 + 4*Li(4,y)*a1*a3*f1 - 6*Li(4,y)*a1*a3*e1 + 4* Li(4,y)*a1*a2*f1 - 6*Li(4,y)*a1*a2*e1 + 6*Li(4,y)*a1*a2*a3 - 2 *Li(4,y)*a1^2*f1 ) + ep^4*F32B(a1,a2,a3,f1,e1,z,ep^4) ) ; .sort ***************************************************************************************** * F32C = Hypergeometri3F2[1+a1*ep,1+a2*ep, 2+d1*ep; 3/2+f1*ep, 2+c1*ep;z] **************************************************************************************** id F32C(a1?,a2?,d1?,f1?,c1?,z?,ep?) = (-2*y/[1-y]^2)*num(1,c1,ep,1)*den(1,d1,ep,1)*num(1,2*f1,ep,1)*[1-y]/[1+y]* ( + ln(y) + ep * ( - 2*ln(y)*ln(1 - y)*f1 - 2*ln(y)*ln(1 + y)*a1 - 2*ln(y)*ln(1 + y)*a2 + 2*ln(y)*ln(1 + y)*c1 - 2*ln(y)*ln(1 + y)*d1 + 2*ln(y)*ln(1 + y)*f1 + 1/2*ln(y)*ln(y)*a1 + 1/2*ln(y)*ln(y)*a2 + ln(y)*ln(y)*c1* [1-y]^-1 - ln(y)*ln(y)*c1 - ln(y)*ln(y)*d1*[1-y]^-1 + ln(y)*ln(y)*d1 - zt(2)*a1 - zt(2)*a2 + zt(2)*c1 - zt(2)*d1 + 3*zt(2)*f1 - 2*Li(2, - y)*a1 - 2*Li(2, - y)*a2 + 2*Li(2, - y)*c1 - 2*Li(2, - y)*d1 + 2* Li(2, - y)*f1 - 2*Li(2,y)*f1 ) + ep^2 * ( 8*S(1,2, - y)*a1*a2 - 8*S(1,2, - y)*a1*c1 + 8*S(1,2, - y)*a1 *d1 - 12*S(1,2, - y)*a1*f1 + 4*S(1,2, - y)*a1^2 - 8*S(1,2, - y)*a2*c1 + 8*S(1,2, - y)*a2*d1 - 12*S(1,2, - y)*a2*f1 + 4*S(1,2, - y)*a2^2 - 8*S(1,2, - y)*c1*d1 + 12*S(1,2, - y)*c1*f1 + 4*S(1,2, - y)*c1^2 - 12* S(1,2, - y)*d1*f1 + 4*S(1,2, - y)*d1^2 + 8*S(1,2, - y)*f1^2 + 2*S(1,2 ,y^2)*a1*f1 + 2*S(1,2,y^2)*a2*f1 - 2*S(1,2,y^2)*c1*f1 + 2*S(1,2,y^2)* d1*f1 - 2*S(1,2,y^2)*f1^2 - 4*S(1,2,y)*a1*f1 - 4*S(1,2,y)*a2*f1 + 4* S(1,2,y)*c1*f1 - 4*S(1,2,y)*d1*f1 + 8*S(1,2,y)*f1^2 + 2*ln(1 - y)*zt( 2)*a1*f1 + 2*ln(1 - y)*zt(2)*a2*f1 - 2*ln(1 - y)*zt(2)*c1*f1 + 2*ln(1 - y)*zt(2)*d1*f1 - 6*ln(1 - y)*zt(2)*f1^2 + 4*ln(1 - y)*Li(2, - y)* a1*f1 + 4*ln(1 - y)*Li(2, - y)*a2*f1 - 4*ln(1 - y)*Li(2, - y)*c1*f1 + 4*ln(1 - y)*Li(2, - y)*d1*f1 - 4*ln(1 - y)*Li(2, - y)*f1^2 + 4*ln( 1 - y)*Li(2,y)*f1^2 + 4*ln(1 + y)*zt(2)*a1*a2 - 4*ln(1 + y)*zt(2)*a1* c1 + 4*ln(1 + y)*zt(2)*a1*d1 - 8*ln(1 + y)*zt(2)*a1*f1 + 2*ln(1 + y)* zt(2)*a1^2 - 4*ln(1 + y)*zt(2)*a2*c1 + 4*ln(1 + y)*zt(2)*a2*d1 - 8* ln(1 + y)*zt(2)*a2*f1 + 2*ln(1 + y)*zt(2)*a2^2 - 4*ln(1 + y)*zt(2)*c1 *d1 + 8*ln(1 + y)*zt(2)*c1*f1 + 2*ln(1 + y)*zt(2)*c1^2 - 8*ln(1 + y)* zt(2)*d1*f1 + 2*ln(1 + y)*zt(2)*d1^2 + 6*ln(1 + y)*zt(2)*f1^2 + 8*ln( 1 + y)*Li(2, - y)*a1*a2 - 8*ln(1 + y)*Li(2, - y)*a1*c1 + 8*ln(1 + y)* Li(2, - y)*a1*d1 - 8*ln(1 + y)*Li(2, - y)*a1*f1 + 4*ln(1 + y)*Li(2, - y)*a1^2 - 8*ln(1 + y)*Li(2, - y)*a2*c1 + 8*ln(1 + y)*Li(2, - y)*a2 *d1 - 8*ln(1 + y)*Li(2, - y)*a2*f1 + 4*ln(1 + y)*Li(2, - y)*a2^2 - 8* ln(1 + y)*Li(2, - y)*c1*d1 + 8*ln(1 + y)*Li(2, - y)*c1*f1 + 4*ln(1 + y)*Li(2, - y)*c1^2 - 8*ln(1 + y)*Li(2, - y)*d1*f1 + 4*ln(1 + y)*Li(2, - y)*d1^2 + 4*ln(1 + y)*Li(2, - y)*f1^2 + 4*ln(1 + y)*Li(2,y)*a1*f1 + 4*ln(1 + y)*Li(2,y)*a2*f1 - 4*ln(1 + y)*Li(2,y)*c1*f1 + 4*ln(1 + y )*Li(2,y)*d1*f1 - 4*ln(1 + y)*Li(2,y)*f1^2 + 2*ln(y)*ln(1 - y)*ln(1 - y)*f1^2 + 4*ln(y)*ln(1 - y)*ln(1 + y)*a1*f1 + 4*ln(y)*ln(1 - y)* ln(1 + y)*a2*f1 - 4*ln(y)*ln(1 - y)*ln(1 + y)*c1*f1 + 4*ln(y)*ln(1 - y)*ln(1 + y)*d1*f1 - 4*ln(y)*ln(1 - y)*ln(1 + y)*f1^2 + 4*ln(y)*ln(1 + y)*ln(1 + y)*a1*a2 - 4*ln(y)*ln(1 + y)*ln(1 + y)*a1*c1 + 4*ln(y)* ln(1 + y)*ln(1 + y)*a1*d1 - 4*ln(y)*ln(1 + y)*ln(1 + y)*a1*f1 + 2*ln( y)*ln(1 + y)*ln(1 + y)*a1^2 - 4*ln(y)*ln(1 + y)*ln(1 + y)*a2*c1 + 4* ln(y)*ln(1 + y)*ln(1 + y)*a2*d1 - 4*ln(y)*ln(1 + y)*ln(1 + y)*a2*f1 + 2*ln(y)*ln(1 + y)*ln(1 + y)*a2^2 - 4*ln(y)*ln(1 + y)*ln(1 + y)*c1* d1 + 4*ln(y)*ln(1 + y)*ln(1 + y)*c1*f1 + 2*ln(y)*ln(1 + y)*ln(1 + y)* c1^2 - 4*ln(y)*ln(1 + y)*ln(1 + y)*d1*f1 + 2*ln(y)*ln(1 + y)*ln(1 + y )*d1^2 + 2*ln(y)*ln(1 + y)*ln(1 + y)*f1^2 - ln(y)*ln(y)*ln(1 - y)*a1* f1 - ln(y)*ln(y)*ln(1 - y)*a2*f1 + 2*ln(y)*ln(y)*ln(1 - y)*c1*d1* [1-y]^-1 - ln(y)*ln(y)*ln(1 - y)*c1*d1 + ln(y)*ln(y)*ln(1 - y)*c1*f1 - 2*ln(y)*ln(y)*ln(1 - y)*c1^2*[1-y]^-1 + ln(y)*ln(y)*ln(1 - y)*c1^2 - ln(y)*ln(y)*ln(1 - y)*d1*f1 - 2*ln(y)*ln(y)*ln(1 + y)*a1*a2 + 2* ln(y)*ln(y)*ln(1 + y)*a1*c1 - 2*ln(y)*ln(y)*ln(1 + y)*a1*d1 + ln(y)* ln(y)*ln(1 + y)*a1*f1 - ln(y)*ln(y)*ln(1 + y)*a1^2 + 2*ln(y)*ln(y)* ln(1 + y)*a2*c1 - 2*ln(y)*ln(y)*ln(1 + y)*a2*d1 + ln(y)*ln(y)*ln(1 + y)*a2*f1 - ln(y)*ln(y)*ln(1 + y)*a2^2 + 2*ln(y)*ln(y)*ln(1 + y)*c1*d1 - ln(y)*ln(y)*ln(1 + y)*c1*f1 - ln(y)*ln(y)*ln(1 + y)*c1^2 + ln(y)* ln(y)*ln(1 + y)*d1*f1 - ln(y)*ln(y)*ln(1 + y)*d1^2 + 1/6*ln(y)*ln(y)* ln(y)*a1*a2 + 1/3*ln(y)*ln(y)*ln(y)*a1*c1*[1-y]^-1 - 1/3*ln(y)*ln(y)* ln(y)*a1*c1 - 1/3*ln(y)*ln(y)*ln(y)*a1*d1*[1-y]^-1 + 1/3*ln(y)*ln(y)* ln(y)*a1*d1 + 1/6*ln(y)*ln(y)*ln(y)*a1^2 + 1/3*ln(y)*ln(y)*ln(y)*a2* c1*[1-y]^-1 - 1/3*ln(y)*ln(y)*ln(y)*a2*c1 - 1/3*ln(y)*ln(y)*ln(y)*a2* d1*[1-y]^-1 + 1/3*ln(y)*ln(y)*ln(y)*a2*d1 + 1/6*ln(y)*ln(y)*ln(y)* a2^2 + 1/3*ln(y)*ln(y)*ln(y)*c1*d1*[1-y]^-1 - 1/3*ln(y)*ln(y)*ln(y)* c1*d1 - 1/3*ln(y)*ln(y)*ln(y)*d1^2*[1-y]^-1 + 1/3*ln(y)*ln(y)*ln(y)* d1^2 - 2*ln(y)*zt(2)*a1*a2 - 2*ln(y)*zt(2)*a1*c1*[1-y]^-1 + 3*ln(y)* zt(2)*a1*c1 + 2*ln(y)*zt(2)*a1*d1*[1-y]^-1 - 3*ln(y)*zt(2)*a1*d1 + 3* ln(y)*zt(2)*a1*f1 - ln(y)*zt(2)*a1^2 - 2*ln(y)*zt(2)*a2*c1*[1-y]^-1 + 3*ln(y)*zt(2)*a2*c1 + 2*ln(y)*zt(2)*a2*d1*[1-y]^-1 - 3*ln(y)*zt(2) *a2*d1 + 3*ln(y)*zt(2)*a2*f1 - ln(y)*zt(2)*a2^2 - 4*ln(y)*zt(2)*c1*d1 *[1-y]^-1 + 4*ln(y)*zt(2)*c1*d1 + 6*ln(y)*zt(2)*c1*f1*[1-y]^-1 - 6* ln(y)*zt(2)*c1*f1 + 2*ln(y)*zt(2)*c1^2*[1-y]^-1 - 2*ln(y)*zt(2)*c1^2 - 6*ln(y)*zt(2)*d1*f1*[1-y]^-1 + 6*ln(y)*zt(2)*d1*f1 + 2*ln(y)*zt(2) *d1^2*[1-y]^-1 - 2*ln(y)*zt(2)*d1^2 + 4*ln(y)*Li(2, - y)*a1*c1* [1-y]^-1 - 2*ln(y)*Li(2, - y)*a1*c1 - 4*ln(y)*Li(2, - y)*a1*d1* [1-y]^-1 + 2*ln(y)*Li(2, - y)*a1*d1 + 4*ln(y)*Li(2, - y)*a2*c1* [1-y]^-1 - 2*ln(y)*Li(2, - y)*a2*c1 - 4*ln(y)*Li(2, - y)*a2*d1* [1-y]^-1 + 2*ln(y)*Li(2, - y)*a2*d1 + 8*ln(y)*Li(2, - y)*c1*d1* [1-y]^-1 - 4*ln(y)*Li(2, - y)*c1*d1 - 4*ln(y)*Li(2, - y)*c1*f1* [1-y]^-1 + 2*ln(y)*Li(2, - y)*c1*f1 - 4*ln(y)*Li(2, - y)*c1^2* [1-y]^-1 + 2*ln(y)*Li(2, - y)*c1^2 + 4*ln(y)*Li(2, - y)*d1*f1* [1-y]^-1 - 2*ln(y)*Li(2, - y)*d1*f1 - 4*ln(y)*Li(2, - y)*d1^2* [1-y]^-1 + 2*ln(y)*Li(2, - y)*d1^2 + 4*ln(y)*Li(2,y)*c1*d1*[1-y]^-1 - 2*ln(y)*Li(2,y)*c1*d1 + 4*ln(y)*Li(2,y)*c1*f1*[1-y]^-1 - 2*ln(y)* Li(2,y)*c1*f1 - 4*ln(y)*Li(2,y)*c1^2*[1-y]^-1 + 2*ln(y)*Li(2,y)*c1^2 - 4*ln(y)*Li(2,y)*d1*f1*[1-y]^-1 + 2*ln(y)*Li(2,y)*d1*f1 - 4*zt(3)* a1*a2 - 6*zt(3)*a1*c1*[1-y]^-1 + 7*zt(3)*a1*c1 + 6*zt(3)*a1*d1* [1-y]^-1 - 7*zt(3)*a1*d1 + 7*zt(3)*a1*f1 - 2*zt(3)*a1^2 - 6*zt(3)*a2* c1*[1-y]^-1 + 7*zt(3)*a2*c1 + 6*zt(3)*a2*d1*[1-y]^-1 - 7*zt(3)*a2*d1 + 7*zt(3)*a2*f1 - 2*zt(3)*a2^2 - 8*zt(3)*c1*d1*[1-y]^-1 + 8*zt(3)*c1 *d1 + 14*zt(3)*c1*f1*[1-y]^-1 - 14*zt(3)*c1*f1 + 2*zt(3)*c1^2* [1-y]^-1 - 3*zt(3)*c1^2 - 14*zt(3)*d1*f1*[1-y]^-1 + 14*zt(3)*d1*f1 + 6*zt(3)*d1^2*[1-y]^-1 - 5*zt(3)*d1^2 - 7*zt(3)*f1^2 - 4*Li(3, - y)*a1 *a2 - 8*Li(3, - y)*a1*c1*[1-y]^-1 + 8*Li(3, - y)*a1*c1 + 8*Li(3, - y) *a1*d1*[1-y]^-1 - 8*Li(3, - y)*a1*d1 + 2*Li(3, - y)*a1*f1 - 2*Li(3, - y)*a1^2 - 8*Li(3, - y)*a2*c1*[1-y]^-1 + 8*Li(3, - y)*a2*c1 + 8*Li( 3, - y)*a2*d1*[1-y]^-1 - 8*Li(3, - y)*a2*d1 + 2*Li(3, - y)*a2*f1 - 2* Li(3, - y)*a2^2 - 16*Li(3, - y)*c1*d1*[1-y]^-1 + 12*Li(3, - y)*c1*d1 + 8*Li(3, - y)*c1*f1*[1-y]^-1 - 6*Li(3, - y)*c1*f1 + 8*Li(3, - y)* c1^2*[1-y]^-1 - 6*Li(3, - y)*c1^2 - 8*Li(3, - y)*d1*f1*[1-y]^-1 + 6* Li(3, - y)*d1*f1 + 8*Li(3, - y)*d1^2*[1-y]^-1 - 6*Li(3, - y)*d1^2 - 2 *Li(3,y)*a1*f1 - 2*Li(3,y)*a2*f1 - 4*Li(3,y)*c1*d1*[1-y]^-1 + 2*Li(3, y)*c1*d1 - 8*Li(3,y)*c1*f1*[1-y]^-1 + 6*Li(3,y)*c1*f1 + 4*Li(3,y)* c1^2*[1-y]^-1 - 2*Li(3,y)*c1^2 + 8*Li(3,y)*d1*f1*[1-y]^-1 - 6*Li(3,y) *d1*f1 ) + ep^3 * ( 48*S(1,2, - y)*ln(1 + y)*a1*a2*c1 - 48*S(1,2, - y)*ln(1 + y) *a1*a2*d1 + 48*S(1,2, - y)*ln(1 + y)*a1*a2*f1 - 24*S(1,2, - y)*ln(1 + y)*a1*a2^2 + 48*S(1,2, - y)*ln(1 + y)*a1*c1*d1 - 48*S(1,2, - y)* ln(1 + y)*a1*c1*f1 - 24*S(1,2, - y)*ln(1 + y)*a1*c1^2 + 48*S(1,2, - y )*ln(1 + y)*a1*d1*f1 - 24*S(1,2, - y)*ln(1 + y)*a1*d1^2 - 16*S(1,2, - y)*ln(1 + y)*a1*f1^2 - 24*S(1,2, - y)*ln(1 + y)*a1^2*a2 + 24*S(1,2 , - y)*ln(1 + y)*a1^2*c1 - 24*S(1,2, - y)*ln(1 + y)*a1^2*d1 + 24*S(1, 2, - y)*ln(1 + y)*a1^2*f1 - 8*S(1,2, - y)*ln(1 + y)*a1^3 + 48*S(1,2, - y)*ln(1 + y)*a2*c1*d1 - 48*S(1,2, - y)*ln(1 + y)*a2*c1*f1 - 24*S(1 ,2, - y)*ln(1 + y)*a2*c1^2 + 48*S(1,2, - y)*ln(1 + y)*a2*d1*f1 - 24* S(1,2, - y)*ln(1 + y)*a2*d1^2 - 16*S(1,2, - y)*ln(1 + y)*a2*f1^2 + 24 *S(1,2, - y)*ln(1 + y)*a2^2*c1 - 24*S(1,2, - y)*ln(1 + y)*a2^2*d1 + 24*S(1,2, - y)*ln(1 + y)*a2^2*f1 - 8*S(1,2, - y)*ln(1 + y)*a2^3 - 48* S(1,2, - y)*ln(1 + y)*c1*d1*f1 + 24*S(1,2, - y)*ln(1 + y)*c1*d1^2 + 16*S(1,2, - y)*ln(1 + y)*c1*f1^2 - 24*S(1,2, - y)*ln(1 + y)*c1^2*d1 + 24*S(1,2, - y)*ln(1 + y)*c1^2*f1 + 8*S(1,2, - y)*ln(1 + y)*c1^3 - 16*S(1,2, - y)*ln(1 + y)*d1*f1^2 + 24*S(1,2, - y)*ln(1 + y)*d1^2*f1 - 8*S(1,2, - y)*ln(1 + y)*d1^3 + 16*S(1,2, - y)*ln(y)*a1*a2*c1* [1-y]^-1 - 8*S(1,2, - y)*ln(y)*a1*a2*c1 - 16*S(1,2, - y)*ln(y)*a1*a2* d1*[1-y]^-1 + 8*S(1,2, - y)*ln(y)*a1*a2*d1 + 24*S(1,2, - y)*ln(y)*a1* c1*d1*[1-y]^-1 - 12*S(1,2, - y)*ln(y)*a1*c1*d1 - 24*S(1,2, - y)*ln(y) *a1*c1*f1*[1-y]^-1 + 12*S(1,2, - y)*ln(y)*a1*c1*f1 - 8*S(1,2, - y)* ln(y)*a1*c1^2*[1-y]^-1 + 4*S(1,2, - y)*ln(y)*a1*c1^2 + 24*S(1,2, - y) *ln(y)*a1*d1*f1*[1-y]^-1 - 12*S(1,2, - y)*ln(y)*a1*d1*f1 - 16*S(1,2, - y)*ln(y)*a1*d1^2*[1-y]^-1 + 8*S(1,2, - y)*ln(y)*a1*d1^2 + 8*S(1,2, - y)*ln(y)*a1^2*c1*[1-y]^-1 - 4*S(1,2, - y)*ln(y)*a1^2*c1 - 8*S(1,2, - y)*ln(y)*a1^2*d1*[1-y]^-1 + 4*S(1,2, - y)*ln(y)*a1^2*d1 + 24*S(1,2 , - y)*ln(y)*a2*c1*d1*[1-y]^-1 - 12*S(1,2, - y)*ln(y)*a2*c1*d1 - 24* S(1,2, - y)*ln(y)*a2*c1*f1*[1-y]^-1 + 12*S(1,2, - y)*ln(y)*a2*c1*f1 - 8*S(1,2, - y)*ln(y)*a2*c1^2*[1-y]^-1 + 4*S(1,2, - y)*ln(y)*a2*c1^2 + 24*S(1,2, - y)*ln(y)*a2*d1*f1*[1-y]^-1 - 12*S(1,2, - y)*ln(y)*a2* d1*f1 - 16*S(1,2, - y)*ln(y)*a2*d1^2*[1-y]^-1 + 8*S(1,2, - y)*ln(y)* a2*d1^2 + 8*S(1,2, - y)*ln(y)*a2^2*c1*[1-y]^-1 - 4*S(1,2, - y)*ln(y)* a2^2*c1 - 8*S(1,2, - y)*ln(y)*a2^2*d1*[1-y]^-1 + 4*S(1,2, - y)*ln(y)* a2^2*d1 - 40*S(1,2, - y)*ln(y)*c1*d1*f1*[1-y]^-1 + 20*S(1,2, - y)*ln( y)*c1*d1*f1 + 16*S(1,2, - y)*ln(y)*c1*d1^2*[1-y]^-1 - 8*S(1,2, - y)* ln(y)*c1*d1^2 + 16*S(1,2, - y)*ln(y)*c1*f1^2*[1-y]^-1 - 8*S(1,2, - y) *ln(y)*c1*f1^2 - 8*S(1,2, - y)*ln(y)*c1^2*d1*[1-y]^-1 + 4*S(1,2, - y) *ln(y)*c1^2*d1 + 16*S(1,2, - y)*ln(y)*c1^2*f1*[1-y]^-1 - 8*S(1,2, - y )*ln(y)*c1^2*f1 - 16*S(1,2, - y)*ln(y)*d1*f1^2*[1-y]^-1 + 8*S(1,2, - y)*ln(y)*d1*f1^2 + 24*S(1,2, - y)*ln(y)*d1^2*f1*[1-y]^-1 - 12*S(1,2, - y)*ln(y)*d1^2*f1 - 8*S(1,2, - y)*ln(y)*d1^3*[1-y]^-1 + 4*S(1,2, - y)*ln(y)*d1^3 - 4*S(1,2,y^2)*ln(1 - y)*a1*a2*f1 + 4*S(1,2,y^2)*ln(1 - y)*a1*c1*f1 - 4*S(1,2,y^2)*ln(1 - y)*a1*d1*f1 + 2*S(1,2,y^2)*ln(1 - y)*a1*f1^2 - 2*S(1,2,y^2)*ln(1 - y)*a1^2*f1 + 4*S(1,2,y^2)*ln(1 - y)*a2*c1*f1 - 4*S(1,2,y^2)*ln(1 - y)*a2*d1*f1 + 2*S(1,2,y^2)*ln(1 - y )*a2*f1^2 - 2*S(1,2,y^2)*ln(1 - y)*a2^2*f1 + 4*S(1,2,y^2)*ln(1 - y)* c1*d1*f1 - 2*S(1,2,y^2)*ln(1 - y)*c1*f1^2 - 2*S(1,2,y^2)*ln(1 - y)* c1^2*f1 + 2*S(1,2,y^2)*ln(1 - y)*d1*f1^2 - 2*S(1,2,y^2)*ln(1 - y)* d1^2*f1 - 4*S(1,2,y^2)*ln(1 + y)*a1*a2*f1 + 4*S(1,2,y^2)*ln(1 + y)*a1 *c1*f1 - 4*S(1,2,y^2)*ln(1 + y)*a1*d1*f1 + 2*S(1,2,y^2)*ln(1 + y)*a1* f1^2 - 2*S(1,2,y^2)*ln(1 + y)*a1^2*f1 + 4*S(1,2,y^2)*ln(1 + y)*a2*c1* f1 - 4*S(1,2,y^2)*ln(1 + y)*a2*d1*f1 + 2*S(1,2,y^2)*ln(1 + y)*a2*f1^2 - 2*S(1,2,y^2)*ln(1 + y)*a2^2*f1 + 4*S(1,2,y^2)*ln(1 + y)*c1*d1*f1 - 2*S(1,2,y^2)*ln(1 + y)*c1*f1^2 - 2*S(1,2,y^2)*ln(1 + y)*c1^2*f1 + 2*S(1,2,y^2)*ln(1 + y)*d1*f1^2 - 2*S(1,2,y^2)*ln(1 + y)*d1^2*f1 + 4* S(1,2,y^2)*ln(y)*a1*c1*d1*[1-y]^-1 - 2*S(1,2,y^2)*ln(y)*a1*c1*d1 + 4* S(1,2,y^2)*ln(y)*a1*c1*f1*[1-y]^-1 - 2*S(1,2,y^2)*ln(y)*a1*c1*f1 - 4* S(1,2,y^2)*ln(y)*a1*c1^2*[1-y]^-1 + 2*S(1,2,y^2)*ln(y)*a1*c1^2 - 4*S( 1,2,y^2)*ln(y)*a1*d1*f1*[1-y]^-1 + 2*S(1,2,y^2)*ln(y)*a1*d1*f1 + 4*S( 1,2,y^2)*ln(y)*a2*c1*d1*[1-y]^-1 - 2*S(1,2,y^2)*ln(y)*a2*c1*d1 + 4*S( 1,2,y^2)*ln(y)*a2*c1*f1*[1-y]^-1 - 2*S(1,2,y^2)*ln(y)*a2*c1*f1 - 4*S( 1,2,y^2)*ln(y)*a2*c1^2*[1-y]^-1 + 2*S(1,2,y^2)*ln(y)*a2*c1^2 - 4*S(1, 2,y^2)*ln(y)*a2*d1*f1*[1-y]^-1 + 2*S(1,2,y^2)*ln(y)*a2*d1*f1 + 4*S(1, 2,y^2)*ln(y)*c1*d1*f1*[1-y]^-1 - 2*S(1,2,y^2)*ln(y)*c1*d1*f1 + 4*S(1, 2,y^2)*ln(y)*c1*d1^2*[1-y]^-1 - 2*S(1,2,y^2)*ln(y)*c1*d1^2 - 4*S(1,2, y^2)*ln(y)*c1*f1^2*[1-y]^-1 + 2*S(1,2,y^2)*ln(y)*c1*f1^2 - 8*S(1,2, y^2)*ln(y)*c1^2*d1*[1-y]^-1 + 4*S(1,2,y^2)*ln(y)*c1^2*d1 + 4*S(1,2, y^2)*ln(y)*c1^3*[1-y]^-1 - 2*S(1,2,y^2)*ln(y)*c1^3 + 4*S(1,2,y^2)*ln( y)*d1*f1^2*[1-y]^-1 - 2*S(1,2,y^2)*ln(y)*d1*f1^2 - 4*S(1,2,y^2)*ln(y) *d1^2*f1*[1-y]^-1 + 2*S(1,2,y^2)*ln(y)*d1^2*f1 + 16*S(1,2,y)*ln(1 - y )*a1*a2*f1 - 16*S(1,2,y)*ln(1 - y)*a1*c1*f1 + 16*S(1,2,y)*ln(1 - y)* a1*d1*f1 - 16*S(1,2,y)*ln(1 - y)*a1*f1^2 + 8*S(1,2,y)*ln(1 - y)*a1^2* f1 - 16*S(1,2,y)*ln(1 - y)*a2*c1*f1 + 16*S(1,2,y)*ln(1 - y)*a2*d1*f1 - 16*S(1,2,y)*ln(1 - y)*a2*f1^2 + 8*S(1,2,y)*ln(1 - y)*a2^2*f1 - 16* S(1,2,y)*ln(1 - y)*c1*d1*f1 + 16*S(1,2,y)*ln(1 - y)*c1*f1^2 + 8*S(1,2 ,y)*ln(1 - y)*c1^2*f1 - 16*S(1,2,y)*ln(1 - y)*d1*f1^2 + 8*S(1,2,y)* ln(1 - y)*d1^2*f1 - 8*S(1,2,y)*ln(y)*a1*c1*d1*[1-y]^-1 + 4*S(1,2,y)* ln(y)*a1*c1*d1 - 8*S(1,2,y)*ln(y)*a1*c1*f1*[1-y]^-1 + 4*S(1,2,y)*ln(y )*a1*c1*f1 + 8*S(1,2,y)*ln(y)*a1*c1^2*[1-y]^-1 - 4*S(1,2,y)*ln(y)*a1* c1^2 + 8*S(1,2,y)*ln(y)*a1*d1*f1*[1-y]^-1 - 4*S(1,2,y)*ln(y)*a1*d1*f1 - 8*S(1,2,y)*ln(y)*a2*c1*d1*[1-y]^-1 + 4*S(1,2,y)*ln(y)*a2*c1*d1 - 8 *S(1,2,y)*ln(y)*a2*c1*f1*[1-y]^-1 + 4*S(1,2,y)*ln(y)*a2*c1*f1 + 8*S(1 ,2,y)*ln(y)*a2*c1^2*[1-y]^-1 - 4*S(1,2,y)*ln(y)*a2*c1^2 + 8*S(1,2,y)* ln(y)*a2*d1*f1*[1-y]^-1 - 4*S(1,2,y)*ln(y)*a2*d1*f1 + 8*S(1,2,y)*ln(y )*c1*d1*f1*[1-y]^-1 - 4*S(1,2,y)*ln(y)*c1*d1*f1 - 8*S(1,2,y)*ln(y)*c1 *d1^2*[1-y]^-1 + 4*S(1,2,y)*ln(y)*c1*d1^2 + 16*S(1,2,y)*ln(y)*c1*f1^2 *[1-y]^-1 - 8*S(1,2,y)*ln(y)*c1*f1^2 + 8*S(1,2,y)*ln(y)*c1^2*d1* [1-y]^-1 - 4*S(1,2,y)*ln(y)*c1^2*d1 - 16*S(1,2,y)*ln(y)*c1^2*f1* [1-y]^-1 + 8*S(1,2,y)*ln(y)*c1^2*f1 - 16*S(1,2,y)*ln(y)*d1*f1^2* [1-y]^-1 + 8*S(1,2,y)*ln(y)*d1*f1^2 + 8*S(1,2,y)*ln(y)*d1^2*f1* [1-y]^-1 - 4*S(1,2,y)*ln(y)*d1^2*f1 + 48*S(1,3, - y)*a1*a2*c1 - 48*S( 1,3, - y)*a1*a2*d1 + 48*S(1,3, - y)*a1*a2*f1 - 24*S(1,3, - y)*a1*a2^2 + 48*S(1,3, - y)*a1*c1*d1 - 48*S(1,3, - y)*a1*c1*f1 - 24*S(1,3, - y) *a1*c1^2 + 48*S(1,3, - y)*a1*d1*f1 - 24*S(1,3, - y)*a1*d1^2 - 16*S(1, 3, - y)*a1*f1^2 - 24*S(1,3, - y)*a1^2*a2 + 24*S(1,3, - y)*a1^2*c1 - 24*S(1,3, - y)*a1^2*d1 + 24*S(1,3, - y)*a1^2*f1 - 8*S(1,3, - y)*a1^3 + 48*S(1,3, - y)*a2*c1*d1 - 48*S(1,3, - y)*a2*c1*f1 - 24*S(1,3, - y) *a2*c1^2 + 48*S(1,3, - y)*a2*d1*f1 - 24*S(1,3, - y)*a2*d1^2 - 16*S(1, 3, - y)*a2*f1^2 + 24*S(1,3, - y)*a2^2*c1 - 24*S(1,3, - y)*a2^2*d1 + 24*S(1,3, - y)*a2^2*f1 - 8*S(1,3, - y)*a2^3 - 48*S(1,3, - y)*c1*d1*f1 + 24*S(1,3, - y)*c1*d1^2 + 16*S(1,3, - y)*c1*f1^2 - 24*S(1,3, - y)* c1^2*d1 + 24*S(1,3, - y)*c1^2*f1 + 8*S(1,3, - y)*c1^3 - 16*S(1,3, - y )*d1*f1^2 + 24*S(1,3, - y)*d1^2*f1 - 8*S(1,3, - y)*d1^3 - 4*S(1,3,y^2 )*a1*a2*f1 + 4*S(1,3,y^2)*a1*c1*f1 - 4*S(1,3,y^2)*a1*d1*f1 + 2*S(1,3, y^2)*a1*f1^2 - 2*S(1,3,y^2)*a1^2*f1 + 4*S(1,3,y^2)*a2*c1*f1 - 4*S(1,3 ,y^2)*a2*d1*f1 + 2*S(1,3,y^2)*a2*f1^2 - 2*S(1,3,y^2)*a2^2*f1 + 4*S(1, 3,y^2)*c1*d1*f1 - 2*S(1,3,y^2)*c1*f1^2 - 2*S(1,3,y^2)*c1^2*f1 + 2*S(1 ,3,y^2)*d1*f1^2 - 2*S(1,3,y^2)*d1^2*f1 + 16*S(1,3,y)*a1*a2*f1 - 16*S( 1,3,y)*a1*c1*f1 + 16*S(1,3,y)*a1*d1*f1 - 16*S(1,3,y)*a1*f1^2 + 8*S(1, 3,y)*a1^2*f1 - 16*S(1,3,y)*a2*c1*f1 + 16*S(1,3,y)*a2*d1*f1 - 16*S(1,3 ,y)*a2*f1^2 + 8*S(1,3,y)*a2^2*f1 - 16*S(1,3,y)*c1*d1*f1 + 16*S(1,3,y) *c1*f1^2 + 8*S(1,3,y)*c1^2*f1 - 16*S(1,3,y)*d1*f1^2 + 8*S(1,3,y)*d1^2 *f1 - 24*S(2,2, - y)*a1*a2*c1 + 24*S(2,2, - y)*a1*a2*d1 - 24*S(2,2, - y)*a1*a2*f1 + 12*S(2,2, - y)*a1*a2^2 + 8*S(2,2, - y)*a1*c1*d1* [1-y]^-1 - 28*S(2,2, - y)*a1*c1*d1 + 24*S(2,2, - y)*a1*c1*f1 - 8*S(2, 2, - y)*a1*c1^2*[1-y]^-1 + 16*S(2,2, - y)*a1*c1^2 - 24*S(2,2, - y)*a1 *d1*f1 + 12*S(2,2, - y)*a1*d1^2 + 8*S(2,2, - y)*a1*f1^2 + 12*S(2,2, - y)*a1^2*a2 - 12*S(2,2, - y)*a1^2*c1 + 12*S(2,2, - y)*a1^2*d1 - 12* S(2,2, - y)*a1^2*f1 + 4*S(2,2, - y)*a1^3 + 8*S(2,2, - y)*a2*c1*d1* [1-y]^-1 - 28*S(2,2, - y)*a2*c1*d1 + 24*S(2,2, - y)*a2*c1*f1 - 8*S(2, 2, - y)*a2*c1^2*[1-y]^-1 + 16*S(2,2, - y)*a2*c1^2 - 24*S(2,2, - y)*a2 *d1*f1 + 12*S(2,2, - y)*a2*d1^2 + 8*S(2,2, - y)*a2*f1^2 - 12*S(2,2, - y)*a2^2*c1 + 12*S(2,2, - y)*a2^2*d1 - 12*S(2,2, - y)*a2^2*f1 + 4* S(2,2, - y)*a2^3 - 8*S(2,2, - y)*c1*d1*f1*[1-y]^-1 + 28*S(2,2, - y)* c1*d1*f1 + 8*S(2,2, - y)*c1*d1^2*[1-y]^-1 - 16*S(2,2, - y)*c1*d1^2 - 8*S(2,2, - y)*c1*f1^2 - 16*S(2,2, - y)*c1^2*d1*[1-y]^-1 + 20*S(2,2, - y)*c1^2*d1 + 8*S(2,2, - y)*c1^2*f1*[1-y]^-1 - 16*S(2,2, - y)*c1^2* f1 + 8*S(2,2, - y)*c1^3*[1-y]^-1 - 8*S(2,2, - y)*c1^3 + 8*S(2,2, - y) *d1*f1^2 - 12*S(2,2, - y)*d1^2*f1 + 4*S(2,2, - y)*d1^3 + 2*S(2,2,y^2) *a1*a2*f1 - 2*S(2,2,y^2)*a1*c1*d1*[1-y]^-1 + S(2,2,y^2)*a1*c1*d1 - 2* S(2,2,y^2)*a1*c1*f1 + 2*S(2,2,y^2)*a1*c1^2*[1-y]^-1 - S(2,2,y^2)*a1* c1^2 + 2*S(2,2,y^2)*a1*d1*f1 - S(2,2,y^2)*a1*f1^2 + S(2,2,y^2)*a1^2* f1 - 2*S(2,2,y^2)*a2*c1*d1*[1-y]^-1 + S(2,2,y^2)*a2*c1*d1 - 2*S(2,2, y^2)*a2*c1*f1 + 2*S(2,2,y^2)*a2*c1^2*[1-y]^-1 - S(2,2,y^2)*a2*c1^2 + 2*S(2,2,y^2)*a2*d1*f1 - S(2,2,y^2)*a2*f1^2 + S(2,2,y^2)*a2^2*f1 + 2* S(2,2,y^2)*c1*d1*f1*[1-y]^-1 - 3*S(2,2,y^2)*c1*d1*f1 - 2*S(2,2,y^2)* c1*d1^2*[1-y]^-1 + S(2,2,y^2)*c1*d1^2 + S(2,2,y^2)*c1*f1^2 + 4*S(2,2, y^2)*c1^2*d1*[1-y]^-1 - 2*S(2,2,y^2)*c1^2*d1 - 2*S(2,2,y^2)*c1^2*f1* [1-y]^-1 + 2*S(2,2,y^2)*c1^2*f1 - 2*S(2,2,y^2)*c1^3*[1-y]^-1 + S(2,2, y^2)*c1^3 - S(2,2,y^2)*d1*f1^2 + S(2,2,y^2)*d1^2*f1 - 8*S(2,2,y)*a1* a2*f1 + 8*S(2,2,y)*a1*c1*d1*[1-y]^-1 - 4*S(2,2,y)*a1*c1*d1 + 8*S(2,2, y)*a1*c1*f1 - 8*S(2,2,y)*a1*c1^2*[1-y]^-1 + 4*S(2,2,y)*a1*c1^2 - 8*S( 2,2,y)*a1*d1*f1 + 8*S(2,2,y)*a1*f1^2 - 4*S(2,2,y)*a1^2*f1 + 8*S(2,2,y )*a2*c1*d1*[1-y]^-1 - 4*S(2,2,y)*a2*c1*d1 + 8*S(2,2,y)*a2*c1*f1 - 8* S(2,2,y)*a2*c1^2*[1-y]^-1 + 4*S(2,2,y)*a2*c1^2 - 8*S(2,2,y)*a2*d1*f1 + 8*S(2,2,y)*a2*f1^2 - 4*S(2,2,y)*a2^2*f1 - 24*S(2,2,y)*c1*d1*f1* [1-y]^-1 + 20*S(2,2,y)*c1*d1*f1 + 8*S(2,2,y)*c1*d1^2*[1-y]^-1 - 4*S(2 ,2,y)*c1*d1^2 - 8*S(2,2,y)*c1*f1^2 - 8*S(2,2,y)*c1^2*d1*[1-y]^-1 + 4* S(2,2,y)*c1^2*d1 + 24*S(2,2,y)*c1^2*f1*[1-y]^-1 - 16*S(2,2,y)*c1^2*f1 + 8*S(2,2,y)*d1*f1^2 - 4*S(2,2,y)*d1^2*f1 - 8*H( - y)*a1*c1*d1* [1-y]^-1 + 4*H( - y)*a1*c1*d1 + 8*H( - y)*a1*c1^2*[1-y]^-1 - 4*H( - y )*a1*c1^2 - 8*H( - y)*a2*c1*d1*[1-y]^-1 + 4*H( - y)*a2*c1*d1 + 8*H( - y)*a2*c1^2*[1-y]^-1 - 4*H( - y)*a2*c1^2 + 8*H( - y)*c1*d1*f1* [1-y]^-1 - 4*H( - y)*c1*d1*f1 - 8*H( - y)*c1*d1^2*[1-y]^-1 + 4*H( - y )*c1*d1^2 + 16*H( - y)*c1^2*d1*[1-y]^-1 - 8*H( - y)*c1^2*d1 - 8*H( - y)*c1^2*f1*[1-y]^-1 + 4*H( - y)*c1^2*f1 - 8*H( - y)*c1^3*[1-y]^-1 + 4 *H( - y)*c1^3 - 6*ln(1 - y)*ln(1 - y)*zt(2)*a1*a2*f1 + 6*ln(1 - y)* ln(1 - y)*zt(2)*a1*c1*f1 - 6*ln(1 - y)*ln(1 - y)*zt(2)*a1*d1*f1 + 7* ln(1 - y)*ln(1 - y)*zt(2)*a1*f1^2 - 3*ln(1 - y)*ln(1 - y)*zt(2)*a1^2* f1 + 6*ln(1 - y)*ln(1 - y)*zt(2)*a2*c1*f1 - 6*ln(1 - y)*ln(1 - y)*zt( 2)*a2*d1*f1 + 7*ln(1 - y)*ln(1 - y)*zt(2)*a2*f1^2 - 3*ln(1 - y)*ln(1 - y)*zt(2)*a2^2*f1 + 6*ln(1 - y)*ln(1 - y)*zt(2)*c1*d1*f1 - 7*ln(1 - y)*ln(1 - y)*zt(2)*c1*f1^2 - 3*ln(1 - y)*ln(1 - y)*zt(2)*c1^2*f1 + 7*ln(1 - y)*ln(1 - y)*zt(2)*d1*f1^2 - 3*ln(1 - y)*ln(1 - y)*zt(2)* d1^2*f1 - 4*ln(1 - y)*ln(1 - y)*Li(2, - y)*a1*a2*f1 + 4*ln(1 - y)*ln( 1 - y)*Li(2, - y)*a1*c1*f1 - 4*ln(1 - y)*ln(1 - y)*Li(2, - y)*a1*d1* f1 + 2*ln(1 - y)*ln(1 - y)*Li(2, - y)*a1*f1^2 - 2*ln(1 - y)*ln(1 - y) *Li(2, - y)*a1^2*f1 + 4*ln(1 - y)*ln(1 - y)*Li(2, - y)*a2*c1*f1 - 4* ln(1 - y)*ln(1 - y)*Li(2, - y)*a2*d1*f1 + 2*ln(1 - y)*ln(1 - y)*Li(2, - y)*a2*f1^2 - 2*ln(1 - y)*ln(1 - y)*Li(2, - y)*a2^2*f1 + 4*ln(1 - y )*ln(1 - y)*Li(2, - y)*c1*d1*f1 - 2*ln(1 - y)*ln(1 - y)*Li(2, - y)*c1 *f1^2 - 2*ln(1 - y)*ln(1 - y)*Li(2, - y)*c1^2*f1 + 2*ln(1 - y)*ln(1 - y)*Li(2, - y)*d1*f1^2 - 2*ln(1 - y)*ln(1 - y)*Li(2, - y)*d1^2*f1 + 4*ln(1 - y)*ln(1 - y)*Li(2,y)*a1*a2*f1 - 4*ln(1 - y)*ln(1 - y)*Li( 2,y)*a1*c1*f1 + 4*ln(1 - y)*ln(1 - y)*Li(2,y)*a1*d1*f1 - 6*ln(1 - y)* ln(1 - y)*Li(2,y)*a1*f1^2 + 2*ln(1 - y)*ln(1 - y)*Li(2,y)*a1^2*f1 - 4 *ln(1 - y)*ln(1 - y)*Li(2,y)*a2*c1*f1 + 4*ln(1 - y)*ln(1 - y)*Li(2,y) *a2*d1*f1 - 6*ln(1 - y)*ln(1 - y)*Li(2,y)*a2*f1^2 + 2*ln(1 - y)*ln(1 - y)*Li(2,y)*a2^2*f1 - 4*ln(1 - y)*ln(1 - y)*Li(2,y)*c1*d1*f1 + 6* ln(1 - y)*ln(1 - y)*Li(2,y)*c1*f1^2 + 2*ln(1 - y)*ln(1 - y)*Li(2,y)* c1^2*f1 - 6*ln(1 - y)*ln(1 - y)*Li(2,y)*d1*f1^2 + 2*ln(1 - y)*ln(1 - y)*Li(2,y)*d1^2*f1 + 4*ln(1 - y)*ln(1 + y)*zt(2)*a1*a2*f1 - 4*ln(1 - y)*ln(1 + y)*zt(2)*a1*c1*f1 + 4*ln(1 - y)*ln(1 + y)*zt(2)*a1*d1*f1 - 2*ln(1 - y)*ln(1 + y)*zt(2)*a1*f1^2 + 2*ln(1 - y)*ln(1 + y)*zt(2)* a1^2*f1 - 4*ln(1 - y)*ln(1 + y)*zt(2)*a2*c1*f1 + 4*ln(1 - y)*ln(1 + y )*zt(2)*a2*d1*f1 - 2*ln(1 - y)*ln(1 + y)*zt(2)*a2*f1^2 + 2*ln(1 - y)* ln(1 + y)*zt(2)*a2^2*f1 - 4*ln(1 - y)*ln(1 + y)*zt(2)*c1*d1*f1 + 2* ln(1 - y)*ln(1 + y)*zt(2)*c1*f1^2 + 2*ln(1 - y)*ln(1 + y)*zt(2)*c1^2* f1 - 2*ln(1 - y)*ln(1 + y)*zt(2)*d1*f1^2 + 2*ln(1 - y)*ln(1 + y)*zt(2 )*d1^2*f1 - 8*ln(1 - y)*ln(1 + y)*Li(2, - y)*a1*a2*f1 + 8*ln(1 - y)* ln(1 + y)*Li(2, - y)*a1*c1*f1 - 8*ln(1 - y)*ln(1 + y)*Li(2, - y)*a1* d1*f1 + 4*ln(1 - y)*ln(1 + y)*Li(2, - y)*a1*f1^2 - 4*ln(1 - y)*ln(1 + y)*Li(2, - y)*a1^2*f1 + 8*ln(1 - y)*ln(1 + y)*Li(2, - y)*a2*c1*f1 - 8*ln(1 - y)*ln(1 + y)*Li(2, - y)*a2*d1*f1 + 4*ln(1 - y)*ln(1 + y)* Li(2, - y)*a2*f1^2 - 4*ln(1 - y)*ln(1 + y)*Li(2, - y)*a2^2*f1 + 8*ln( 1 - y)*ln(1 + y)*Li(2, - y)*c1*d1*f1 - 4*ln(1 - y)*ln(1 + y)*Li(2, - y)*c1*f1^2 - 4*ln(1 - y)*ln(1 + y)*Li(2, - y)*c1^2*f1 + 4*ln(1 - y)* ln(1 + y)*Li(2, - y)*d1*f1^2 - 4*ln(1 - y)*ln(1 + y)*Li(2, - y)*d1^2* f1 - 8*ln(1 - y)*ln(1 + y)*Li(2,y)*a1*a2*f1 + 8*ln(1 - y)*ln(1 + y)* Li(2,y)*a1*c1*f1 - 8*ln(1 - y)*ln(1 + y)*Li(2,y)*a1*d1*f1 + 4*ln(1 - y)*ln(1 + y)*Li(2,y)*a1*f1^2 - 4*ln(1 - y)*ln(1 + y)*Li(2,y)*a1^2*f1 + 8*ln(1 - y)*ln(1 + y)*Li(2,y)*a2*c1*f1 - 8*ln(1 - y)*ln(1 + y)*Li( 2,y)*a2*d1*f1 + 4*ln(1 - y)*ln(1 + y)*Li(2,y)*a2*f1^2 - 4*ln(1 - y)* ln(1 + y)*Li(2,y)*a2^2*f1 + 8*ln(1 - y)*ln(1 + y)*Li(2,y)*c1*d1*f1 - 4*ln(1 - y)*ln(1 + y)*Li(2,y)*c1*f1^2 - 4*ln(1 - y)*ln(1 + y)*Li(2,y) *c1^2*f1 + 4*ln(1 - y)*ln(1 + y)*Li(2,y)*d1*f1^2 - 4*ln(1 - y)*ln(1 + y)*Li(2,y)*d1^2*f1 - 6*ln(1 - y)*zt(3)*a1*a2*f1 - 12*ln(1 - y)*zt( 3)*a1*c1*d1*[1-y]^-1 + 6*ln(1 - y)*zt(3)*a1*c1*d1 + 6*ln(1 - y)*zt(3) *a1*c1*f1 + 12*ln(1 - y)*zt(3)*a1*c1^2*[1-y]^-1 - 6*ln(1 - y)*zt(3)* a1*c1^2 - 6*ln(1 - y)*zt(3)*a1*d1*f1 + 7*ln(1 - y)*zt(3)*a1*f1^2 - 3* ln(1 - y)*zt(3)*a1^2*f1 - 12*ln(1 - y)*zt(3)*a2*c1*d1*[1-y]^-1 + 6* ln(1 - y)*zt(3)*a2*c1*d1 + 6*ln(1 - y)*zt(3)*a2*c1*f1 + 12*ln(1 - y)* zt(3)*a2*c1^2*[1-y]^-1 - 6*ln(1 - y)*zt(3)*a2*c1^2 - 6*ln(1 - y)*zt(3 )*a2*d1*f1 + 7*ln(1 - y)*zt(3)*a2*f1^2 - 3*ln(1 - y)*zt(3)*a2^2*f1 + 28*ln(1 - y)*zt(3)*c1*d1*f1*[1-y]^-1 - 8*ln(1 - y)*zt(3)*c1*d1*f1 - 12*ln(1 - y)*zt(3)*c1*d1^2*[1-y]^-1 + 6*ln(1 - y)*zt(3)*c1*d1^2 - 7* ln(1 - y)*zt(3)*c1*f1^2 + 16*ln(1 - y)*zt(3)*c1^2*d1*[1-y]^-1 - 8*ln( 1 - y)*zt(3)*c1^2*d1 - 28*ln(1 - y)*zt(3)*c1^2*f1*[1-y]^-1 + 11*ln(1 - y)*zt(3)*c1^2*f1 - 4*ln(1 - y)*zt(3)*c1^3*[1-y]^-1 + 2*ln(1 - y)* zt(3)*c1^3 + 7*ln(1 - y)*zt(3)*d1*f1^2 - 3*ln(1 - y)*zt(3)*d1^2*f1 + 8*ln(1 - y)*Li(3, - y)*a1*a2*f1 - 8*ln(1 - y)*Li(3, - y)*a1*c1*d1* [1-y]^-1 + 4*ln(1 - y)*Li(3, - y)*a1*c1*d1 - 8*ln(1 - y)*Li(3, - y)* a1*c1*f1 + 8*ln(1 - y)*Li(3, - y)*a1*c1^2*[1-y]^-1 - 4*ln(1 - y)*Li(3 , - y)*a1*c1^2 + 8*ln(1 - y)*Li(3, - y)*a1*d1*f1 - 4*ln(1 - y)*Li(3, - y)*a1*f1^2 + 4*ln(1 - y)*Li(3, - y)*a1^2*f1 - 8*ln(1 - y)*Li(3, - y)*a2*c1*d1*[1-y]^-1 + 4*ln(1 - y)*Li(3, - y)*a2*c1*d1 - 8*ln(1 - y)* Li(3, - y)*a2*c1*f1 + 8*ln(1 - y)*Li(3, - y)*a2*c1^2*[1-y]^-1 - 4*ln( 1 - y)*Li(3, - y)*a2*c1^2 + 8*ln(1 - y)*Li(3, - y)*a2*d1*f1 - 4*ln(1 - y)*Li(3, - y)*a2*f1^2 + 4*ln(1 - y)*Li(3, - y)*a2^2*f1 + 8*ln(1 - y)*Li(3, - y)*c1*d1*f1*[1-y]^-1 - 12*ln(1 - y)*Li(3, - y)*c1*d1*f1 - 8*ln(1 - y)*Li(3, - y)*c1*d1^2*[1-y]^-1 + 4*ln(1 - y)*Li(3, - y)*c1* d1^2 + 4*ln(1 - y)*Li(3, - y)*c1*f1^2 + 16*ln(1 - y)*Li(3, - y)*c1^2* d1*[1-y]^-1 - 8*ln(1 - y)*Li(3, - y)*c1^2*d1 - 8*ln(1 - y)*Li(3, - y) *c1^2*f1*[1-y]^-1 + 8*ln(1 - y)*Li(3, - y)*c1^2*f1 - 8*ln(1 - y)*Li(3 , - y)*c1^3*[1-y]^-1 + 4*ln(1 - y)*Li(3, - y)*c1^3 - 4*ln(1 - y)*Li(3 , - y)*d1*f1^2 + 4*ln(1 - y)*Li(3, - y)*d1^2*f1 + 4*ln(1 - y)*Li(3,y) *a1*f1^2 + 4*ln(1 - y)*Li(3,y)*a2*f1^2 - 16*ln(1 - y)*Li(3,y)*c1*d1* f1*[1-y]^-1 + 8*ln(1 - y)*Li(3,y)*c1*d1*f1 - 4*ln(1 - y)*Li(3,y)*c1* f1^2 + 8*ln(1 - y)*Li(3,y)*c1^2*d1*[1-y]^-1 - 4*ln(1 - y)*Li(3,y)* c1^2*d1 + 16*ln(1 - y)*Li(3,y)*c1^2*f1*[1-y]^-1 - 8*ln(1 - y)*Li(3,y) *c1^2*f1 - 8*ln(1 - y)*Li(3,y)*c1^3*[1-y]^-1 + 4*ln(1 - y)*Li(3,y)* c1^3 + 4*ln(1 - y)*Li(3,y)*d1*f1^2 + 12*ln(1 + y)*ln(1 + y)*zt(2)*a1* a2*c1 - 12*ln(1 + y)*ln(1 + y)*zt(2)*a1*a2*d1 + 14*ln(1 + y)*ln(1 + y )*zt(2)*a1*a2*f1 - 6*ln(1 + y)*ln(1 + y)*zt(2)*a1*a2^2 + 12*ln(1 + y) *ln(1 + y)*zt(2)*a1*c1*d1 - 14*ln(1 + y)*ln(1 + y)*zt(2)*a1*c1*f1 - 6 *ln(1 + y)*ln(1 + y)*zt(2)*a1*c1^2 + 14*ln(1 + y)*ln(1 + y)*zt(2)*a1* d1*f1 - 6*ln(1 + y)*ln(1 + y)*zt(2)*a1*d1^2 - 5*ln(1 + y)*ln(1 + y)* zt(2)*a1*f1^2 - 6*ln(1 + y)*ln(1 + y)*zt(2)*a1^2*a2 + 6*ln(1 + y)*ln( 1 + y)*zt(2)*a1^2*c1 - 6*ln(1 + y)*ln(1 + y)*zt(2)*a1^2*d1 + 7*ln(1 + y)*ln(1 + y)*zt(2)*a1^2*f1 - 2*ln(1 + y)*ln(1 + y)*zt(2)*a1^3 + 12 *ln(1 + y)*ln(1 + y)*zt(2)*a2*c1*d1 - 14*ln(1 + y)*ln(1 + y)*zt(2)*a2 *c1*f1 - 6*ln(1 + y)*ln(1 + y)*zt(2)*a2*c1^2 + 14*ln(1 + y)*ln(1 + y) *zt(2)*a2*d1*f1 - 6*ln(1 + y)*ln(1 + y)*zt(2)*a2*d1^2 - 5*ln(1 + y)* ln(1 + y)*zt(2)*a2*f1^2 + 6*ln(1 + y)*ln(1 + y)*zt(2)*a2^2*c1 - 6*ln( 1 + y)*ln(1 + y)*zt(2)*a2^2*d1 + 7*ln(1 + y)*ln(1 + y)*zt(2)*a2^2*f1 - 2*ln(1 + y)*ln(1 + y)*zt(2)*a2^3 - 14*ln(1 + y)*ln(1 + y)*zt(2)*c1 *d1*f1 + 6*ln(1 + y)*ln(1 + y)*zt(2)*c1*d1^2 + 5*ln(1 + y)*ln(1 + y)* zt(2)*c1*f1^2 - 6*ln(1 + y)*ln(1 + y)*zt(2)*c1^2*d1 + 7*ln(1 + y)*ln( 1 + y)*zt(2)*c1^2*f1 + 2*ln(1 + y)*ln(1 + y)*zt(2)*c1^3 - 5*ln(1 + y) *ln(1 + y)*zt(2)*d1*f1^2 + 7*ln(1 + y)*ln(1 + y)*zt(2)*d1^2*f1 - 2* ln(1 + y)*ln(1 + y)*zt(2)*d1^3 + 24*ln(1 + y)*ln(1 + y)*Li(2, - y)*a1 *a2*c1 - 24*ln(1 + y)*ln(1 + y)*Li(2, - y)*a1*a2*d1 + 20*ln(1 + y)* ln(1 + y)*Li(2, - y)*a1*a2*f1 - 12*ln(1 + y)*ln(1 + y)*Li(2, - y)*a1* a2^2 + 24*ln(1 + y)*ln(1 + y)*Li(2, - y)*a1*c1*d1 - 20*ln(1 + y)*ln(1 + y)*Li(2, - y)*a1*c1*f1 - 12*ln(1 + y)*ln(1 + y)*Li(2, - y)*a1*c1^2 + 20*ln(1 + y)*ln(1 + y)*Li(2, - y)*a1*d1*f1 - 12*ln(1 + y)*ln(1 + y )*Li(2, - y)*a1*d1^2 - 6*ln(1 + y)*ln(1 + y)*Li(2, - y)*a1*f1^2 - 12* ln(1 + y)*ln(1 + y)*Li(2, - y)*a1^2*a2 + 12*ln(1 + y)*ln(1 + y)*Li(2, - y)*a1^2*c1 - 12*ln(1 + y)*ln(1 + y)*Li(2, - y)*a1^2*d1 + 10*ln(1 + y)*ln(1 + y)*Li(2, - y)*a1^2*f1 - 4*ln(1 + y)*ln(1 + y)*Li(2, - y) *a1^3 + 24*ln(1 + y)*ln(1 + y)*Li(2, - y)*a2*c1*d1 - 20*ln(1 + y)*ln( 1 + y)*Li(2, - y)*a2*c1*f1 - 12*ln(1 + y)*ln(1 + y)*Li(2, - y)*a2* c1^2 + 20*ln(1 + y)*ln(1 + y)*Li(2, - y)*a2*d1*f1 - 12*ln(1 + y)*ln(1 + y)*Li(2, - y)*a2*d1^2 - 6*ln(1 + y)*ln(1 + y)*Li(2, - y)*a2*f1^2 + 12*ln(1 + y)*ln(1 + y)*Li(2, - y)*a2^2*c1 - 12*ln(1 + y)*ln(1 + y) *Li(2, - y)*a2^2*d1 + 10*ln(1 + y)*ln(1 + y)*Li(2, - y)*a2^2*f1 - 4* ln(1 + y)*ln(1 + y)*Li(2, - y)*a2^3 - 20*ln(1 + y)*ln(1 + y)*Li(2, - y)*c1*d1*f1 + 12*ln(1 + y)*ln(1 + y)*Li(2, - y)*c1*d1^2 + 6*ln(1 + y) *ln(1 + y)*Li(2, - y)*c1*f1^2 - 12*ln(1 + y)*ln(1 + y)*Li(2, - y)* c1^2*d1 + 10*ln(1 + y)*ln(1 + y)*Li(2, - y)*c1^2*f1 + 4*ln(1 + y)*ln( 1 + y)*Li(2, - y)*c1^3 - 6*ln(1 + y)*ln(1 + y)*Li(2, - y)*d1*f1^2 + 10*ln(1 + y)*ln(1 + y)*Li(2, - y)*d1^2*f1 - 4*ln(1 + y)*ln(1 + y)*Li( 2, - y)*d1^3 - 4*ln(1 + y)*ln(1 + y)*Li(2,y)*a1*a2*f1 + 4*ln(1 + y)* ln(1 + y)*Li(2,y)*a1*c1*f1 - 4*ln(1 + y)*ln(1 + y)*Li(2,y)*a1*d1*f1 + 2*ln(1 + y)*ln(1 + y)*Li(2,y)*a1*f1^2 - 2*ln(1 + y)*ln(1 + y)*Li(2 ,y)*a1^2*f1 + 4*ln(1 + y)*ln(1 + y)*Li(2,y)*a2*c1*f1 - 4*ln(1 + y)* ln(1 + y)*Li(2,y)*a2*d1*f1 + 2*ln(1 + y)*ln(1 + y)*Li(2,y)*a2*f1^2 - 2*ln(1 + y)*ln(1 + y)*Li(2,y)*a2^2*f1 + 4*ln(1 + y)*ln(1 + y)*Li(2,y) *c1*d1*f1 - 2*ln(1 + y)*ln(1 + y)*Li(2,y)*c1*f1^2 - 2*ln(1 + y)*ln(1 + y)*Li(2,y)*c1^2*f1 + 2*ln(1 + y)*ln(1 + y)*Li(2,y)*d1*f1^2 - 2*ln( 1 + y)*ln(1 + y)*Li(2,y)*d1^2*f1 - 24*ln(1 + y)*zt(3)*a1*a2*c1 + 24* ln(1 + y)*zt(3)*a1*a2*d1 - 22*ln(1 + y)*zt(3)*a1*a2*f1 + 12*ln(1 + y) *zt(3)*a1*a2^2 - 24*ln(1 + y)*zt(3)*a1*c1*d1 + 22*ln(1 + y)*zt(3)*a1* c1*f1 + 12*ln(1 + y)*zt(3)*a1*c1^2 - 22*ln(1 + y)*zt(3)*a1*d1*f1 + 12 *ln(1 + y)*zt(3)*a1*d1^2 + 7*ln(1 + y)*zt(3)*a1*f1^2 + 12*ln(1 + y)* zt(3)*a1^2*a2 - 12*ln(1 + y)*zt(3)*a1^2*c1 + 12*ln(1 + y)*zt(3)*a1^2* d1 - 11*ln(1 + y)*zt(3)*a1^2*f1 + 4*ln(1 + y)*zt(3)*a1^3 - 24*ln(1 + y)*zt(3)*a2*c1*d1 + 22*ln(1 + y)*zt(3)*a2*c1*f1 + 12*ln(1 + y)*zt(3)* a2*c1^2 - 22*ln(1 + y)*zt(3)*a2*d1*f1 + 12*ln(1 + y)*zt(3)*a2*d1^2 + 7*ln(1 + y)*zt(3)*a2*f1^2 - 12*ln(1 + y)*zt(3)*a2^2*c1 + 12*ln(1 + y) *zt(3)*a2^2*d1 - 11*ln(1 + y)*zt(3)*a2^2*f1 + 4*ln(1 + y)*zt(3)*a2^3 + 22*ln(1 + y)*zt(3)*c1*d1*f1 - 12*ln(1 + y)*zt(3)*c1*d1^2 - 7*ln(1 + y)*zt(3)*c1*f1^2 + 12*ln(1 + y)*zt(3)*c1^2*d1 - 11*ln(1 + y)*zt(3) *c1^2*f1 - 4*ln(1 + y)*zt(3)*c1^3 + 7*ln(1 + y)*zt(3)*d1*f1^2 - 11* ln(1 + y)*zt(3)*d1^2*f1 + 4*ln(1 + y)*zt(3)*d1^3 - 24*ln(1 + y)*Li(3, - y)*a1*a2*c1 + 24*ln(1 + y)*Li(3, - y)*a1*a2*d1 - 16*ln(1 + y)*Li(3 , - y)*a1*a2*f1 + 12*ln(1 + y)*Li(3, - y)*a1*a2^2 - 24*ln(1 + y)*Li(3 , - y)*a1*c1*d1 + 16*ln(1 + y)*Li(3, - y)*a1*c1*f1 + 12*ln(1 + y)*Li( 3, - y)*a1*c1^2 - 16*ln(1 + y)*Li(3, - y)*a1*d1*f1 + 12*ln(1 + y)*Li( 3, - y)*a1*d1^2 + 4*ln(1 + y)*Li(3, - y)*a1*f1^2 + 12*ln(1 + y)*Li(3, - y)*a1^2*a2 - 12*ln(1 + y)*Li(3, - y)*a1^2*c1 + 12*ln(1 + y)*Li(3, - y)*a1^2*d1 - 8*ln(1 + y)*Li(3, - y)*a1^2*f1 + 4*ln(1 + y)*Li(3, - y)*a1^3 - 24*ln(1 + y)*Li(3, - y)*a2*c1*d1 + 16*ln(1 + y)*Li(3, - y)* a2*c1*f1 + 12*ln(1 + y)*Li(3, - y)*a2*c1^2 - 16*ln(1 + y)*Li(3, - y)* a2*d1*f1 + 12*ln(1 + y)*Li(3, - y)*a2*d1^2 + 4*ln(1 + y)*Li(3, - y)* a2*f1^2 - 12*ln(1 + y)*Li(3, - y)*a2^2*c1 + 12*ln(1 + y)*Li(3, - y)* a2^2*d1 - 8*ln(1 + y)*Li(3, - y)*a2^2*f1 + 4*ln(1 + y)*Li(3, - y)* a2^3 + 16*ln(1 + y)*Li(3, - y)*c1*d1*f1 - 12*ln(1 + y)*Li(3, - y)*c1* d1^2 - 4*ln(1 + y)*Li(3, - y)*c1*f1^2 + 12*ln(1 + y)*Li(3, - y)*c1^2* d1 - 8*ln(1 + y)*Li(3, - y)*c1^2*f1 - 4*ln(1 + y)*Li(3, - y)*c1^3 + 4 *ln(1 + y)*Li(3, - y)*d1*f1^2 - 8*ln(1 + y)*Li(3, - y)*d1^2*f1 + 4* ln(1 + y)*Li(3, - y)*d1^3 + 8*ln(1 + y)*Li(3,y)*a1*a2*f1 - 8*ln(1 + y )*Li(3,y)*a1*c1*f1 + 8*ln(1 + y)*Li(3,y)*a1*d1*f1 - 4*ln(1 + y)*Li(3, y)*a1*f1^2 + 4*ln(1 + y)*Li(3,y)*a1^2*f1 - 8*ln(1 + y)*Li(3,y)*a2*c1* f1 + 8*ln(1 + y)*Li(3,y)*a2*d1*f1 - 4*ln(1 + y)*Li(3,y)*a2*f1^2 + 4* ln(1 + y)*Li(3,y)*a2^2*f1 - 8*ln(1 + y)*Li(3,y)*c1*d1*f1 + 4*ln(1 + y )*Li(3,y)*c1*f1^2 + 4*ln(1 + y)*Li(3,y)*c1^2*f1 - 4*ln(1 + y)*Li(3,y) *d1*f1^2 + 4*ln(1 + y)*Li(3,y)*d1^2*f1 + 4/3*ln(y)*ln(1 - y)*ln(1 - y )*ln(1 - y)*a1*a2*f1 - 4/3*ln(y)*ln(1 - y)*ln(1 - y)*ln(1 - y)*a1*c1* f1 + 4/3*ln(y)*ln(1 - y)*ln(1 - y)*ln(1 - y)*a1*d1*f1 - 2*ln(y)*ln(1 - y)*ln(1 - y)*ln(1 - y)*a1*f1^2 + 2/3*ln(y)*ln(1 - y)*ln(1 - y)*ln( 1 - y)*a1^2*f1 - 4/3*ln(y)*ln(1 - y)*ln(1 - y)*ln(1 - y)*a2*c1*f1 + 4/ 3*ln(y)*ln(1 - y)*ln(1 - y)*ln(1 - y)*a2*d1*f1 - 2*ln(y)*ln(1 - y)* ln(1 - y)*ln(1 - y)*a2*f1^2 + 2/3*ln(y)*ln(1 - y)*ln(1 - y)*ln(1 - y) *a2^2*f1 - 4/3*ln(y)*ln(1 - y)*ln(1 - y)*ln(1 - y)*c1*d1*f1 + 2*ln(y) *ln(1 - y)*ln(1 - y)*ln(1 - y)*c1*f1^2 + 2/3*ln(y)*ln(1 - y)*ln(1 - y )*ln(1 - y)*c1^2*f1 - 2*ln(y)*ln(1 - y)*ln(1 - y)*ln(1 - y)*d1*f1^2 + 2/3*ln(y)*ln(1 - y)*ln(1 - y)*ln(1 - y)*d1^2*f1 - 4*ln(y)*ln(1 - y )*ln(1 - y)*ln(1 + y)*a1*a2*f1 + 4*ln(y)*ln(1 - y)*ln(1 - y)*ln(1 + y )*a1*c1*f1 - 4*ln(y)*ln(1 - y)*ln(1 - y)*ln(1 + y)*a1*d1*f1 + 2*ln(y) *ln(1 - y)*ln(1 - y)*ln(1 + y)*a1*f1^2 - 2*ln(y)*ln(1 - y)*ln(1 - y)* ln(1 + y)*a1^2*f1 + 4*ln(y)*ln(1 - y)*ln(1 - y)*ln(1 + y)*a2*c1*f1 - 4*ln(y)*ln(1 - y)*ln(1 - y)*ln(1 + y)*a2*d1*f1 + 2*ln(y)*ln(1 - y)* ln(1 - y)*ln(1 + y)*a2*f1^2 - 2*ln(y)*ln(1 - y)*ln(1 - y)*ln(1 + y)* a2^2*f1 + 4*ln(y)*ln(1 - y)*ln(1 - y)*ln(1 + y)*c1*d1*f1 - 2*ln(y)* ln(1 - y)*ln(1 - y)*ln(1 + y)*c1*f1^2 - 2*ln(y)*ln(1 - y)*ln(1 - y)* ln(1 + y)*c1^2*f1 + 2*ln(y)*ln(1 - y)*ln(1 - y)*ln(1 + y)*d1*f1^2 - 2 *ln(y)*ln(1 - y)*ln(1 - y)*ln(1 + y)*d1^2*f1 - 4*ln(y)*ln(1 - y)*ln(1 + y)*ln(1 + y)*a1*a2*f1 + 4*ln(y)*ln(1 - y)*ln(1 + y)*ln(1 + y)*a1* c1*f1 - 4*ln(y)*ln(1 - y)*ln(1 + y)*ln(1 + y)*a1*d1*f1 + 2*ln(y)*ln(1 - y)*ln(1 + y)*ln(1 + y)*a1*f1^2 - 2*ln(y)*ln(1 - y)*ln(1 + y)*ln(1 + y)*a1^2*f1 + 4*ln(y)*ln(1 - y)*ln(1 + y)*ln(1 + y)*a2*c1*f1 - 4* ln(y)*ln(1 - y)*ln(1 + y)*ln(1 + y)*a2*d1*f1 + 2*ln(y)*ln(1 - y)*ln(1 + y)*ln(1 + y)*a2*f1^2 - 2*ln(y)*ln(1 - y)*ln(1 + y)*ln(1 + y)*a2^2* f1 + 4*ln(y)*ln(1 - y)*ln(1 + y)*ln(1 + y)*c1*d1*f1 - 2*ln(y)*ln(1 - y)*ln(1 + y)*ln(1 + y)*c1*f1^2 - 2*ln(y)*ln(1 - y)*ln(1 + y)*ln(1 + y )*c1^2*f1 + 2*ln(y)*ln(1 - y)*ln(1 + y)*ln(1 + y)*d1*f1^2 - 2*ln(y)* ln(1 - y)*ln(1 + y)*ln(1 + y)*d1^2*f1 + 4*ln(y)*ln(1 - y)*zt(2)*a1*a2 *f1 - 4*ln(y)*ln(1 - y)*zt(2)*a1*c1*d1*[1-y]^-1 + 2*ln(y)*ln(1 - y)* zt(2)*a1*c1*d1 - 4*ln(y)*ln(1 - y)*zt(2)*a1*c1*f1 + 4*ln(y)*ln(1 - y) *zt(2)*a1*c1^2*[1-y]^-1 - 2*ln(y)*ln(1 - y)*zt(2)*a1*c1^2 + 4*ln(y)* ln(1 - y)*zt(2)*a1*d1*f1 - 6*ln(y)*ln(1 - y)*zt(2)*a1*f1^2 + 2*ln(y)* ln(1 - y)*zt(2)*a1^2*f1 - 4*ln(y)*ln(1 - y)*zt(2)*a2*c1*d1*[1-y]^-1 + 2*ln(y)*ln(1 - y)*zt(2)*a2*c1*d1 - 4*ln(y)*ln(1 - y)*zt(2)*a2*c1* f1 + 4*ln(y)*ln(1 - y)*zt(2)*a2*c1^2*[1-y]^-1 - 2*ln(y)*ln(1 - y)*zt( 2)*a2*c1^2 + 4*ln(y)*ln(1 - y)*zt(2)*a2*d1*f1 - 6*ln(y)*ln(1 - y)*zt( 2)*a2*f1^2 + 2*ln(y)*ln(1 - y)*zt(2)*a2^2*f1 + 12*ln(y)*ln(1 - y)*zt( 2)*c1*d1*f1*[1-y]^-1 - 10*ln(y)*ln(1 - y)*zt(2)*c1*d1*f1 - 4*ln(y)* ln(1 - y)*zt(2)*c1*d1^2*[1-y]^-1 + 2*ln(y)*ln(1 - y)*zt(2)*c1*d1^2 + 6*ln(y)*ln(1 - y)*zt(2)*c1*f1^2 + 8*ln(y)*ln(1 - y)*zt(2)*c1^2*d1* [1-y]^-1 - 4*ln(y)*ln(1 - y)*zt(2)*c1^2*d1 - 12*ln(y)*ln(1 - y)*zt(2) *c1^2*f1*[1-y]^-1 + 8*ln(y)*ln(1 - y)*zt(2)*c1^2*f1 - 4*ln(y)*ln(1 - y)*zt(2)*c1^3*[1-y]^-1 + 2*ln(y)*ln(1 - y)*zt(2)*c1^3 - 6*ln(y)*ln(1 - y)*zt(2)*d1*f1^2 + 2*ln(y)*ln(1 - y)*zt(2)*d1^2*f1 + 8*ln(y)*ln(1 - y)*Li(2, - y)*a1*c1*d1*[1-y]^-1 - 4*ln(y)*ln(1 - y)*Li(2, - y)*a1* c1*d1 - 8*ln(y)*ln(1 - y)*Li(2, - y)*a1*c1^2*[1-y]^-1 + 4*ln(y)*ln(1 - y)*Li(2, - y)*a1*c1^2 + 8*ln(y)*ln(1 - y)*Li(2, - y)*a2*c1*d1* [1-y]^-1 - 4*ln(y)*ln(1 - y)*Li(2, - y)*a2*c1*d1 - 8*ln(y)*ln(1 - y)* Li(2, - y)*a2*c1^2*[1-y]^-1 + 4*ln(y)*ln(1 - y)*Li(2, - y)*a2*c1^2 - 8*ln(y)*ln(1 - y)*Li(2, - y)*c1*d1*f1*[1-y]^-1 + 4*ln(y)*ln(1 - y)* Li(2, - y)*c1*d1*f1 + 8*ln(y)*ln(1 - y)*Li(2, - y)*c1*d1^2*[1-y]^-1 - 4*ln(y)*ln(1 - y)*Li(2, - y)*c1*d1^2 - 16*ln(y)*ln(1 - y)*Li(2, - y)*c1^2*d1*[1-y]^-1 + 8*ln(y)*ln(1 - y)*Li(2, - y)*c1^2*d1 + 8*ln(y)* ln(1 - y)*Li(2, - y)*c1^2*f1*[1-y]^-1 - 4*ln(y)*ln(1 - y)*Li(2, - y)* c1^2*f1 + 8*ln(y)*ln(1 - y)*Li(2, - y)*c1^3*[1-y]^-1 - 4*ln(y)*ln(1 - y)*Li(2, - y)*c1^3 + 8*ln(y)*ln(1 - y)*Li(2,y)*c1*d1*f1*[1-y]^-1 - 4*ln(y)*ln(1 - y)*Li(2,y)*c1*d1*f1 - 8*ln(y)*ln(1 - y)*Li(2,y)* c1^2*d1*[1-y]^-1 + 4*ln(y)*ln(1 - y)*Li(2,y)*c1^2*d1 - 8*ln(y)*ln(1 - y)*Li(2,y)*c1^2*f1*[1-y]^-1 + 4*ln(y)*ln(1 - y)*Li(2,y)*c1^2*f1 + 8*ln(y)*ln(1 - y)*Li(2,y)*c1^3*[1-y]^-1 - 4*ln(y)*ln(1 - y)*Li(2,y)* c1^3 + 8*ln(y)*ln(1 + y)*ln(1 + y)*ln(1 + y)*a1*a2*c1 - 8*ln(y)*ln(1 + y)*ln(1 + y)*ln(1 + y)*a1*a2*d1 + 20/3*ln(y)*ln(1 + y)*ln(1 + y)* ln(1 + y)*a1*a2*f1 - 4*ln(y)*ln(1 + y)*ln(1 + y)*ln(1 + y)*a1*a2^2 + 8*ln(y)*ln(1 + y)*ln(1 + y)*ln(1 + y)*a1*c1*d1 - 20/3*ln(y)*ln(1 + y) *ln(1 + y)*ln(1 + y)*a1*c1*f1 - 4*ln(y)*ln(1 + y)*ln(1 + y)*ln(1 + y) *a1*c1^2 + 20/3*ln(y)*ln(1 + y)*ln(1 + y)*ln(1 + y)*a1*d1*f1 - 4*ln(y )*ln(1 + y)*ln(1 + y)*ln(1 + y)*a1*d1^2 - 2*ln(y)*ln(1 + y)*ln(1 + y) *ln(1 + y)*a1*f1^2 - 4*ln(y)*ln(1 + y)*ln(1 + y)*ln(1 + y)*a1^2*a2 + 4*ln(y)*ln(1 + y)*ln(1 + y)*ln(1 + y)*a1^2*c1 - 4*ln(y)*ln(1 + y)*ln( 1 + y)*ln(1 + y)*a1^2*d1 + 10/3*ln(y)*ln(1 + y)*ln(1 + y)*ln(1 + y)* a1^2*f1 - 4/3*ln(y)*ln(1 + y)*ln(1 + y)*ln(1 + y)*a1^3 + 8*ln(y)*ln(1 + y)*ln(1 + y)*ln(1 + y)*a2*c1*d1 - 20/3*ln(y)*ln(1 + y)*ln(1 + y)* ln(1 + y)*a2*c1*f1 - 4*ln(y)*ln(1 + y)*ln(1 + y)*ln(1 + y)*a2*c1^2 + 20/3*ln(y)*ln(1 + y)*ln(1 + y)*ln(1 + y)*a2*d1*f1 - 4*ln(y)*ln(1 + y) *ln(1 + y)*ln(1 + y)*a2*d1^2 - 2*ln(y)*ln(1 + y)*ln(1 + y)*ln(1 + y)* a2*f1^2 + 4*ln(y)*ln(1 + y)*ln(1 + y)*ln(1 + y)*a2^2*c1 - 4*ln(y)*ln( 1 + y)*ln(1 + y)*ln(1 + y)*a2^2*d1 + 10/3*ln(y)*ln(1 + y)*ln(1 + y)* ln(1 + y)*a2^2*f1 - 4/3*ln(y)*ln(1 + y)*ln(1 + y)*ln(1 + y)*a2^3 - 20/ 3*ln(y)*ln(1 + y)*ln(1 + y)*ln(1 + y)*c1*d1*f1 + 4*ln(y)*ln(1 + y)* ln(1 + y)*ln(1 + y)*c1*d1^2 + 2*ln(y)*ln(1 + y)*ln(1 + y)*ln(1 + y)* c1*f1^2 - 4*ln(y)*ln(1 + y)*ln(1 + y)*ln(1 + y)*c1^2*d1 + 10/3*ln(y)* ln(1 + y)*ln(1 + y)*ln(1 + y)*c1^2*f1 + 4/3*ln(y)*ln(1 + y)*ln(1 + y) *ln(1 + y)*c1^3 - 2*ln(y)*ln(1 + y)*ln(1 + y)*ln(1 + y)*d1*f1^2 + 10/ 3*ln(y)*ln(1 + y)*ln(1 + y)*ln(1 + y)*d1^2*f1 - 4/3*ln(y)*ln(1 + y)* ln(1 + y)*ln(1 + y)*d1^3 - 12*ln(y)*ln(1 + y)*zt(2)*a1*a2*c1 + 12*ln( y)*ln(1 + y)*zt(2)*a1*a2*d1 - 16*ln(y)*ln(1 + y)*zt(2)*a1*a2*f1 + 6* ln(y)*ln(1 + y)*zt(2)*a1*a2^2 - 12*ln(y)*ln(1 + y)*zt(2)*a1*c1*d1 + 16*ln(y)*ln(1 + y)*zt(2)*a1*c1*f1 + 6*ln(y)*ln(1 + y)*zt(2)*a1*c1^2 - 16*ln(y)*ln(1 + y)*zt(2)*a1*d1*f1 + 6*ln(y)*ln(1 + y)*zt(2)*a1* d1^2 + 6*ln(y)*ln(1 + y)*zt(2)*a1*f1^2 + 6*ln(y)*ln(1 + y)*zt(2)*a1^2 *a2 - 6*ln(y)*ln(1 + y)*zt(2)*a1^2*c1 + 6*ln(y)*ln(1 + y)*zt(2)*a1^2* d1 - 8*ln(y)*ln(1 + y)*zt(2)*a1^2*f1 + 2*ln(y)*ln(1 + y)*zt(2)*a1^3 - 12*ln(y)*ln(1 + y)*zt(2)*a2*c1*d1 + 16*ln(y)*ln(1 + y)*zt(2)*a2*c1 *f1 + 6*ln(y)*ln(1 + y)*zt(2)*a2*c1^2 - 16*ln(y)*ln(1 + y)*zt(2)*a2* d1*f1 + 6*ln(y)*ln(1 + y)*zt(2)*a2*d1^2 + 6*ln(y)*ln(1 + y)*zt(2)*a2* f1^2 - 6*ln(y)*ln(1 + y)*zt(2)*a2^2*c1 + 6*ln(y)*ln(1 + y)*zt(2)*a2^2 *d1 - 8*ln(y)*ln(1 + y)*zt(2)*a2^2*f1 + 2*ln(y)*ln(1 + y)*zt(2)*a2^3 + 16*ln(y)*ln(1 + y)*zt(2)*c1*d1*f1 - 6*ln(y)*ln(1 + y)*zt(2)*c1* d1^2 - 6*ln(y)*ln(1 + y)*zt(2)*c1*f1^2 + 6*ln(y)*ln(1 + y)*zt(2)*c1^2 *d1 - 8*ln(y)*ln(1 + y)*zt(2)*c1^2*f1 - 2*ln(y)*ln(1 + y)*zt(2)*c1^3 + 6*ln(y)*ln(1 + y)*zt(2)*d1*f1^2 - 8*ln(y)*ln(1 + y)*zt(2)*d1^2*f1 + 2*ln(y)*ln(1 + y)*zt(2)*d1^3 + ln(y)*ln(y)*ln(1 - y)*ln(1 - y)*a1* f1^2 + ln(y)*ln(y)*ln(1 - y)*ln(1 - y)*a2*f1^2 - ln(y)*ln(y)*ln(1 - y )*ln(1 - y)*c1*f1^2 - 2*ln(y)*ln(y)*ln(1 - y)*ln(1 - y)*c1^2*d1* [1-y]^-1 + ln(y)*ln(y)*ln(1 - y)*ln(1 - y)*c1^2*d1 + 2*ln(y)*ln(y)* ln(1 - y)*ln(1 - y)*c1^3*[1-y]^-1 - ln(y)*ln(y)*ln(1 - y)*ln(1 - y)* c1^3 + ln(y)*ln(y)*ln(1 - y)*ln(1 - y)*d1*f1^2 + 4*ln(y)*ln(y)*ln(1 - y)*ln(1 + y)*a1*a2*f1 - 4*ln(y)*ln(y)*ln(1 - y)*ln(1 + y)*a1*c1*f1 + 4*ln(y)*ln(y)*ln(1 - y)*ln(1 + y)*a1*d1*f1 - 2*ln(y)*ln(y)*ln(1 - y)*ln(1 + y)*a1*f1^2 + 2*ln(y)*ln(y)*ln(1 - y)*ln(1 + y)*a1^2*f1 - 4* ln(y)*ln(y)*ln(1 - y)*ln(1 + y)*a2*c1*f1 + 4*ln(y)*ln(y)*ln(1 - y)* ln(1 + y)*a2*d1*f1 - 2*ln(y)*ln(y)*ln(1 - y)*ln(1 + y)*a2*f1^2 + 2* ln(y)*ln(y)*ln(1 - y)*ln(1 + y)*a2^2*f1 - 4*ln(y)*ln(y)*ln(1 - y)*ln( 1 + y)*c1*d1*f1 + 2*ln(y)*ln(y)*ln(1 - y)*ln(1 + y)*c1*f1^2 + 2*ln(y) *ln(y)*ln(1 - y)*ln(1 + y)*c1^2*f1 - 2*ln(y)*ln(y)*ln(1 - y)*ln(1 + y )*d1*f1^2 + 2*ln(y)*ln(y)*ln(1 - y)*ln(1 + y)*d1^2*f1 - 6*ln(y)*ln(y) *ln(1 + y)*ln(1 + y)*a1*a2*c1 + 6*ln(y)*ln(y)*ln(1 + y)*ln(1 + y)*a1* a2*d1 - 4*ln(y)*ln(y)*ln(1 + y)*ln(1 + y)*a1*a2*f1 + 3*ln(y)*ln(y)* ln(1 + y)*ln(1 + y)*a1*a2^2 - 6*ln(y)*ln(y)*ln(1 + y)*ln(1 + y)*a1*c1 *d1 + 4*ln(y)*ln(y)*ln(1 + y)*ln(1 + y)*a1*c1*f1 + 3*ln(y)*ln(y)*ln(1 + y)*ln(1 + y)*a1*c1^2 - 4*ln(y)*ln(y)*ln(1 + y)*ln(1 + y)*a1*d1*f1 + 3*ln(y)*ln(y)*ln(1 + y)*ln(1 + y)*a1*d1^2 + ln(y)*ln(y)*ln(1 + y)* ln(1 + y)*a1*f1^2 + 3*ln(y)*ln(y)*ln(1 + y)*ln(1 + y)*a1^2*a2 - 3*ln( y)*ln(y)*ln(1 + y)*ln(1 + y)*a1^2*c1 + 3*ln(y)*ln(y)*ln(1 + y)*ln(1 + y)*a1^2*d1 - 2*ln(y)*ln(y)*ln(1 + y)*ln(1 + y)*a1^2*f1 + ln(y)*ln( y)*ln(1 + y)*ln(1 + y)*a1^3 - 6*ln(y)*ln(y)*ln(1 + y)*ln(1 + y)*a2*c1 *d1 + 4*ln(y)*ln(y)*ln(1 + y)*ln(1 + y)*a2*c1*f1 + 3*ln(y)*ln(y)*ln(1 + y)*ln(1 + y)*a2*c1^2 - 4*ln(y)*ln(y)*ln(1 + y)*ln(1 + y)*a2*d1*f1 + 3*ln(y)*ln(y)*ln(1 + y)*ln(1 + y)*a2*d1^2 + ln(y)*ln(y)*ln(1 + y)* ln(1 + y)*a2*f1^2 - 3*ln(y)*ln(y)*ln(1 + y)*ln(1 + y)*a2^2*c1 + 3*ln( y)*ln(y)*ln(1 + y)*ln(1 + y)*a2^2*d1 - 2*ln(y)*ln(y)*ln(1 + y)*ln(1 + y)*a2^2*f1 + ln(y)*ln(y)*ln(1 + y)*ln(1 + y)*a2^3 + 4*ln(y)*ln(y)* ln(1 + y)*ln(1 + y)*c1*d1*f1 - 3*ln(y)*ln(y)*ln(1 + y)*ln(1 + y)*c1* d1^2 - ln(y)*ln(y)*ln(1 + y)*ln(1 + y)*c1*f1^2 + 3*ln(y)*ln(y)*ln(1 + y)*ln(1 + y)*c1^2*d1 - 2*ln(y)*ln(y)*ln(1 + y)*ln(1 + y)*c1^2*f1 - ln(y)*ln(y)*ln(1 + y)*ln(1 + y)*c1^3 + ln(y)*ln(y)*ln(1 + y)*ln(1 + y)*d1*f1^2 - 2*ln(y)*ln(y)*ln(1 + y)*ln(1 + y)*d1^2*f1 + ln(y)*ln( y)*ln(1 + y)*ln(1 + y)*d1^3 - 1/3*ln(y)*ln(y)*ln(y)*ln(1 - y)*a1*a2* f1 + 2/3*ln(y)*ln(y)*ln(y)*ln(1 - y)*a1*c1*d1*[1-y]^-1 - 1/3*ln(y)* ln(y)*ln(y)*ln(1 - y)*a1*c1*d1 + 1/3*ln(y)*ln(y)*ln(y)*ln(1 - y)*a1* c1*f1 - 2/3*ln(y)*ln(y)*ln(y)*ln(1 - y)*a1*c1^2*[1-y]^-1 + 1/3*ln(y)* ln(y)*ln(y)*ln(1 - y)*a1*c1^2 - 1/3*ln(y)*ln(y)*ln(y)*ln(1 - y)*a1*d1 *f1 - 1/3*ln(y)*ln(y)*ln(y)*ln(1 - y)*a1^2*f1 + 2/3*ln(y)*ln(y)*ln(y) *ln(1 - y)*a2*c1*d1*[1-y]^-1 - 1/3*ln(y)*ln(y)*ln(y)*ln(1 - y)*a2*c1* d1 + 1/3*ln(y)*ln(y)*ln(y)*ln(1 - y)*a2*c1*f1 - 2/3*ln(y)*ln(y)*ln(y) *ln(1 - y)*a2*c1^2*[1-y]^-1 + 1/3*ln(y)*ln(y)*ln(y)*ln(1 - y)*a2*c1^2 - 1/3*ln(y)*ln(y)*ln(y)*ln(1 - y)*a2*d1*f1 - 1/3*ln(y)*ln(y)*ln(y)* ln(1 - y)*a2^2*f1 + 1/3*ln(y)*ln(y)*ln(y)*ln(1 - y)*c1*d1*f1 + 2/3* ln(y)*ln(y)*ln(y)*ln(1 - y)*c1*d1^2*[1-y]^-1 - 1/3*ln(y)*ln(y)*ln(y)* ln(1 - y)*c1*d1^2 - 2/3*ln(y)*ln(y)*ln(y)*ln(1 - y)*c1^2*d1*[1-y]^-1 + 1/3*ln(y)*ln(y)*ln(y)*ln(1 - y)*c1^2*d1 - 1/3*ln(y)*ln(y)*ln(y)* ln(1 - y)*d1^2*f1 + ln(y)*ln(y)*ln(y)*ln(1 + y)*a1*a2*c1 - ln(y)*ln(y )*ln(y)*ln(1 + y)*a1*a2*d1 + 1/3*ln(y)*ln(y)*ln(y)*ln(1 + y)*a1*a2*f1 - 2/3*ln(y)*ln(y)*ln(y)*ln(1 + y)*a1*a2^2 + ln(y)*ln(y)*ln(y)*ln(1 + y)*a1*c1*d1 - 1/3*ln(y)*ln(y)*ln(y)*ln(1 + y)*a1*c1*f1 - 1/3*ln(y) *ln(y)*ln(y)*ln(1 + y)*a1*c1^2 + 1/3*ln(y)*ln(y)*ln(y)*ln(1 + y)*a1* d1*f1 - 2/3*ln(y)*ln(y)*ln(y)*ln(1 + y)*a1*d1^2 - 2/3*ln(y)*ln(y)*ln( y)*ln(1 + y)*a1^2*a2 + 2/3*ln(y)*ln(y)*ln(y)*ln(1 + y)*a1^2*c1 - 2/3* ln(y)*ln(y)*ln(y)*ln(1 + y)*a1^2*d1 + 1/3*ln(y)*ln(y)*ln(y)*ln(1 + y) *a1^2*f1 - 1/3*ln(y)*ln(y)*ln(y)*ln(1 + y)*a1^3 + ln(y)*ln(y)*ln(y)* ln(1 + y)*a2*c1*d1 - 1/3*ln(y)*ln(y)*ln(y)*ln(1 + y)*a2*c1*f1 - 1/3* ln(y)*ln(y)*ln(y)*ln(1 + y)*a2*c1^2 + 1/3*ln(y)*ln(y)*ln(y)*ln(1 + y) *a2*d1*f1 - 2/3*ln(y)*ln(y)*ln(y)*ln(1 + y)*a2*d1^2 + 2/3*ln(y)*ln(y) *ln(y)*ln(1 + y)*a2^2*c1 - 2/3*ln(y)*ln(y)*ln(y)*ln(1 + y)*a2^2*d1 + 1/3*ln(y)*ln(y)*ln(y)*ln(1 + y)*a2^2*f1 - 1/3*ln(y)*ln(y)*ln(y)*ln(1 + y)*a2^3 - 1/3*ln(y)*ln(y)*ln(y)*ln(1 + y)*c1*d1*f1 + 2/3*ln(y)*ln( y)*ln(y)*ln(1 + y)*c1*d1^2 - 1/3*ln(y)*ln(y)*ln(y)*ln(1 + y)*c1^2*d1 + 1/3*ln(y)*ln(y)*ln(y)*ln(1 + y)*d1^2*f1 - 1/3*ln(y)*ln(y)*ln(y)* ln(1 + y)*d1^3 + 1/12*ln(y)*ln(y)*ln(y)*ln(y)*a1*a2*c1*[1-y]^-1 - 1/ 12*ln(y)*ln(y)*ln(y)*ln(y)*a1*a2*c1 - 1/12*ln(y)*ln(y)*ln(y)*ln(y)*a1 *a2*d1*[1-y]^-1 + 1/12*ln(y)*ln(y)*ln(y)*ln(y)*a1*a2*d1 + 1/24*ln(y)* ln(y)*ln(y)*ln(y)*a1*a2^2 + 1/12*ln(y)*ln(y)*ln(y)*ln(y)*a1*c1*d1* [1-y]^-1 - 1/12*ln(y)*ln(y)*ln(y)*ln(y)*a1*c1*d1 - 1/12*ln(y)*ln(y)* ln(y)*ln(y)*a1*d1^2*[1-y]^-1 + 1/12*ln(y)*ln(y)*ln(y)*ln(y)*a1*d1^2 + 1/24*ln(y)*ln(y)*ln(y)*ln(y)*a1^2*a2 + 1/12*ln(y)*ln(y)*ln(y)*ln(y )*a1^2*c1*[1-y]^-1 - 1/12*ln(y)*ln(y)*ln(y)*ln(y)*a1^2*c1 - 1/12*ln(y )*ln(y)*ln(y)*ln(y)*a1^2*d1*[1-y]^-1 + 1/12*ln(y)*ln(y)*ln(y)*ln(y)* a1^2*d1 + 1/24*ln(y)*ln(y)*ln(y)*ln(y)*a1^3 + 1/12*ln(y)*ln(y)*ln(y)* ln(y)*a2*c1*d1*[1-y]^-1 - 1/12*ln(y)*ln(y)*ln(y)*ln(y)*a2*c1*d1 - 1/ 12*ln(y)*ln(y)*ln(y)*ln(y)*a2*d1^2*[1-y]^-1 + 1/12*ln(y)*ln(y)*ln(y)* ln(y)*a2*d1^2 + 1/12*ln(y)*ln(y)*ln(y)*ln(y)*a2^2*c1*[1-y]^-1 - 1/12* ln(y)*ln(y)*ln(y)*ln(y)*a2^2*c1 - 1/12*ln(y)*ln(y)*ln(y)*ln(y)*a2^2* d1*[1-y]^-1 + 1/12*ln(y)*ln(y)*ln(y)*ln(y)*a2^2*d1 + 1/24*ln(y)*ln(y) *ln(y)*ln(y)*a2^3 + 1/12*ln(y)*ln(y)*ln(y)*ln(y)*c1*d1^2*[1-y]^-1 - 1/ 12*ln(y)*ln(y)*ln(y)*ln(y)*c1*d1^2 - 1/12*ln(y)*ln(y)*ln(y)*ln(y)* d1^3*[1-y]^-1 + 1/12*ln(y)*ln(y)*ln(y)*ln(y)*d1^3 - 2*ln(y)*ln(y)*zt( 2)*a1*a2*c1*[1-y]^-1 + 5/2*ln(y)*ln(y)*zt(2)*a1*a2*c1 + 2*ln(y)*ln(y) *zt(2)*a1*a2*d1*[1-y]^-1 - 5/2*ln(y)*ln(y)*zt(2)*a1*a2*d1 + 3/2*ln(y) *ln(y)*zt(2)*a1*a2*f1 - ln(y)*ln(y)*zt(2)*a1*a2^2 - 3*ln(y)*ln(y)*zt( 2)*a1*c1*d1*[1-y]^-1 + 3*ln(y)*ln(y)*zt(2)*a1*c1*d1 + 3*ln(y)*ln(y)* zt(2)*a1*c1*f1*[1-y]^-1 - 3*ln(y)*ln(y)*zt(2)*a1*c1*f1 + ln(y)*ln(y)* zt(2)*a1*c1^2*[1-y]^-1 - ln(y)*ln(y)*zt(2)*a1*c1^2 - 3*ln(y)*ln(y)* zt(2)*a1*d1*f1*[1-y]^-1 + 3*ln(y)*ln(y)*zt(2)*a1*d1*f1 + 2*ln(y)*ln(y )*zt(2)*a1*d1^2*[1-y]^-1 - 2*ln(y)*ln(y)*zt(2)*a1*d1^2 - ln(y)*ln(y)* zt(2)*a1^2*a2 - ln(y)*ln(y)*zt(2)*a1^2*c1*[1-y]^-1 + 3/2*ln(y)*ln(y)* zt(2)*a1^2*c1 + ln(y)*ln(y)*zt(2)*a1^2*d1*[1-y]^-1 - 3/2*ln(y)*ln(y)* zt(2)*a1^2*d1 + 3/2*ln(y)*ln(y)*zt(2)*a1^2*f1 - 1/2*ln(y)*ln(y)*zt(2) *a1^3 - 3*ln(y)*ln(y)*zt(2)*a2*c1*d1*[1-y]^-1 + 3*ln(y)*ln(y)*zt(2)* a2*c1*d1 + 3*ln(y)*ln(y)*zt(2)*a2*c1*f1*[1-y]^-1 - 3*ln(y)*ln(y)*zt(2 )*a2*c1*f1 + ln(y)*ln(y)*zt(2)*a2*c1^2*[1-y]^-1 - ln(y)*ln(y)*zt(2)* a2*c1^2 - 3*ln(y)*ln(y)*zt(2)*a2*d1*f1*[1-y]^-1 + 3*ln(y)*ln(y)*zt(2) *a2*d1*f1 + 2*ln(y)*ln(y)*zt(2)*a2*d1^2*[1-y]^-1 - 2*ln(y)*ln(y)*zt(2 )*a2*d1^2 - ln(y)*ln(y)*zt(2)*a2^2*c1*[1-y]^-1 + 3/2*ln(y)*ln(y)*zt(2 )*a2^2*c1 + ln(y)*ln(y)*zt(2)*a2^2*d1*[1-y]^-1 - 3/2*ln(y)*ln(y)*zt(2 )*a2^2*d1 + 3/2*ln(y)*ln(y)*zt(2)*a2^2*f1 - 1/2*ln(y)*ln(y)*zt(2)* a2^3 + 3*ln(y)*ln(y)*zt(2)*c1*d1*f1*[1-y]^-1 - 3*ln(y)*ln(y)*zt(2)*c1 *d1*f1 - 2*ln(y)*ln(y)*zt(2)*c1*d1^2*[1-y]^-1 + 2*ln(y)*ln(y)*zt(2)* c1*d1^2 + ln(y)*ln(y)*zt(2)*c1^2*d1*[1-y]^-1 - ln(y)*ln(y)*zt(2)*c1^2 *d1 - 3*ln(y)*ln(y)*zt(2)*d1^2*f1*[1-y]^-1 + 3*ln(y)*ln(y)*zt(2)*d1^2 *f1 + ln(y)*ln(y)*zt(2)*d1^3*[1-y]^-1 - ln(y)*ln(y)*zt(2)*d1^3 + 4* ln(y)*ln(y)*Li(2, - y)*a1*a2*c1*[1-y]^-1 - 5*ln(y)*ln(y)*Li(2, - y)* a1*a2*c1 - 4*ln(y)*ln(y)*Li(2, - y)*a1*a2*d1*[1-y]^-1 + 5*ln(y)*ln(y) *Li(2, - y)*a1*a2*d1 - ln(y)*ln(y)*Li(2, - y)*a1*a2*f1 + ln(y)*ln(y)* Li(2, - y)*a1*a2^2 + 8*ln(y)*ln(y)*Li(2, - y)*a1*c1*d1*[1-y]^-1 - 7* ln(y)*ln(y)*Li(2, - y)*a1*c1*d1 - 2*ln(y)*ln(y)*Li(2, - y)*a1*c1*f1* [1-y]^-1 + 2*ln(y)*ln(y)*Li(2, - y)*a1*c1*f1 - 4*ln(y)*ln(y)*Li(2, - y)*a1*c1^2*[1-y]^-1 + 4*ln(y)*ln(y)*Li(2, - y)*a1*c1^2 + 2*ln(y)*ln(y )*Li(2, - y)*a1*d1*f1*[1-y]^-1 - 2*ln(y)*ln(y)*Li(2, - y)*a1*d1*f1 - 4*ln(y)*ln(y)*Li(2, - y)*a1*d1^2*[1-y]^-1 + 3*ln(y)*ln(y)*Li(2, - y)* a1*d1^2 + ln(y)*ln(y)*Li(2, - y)*a1^2*a2 + 2*ln(y)*ln(y)*Li(2, - y)* a1^2*c1*[1-y]^-1 - 2*ln(y)*ln(y)*Li(2, - y)*a1^2*c1 - 2*ln(y)*ln(y)* Li(2, - y)*a1^2*d1*[1-y]^-1 + 2*ln(y)*ln(y)*Li(2, - y)*a1^2*d1 + 8* ln(y)*ln(y)*Li(2, - y)*a2*c1*d1*[1-y]^-1 - 7*ln(y)*ln(y)*Li(2, - y)* a2*c1*d1 - 2*ln(y)*ln(y)*Li(2, - y)*a2*c1*f1*[1-y]^-1 + 2*ln(y)*ln(y) *Li(2, - y)*a2*c1*f1 - 4*ln(y)*ln(y)*Li(2, - y)*a2*c1^2*[1-y]^-1 + 4* ln(y)*ln(y)*Li(2, - y)*a2*c1^2 + 2*ln(y)*ln(y)*Li(2, - y)*a2*d1*f1* [1-y]^-1 - 2*ln(y)*ln(y)*Li(2, - y)*a2*d1*f1 - 4*ln(y)*ln(y)*Li(2, - y)*a2*d1^2*[1-y]^-1 + 3*ln(y)*ln(y)*Li(2, - y)*a2*d1^2 + 2*ln(y)*ln(y )*Li(2, - y)*a2^2*c1*[1-y]^-1 - 2*ln(y)*ln(y)*Li(2, - y)*a2^2*c1 - 2* ln(y)*ln(y)*Li(2, - y)*a2^2*d1*[1-y]^-1 + 2*ln(y)*ln(y)*Li(2, - y)* a2^2*d1 - 4*ln(y)*ln(y)*Li(2, - y)*c1*d1*f1*[1-y]^-1 + 3*ln(y)*ln(y)* Li(2, - y)*c1*d1*f1 + 6*ln(y)*ln(y)*Li(2, - y)*c1*d1^2*[1-y]^-1 - 4* ln(y)*ln(y)*Li(2, - y)*c1*d1^2 - 6*ln(y)*ln(y)*Li(2, - y)*c1^2*d1* [1-y]^-1 + 5*ln(y)*ln(y)*Li(2, - y)*c1^2*d1 + 2*ln(y)*ln(y)*Li(2, - y )*c1^2*f1*[1-y]^-1 - 2*ln(y)*ln(y)*Li(2, - y)*c1^2*f1 + 2*ln(y)*ln(y) *Li(2, - y)*c1^3*[1-y]^-1 - 2*ln(y)*ln(y)*Li(2, - y)*c1^3 + 2*ln(y)* ln(y)*Li(2, - y)*d1^2*f1*[1-y]^-1 - ln(y)*ln(y)*Li(2, - y)*d1^2*f1 - 2*ln(y)*ln(y)*Li(2, - y)*d1^3*[1-y]^-1 + ln(y)*ln(y)*Li(2, - y)*d1^3 - ln(y)*ln(y)*Li(2,y)*a1*a2*c1 + ln(y)*ln(y)*Li(2,y)*a1*a2*d1 + ln(y )*ln(y)*Li(2,y)*a1*a2*f1 + 2*ln(y)*ln(y)*Li(2,y)*a1*c1*d1*[1-y]^-1 - 2*ln(y)*ln(y)*Li(2,y)*a1*c1*d1 + 2*ln(y)*ln(y)*Li(2,y)*a1*c1*f1* [1-y]^-1 - 2*ln(y)*ln(y)*Li(2,y)*a1*c1*f1 - 2*ln(y)*ln(y)*Li(2,y)*a1* c1^2*[1-y]^-1 + 2*ln(y)*ln(y)*Li(2,y)*a1*c1^2 - 2*ln(y)*ln(y)*Li(2,y) *a1*d1*f1*[1-y]^-1 + 2*ln(y)*ln(y)*Li(2,y)*a1*d1*f1 + 2*ln(y)*ln(y)* Li(2,y)*a2*c1*d1*[1-y]^-1 - 2*ln(y)*ln(y)*Li(2,y)*a2*c1*d1 + 2*ln(y)* ln(y)*Li(2,y)*a2*c1*f1*[1-y]^-1 - 2*ln(y)*ln(y)*Li(2,y)*a2*c1*f1 - 2* ln(y)*ln(y)*Li(2,y)*a2*c1^2*[1-y]^-1 + 2*ln(y)*ln(y)*Li(2,y)*a2*c1^2 - 2*ln(y)*ln(y)*Li(2,y)*a2*d1*f1*[1-y]^-1 + 2*ln(y)*ln(y)*Li(2,y)*a2 *d1*f1 + 4*ln(y)*ln(y)*Li(2,y)*c1*d1*f1*[1-y]^-1 - 3*ln(y)*ln(y)*Li(2 ,y)*c1*d1*f1 + 2*ln(y)*ln(y)*Li(2,y)*c1*d1^2*[1-y]^-1 - ln(y)*ln(y)* Li(2,y)*c1*d1^2 - 4*ln(y)*ln(y)*Li(2,y)*c1^2*d1*[1-y]^-1 + 3*ln(y)* ln(y)*Li(2,y)*c1^2*d1 - 2*ln(y)*ln(y)*Li(2,y)*c1^2*f1*[1-y]^-1 + 2* ln(y)*ln(y)*Li(2,y)*c1^2*f1 + 2*ln(y)*ln(y)*Li(2,y)*c1^3*[1-y]^-1 - 2 *ln(y)*ln(y)*Li(2,y)*c1^3 - 2*ln(y)*ln(y)*Li(2,y)*d1^2*f1*[1-y]^-1 + ln(y)*ln(y)*Li(2,y)*d1^2*f1 - 8*ln(y)*zt(3)*a1*a2*c1*[1-y]^-1 + 9*ln( y)*zt(3)*a1*a2*c1 + 8*ln(y)*zt(3)*a1*a2*d1*[1-y]^-1 - 9*ln(y)*zt(3)* a1*a2*d1 + 7*ln(y)*zt(3)*a1*a2*f1 - 3*ln(y)*zt(3)*a1*a2^2 - 10*ln(y)* zt(3)*a1*c1*d1*[1-y]^-1 + 10*ln(y)*zt(3)*a1*c1*d1 + 14*ln(y)*zt(3)*a1 *c1*f1*[1-y]^-1 - 14*ln(y)*zt(3)*a1*c1*f1 + 2*ln(y)*zt(3)*a1*c1^2* [1-y]^-1 - 3*ln(y)*zt(3)*a1*c1^2 - 14*ln(y)*zt(3)*a1*d1*f1*[1-y]^-1 + 14*ln(y)*zt(3)*a1*d1*f1 + 8*ln(y)*zt(3)*a1*d1^2*[1-y]^-1 - 7*ln(y) *zt(3)*a1*d1^2 - 7*ln(y)*zt(3)*a1*f1^2 - 3*ln(y)*zt(3)*a1^2*a2 - 4* ln(y)*zt(3)*a1^2*c1*[1-y]^-1 + 5*ln(y)*zt(3)*a1^2*c1 + 4*ln(y)*zt(3)* a1^2*d1*[1-y]^-1 - 5*ln(y)*zt(3)*a1^2*d1 + 7*ln(y)*zt(3)*a1^2*f1 - 2* ln(y)*zt(3)*a1^3 - 10*ln(y)*zt(3)*a2*c1*d1*[1-y]^-1 + 10*ln(y)*zt(3)* a2*c1*d1 + 14*ln(y)*zt(3)*a2*c1*f1*[1-y]^-1 - 14*ln(y)*zt(3)*a2*c1*f1 + 2*ln(y)*zt(3)*a2*c1^2*[1-y]^-1 - 3*ln(y)*zt(3)*a2*c1^2 - 14*ln(y)* zt(3)*a2*d1*f1*[1-y]^-1 + 14*ln(y)*zt(3)*a2*d1*f1 + 8*ln(y)*zt(3)*a2* d1^2*[1-y]^-1 - 7*ln(y)*zt(3)*a2*d1^2 - 7*ln(y)*zt(3)*a2*f1^2 - 4*ln( y)*zt(3)*a2^2*c1*[1-y]^-1 + 5*ln(y)*zt(3)*a2^2*c1 + 4*ln(y)*zt(3)* a2^2*d1*[1-y]^-1 - 5*ln(y)*zt(3)*a2^2*d1 + 7*ln(y)*zt(3)*a2^2*f1 - 2* ln(y)*zt(3)*a2^3 + 14*ln(y)*zt(3)*c1*d1*f1*[1-y]^-1 - 14*ln(y)*zt(3)* c1*d1*f1 - 6*ln(y)*zt(3)*c1*d1^2*[1-y]^-1 + 6*ln(y)*zt(3)*c1*d1^2 - 14*ln(y)*zt(3)*c1*f1^2*[1-y]^-1 + 14*ln(y)*zt(3)*c1*f1^2 + 4*ln(y)* zt(3)*c1^2*d1*[1-y]^-1 - 4*ln(y)*zt(3)*c1^2*d1 - 2*ln(y)*zt(3)*c1^3* [1-y]^-1 + 2*ln(y)*zt(3)*c1^3 + 14*ln(y)*zt(3)*d1*f1^2*[1-y]^-1 - 14* ln(y)*zt(3)*d1*f1^2 - 14*ln(y)*zt(3)*d1^2*f1*[1-y]^-1 + 14*ln(y)*zt(3 )*d1^2*f1 + 4*ln(y)*zt(3)*d1^3*[1-y]^-1 - 4*ln(y)*zt(3)*d1^3 - 8*ln(y )*Li(3, - y)*a1*a2*c1*[1-y]^-1 + 16*ln(y)*Li(3, - y)*a1*a2*c1 + 8*ln( y)*Li(3, - y)*a1*a2*d1*[1-y]^-1 - 16*ln(y)*Li(3, - y)*a1*a2*d1 + 4* ln(y)*Li(3, - y)*a1*a2*f1 - 4*ln(y)*Li(3, - y)*a1*a2^2 - 20*ln(y)*Li( 3, - y)*a1*c1*d1*[1-y]^-1 + 22*ln(y)*Li(3, - y)*a1*c1*d1 + 4*ln(y)* Li(3, - y)*a1*c1*f1*[1-y]^-1 - 6*ln(y)*Li(3, - y)*a1*c1*f1 + 12*ln(y) *Li(3, - y)*a1*c1^2*[1-y]^-1 - 14*ln(y)*Li(3, - y)*a1*c1^2 - 4*ln(y)* Li(3, - y)*a1*d1*f1*[1-y]^-1 + 6*ln(y)*Li(3, - y)*a1*d1*f1 + 8*ln(y)* Li(3, - y)*a1*d1^2*[1-y]^-1 - 8*ln(y)*Li(3, - y)*a1*d1^2 - 4*ln(y)* Li(3, - y)*a1^2*a2 - 4*ln(y)*Li(3, - y)*a1^2*c1*[1-y]^-1 + 6*ln(y)* Li(3, - y)*a1^2*c1 + 4*ln(y)*Li(3, - y)*a1^2*d1*[1-y]^-1 - 6*ln(y)* Li(3, - y)*a1^2*d1 - 20*ln(y)*Li(3, - y)*a2*c1*d1*[1-y]^-1 + 22*ln(y) *Li(3, - y)*a2*c1*d1 + 4*ln(y)*Li(3, - y)*a2*c1*f1*[1-y]^-1 - 6*ln(y) *Li(3, - y)*a2*c1*f1 + 12*ln(y)*Li(3, - y)*a2*c1^2*[1-y]^-1 - 14*ln(y )*Li(3, - y)*a2*c1^2 - 4*ln(y)*Li(3, - y)*a2*d1*f1*[1-y]^-1 + 6*ln(y) *Li(3, - y)*a2*d1*f1 + 8*ln(y)*Li(3, - y)*a2*d1^2*[1-y]^-1 - 8*ln(y)* Li(3, - y)*a2*d1^2 - 4*ln(y)*Li(3, - y)*a2^2*c1*[1-y]^-1 + 6*ln(y)* Li(3, - y)*a2^2*c1 + 4*ln(y)*Li(3, - y)*a2^2*d1*[1-y]^-1 - 6*ln(y)* Li(3, - y)*a2^2*d1 + 12*ln(y)*Li(3, - y)*c1*d1*f1*[1-y]^-1 - 10*ln(y) *Li(3, - y)*c1*d1*f1 - 16*ln(y)*Li(3, - y)*c1*d1^2*[1-y]^-1 + 12*ln(y )*Li(3, - y)*c1*d1^2 + 20*ln(y)*Li(3, - y)*c1^2*d1*[1-y]^-1 - 18*ln(y )*Li(3, - y)*c1^2*d1 - 8*ln(y)*Li(3, - y)*c1^2*f1*[1-y]^-1 + 8*ln(y)* Li(3, - y)*c1^2*f1 - 8*ln(y)*Li(3, - y)*c1^3*[1-y]^-1 + 8*ln(y)*Li(3, - y)*c1^3 - 4*ln(y)*Li(3, - y)*d1^2*f1*[1-y]^-1 + 2*ln(y)*Li(3, - y) *d1^2*f1 + 4*ln(y)*Li(3, - y)*d1^3*[1-y]^-1 - 2*ln(y)*Li(3, - y)*d1^3 + 4*ln(y)*Li(3,y)*a1*a2*c1 - 4*ln(y)*Li(3,y)*a1*a2*d1 - 4*ln(y)*Li(3 ,y)*a1*a2*f1 - 4*ln(y)*Li(3,y)*a1*c1*d1*[1-y]^-1 + 6*ln(y)*Li(3,y)*a1 *c1*d1 - 4*ln(y)*Li(3,y)*a1*c1*f1*[1-y]^-1 + 6*ln(y)*Li(3,y)*a1*c1*f1 + 4*ln(y)*Li(3,y)*a1*c1^2*[1-y]^-1 - 6*ln(y)*Li(3,y)*a1*c1^2 + 4*ln( y)*Li(3,y)*a1*d1*f1*[1-y]^-1 - 6*ln(y)*Li(3,y)*a1*d1*f1 - 4*ln(y)*Li( 3,y)*a2*c1*d1*[1-y]^-1 + 6*ln(y)*Li(3,y)*a2*c1*d1 - 4*ln(y)*Li(3,y)* a2*c1*f1*[1-y]^-1 + 6*ln(y)*Li(3,y)*a2*c1*f1 + 4*ln(y)*Li(3,y)*a2* c1^2*[1-y]^-1 - 6*ln(y)*Li(3,y)*a2*c1^2 + 4*ln(y)*Li(3,y)*a2*d1*f1* [1-y]^-1 - 6*ln(y)*Li(3,y)*a2*d1*f1 - 12*ln(y)*Li(3,y)*c1*d1*f1* [1-y]^-1 + 10*ln(y)*Li(3,y)*c1*d1*f1 - 4*ln(y)*Li(3,y)*c1*d1^2* [1-y]^-1 + 2*ln(y)*Li(3,y)*c1*d1^2 + 12*ln(y)*Li(3,y)*c1^2*d1* [1-y]^-1 - 10*ln(y)*Li(3,y)*c1^2*d1 + 8*ln(y)*Li(3,y)*c1^2*f1* [1-y]^-1 - 8*ln(y)*Li(3,y)*c1^2*f1 - 8*ln(y)*Li(3,y)*c1^3*[1-y]^-1 + 8*ln(y)*Li(3,y)*c1^3 + 4*ln(y)*Li(3,y)*d1^2*f1*[1-y]^-1 - 2*ln(y)*Li( 3,y)*d1^2*f1 - 8*zt(2)*Li(2, - y)*a1*a2*c1*[1-y]^-1 + 4*zt(2)*Li(2, - y)*a1*a2*c1 + 8*zt(2)*Li(2, - y)*a1*a2*d1*[1-y]^-1 - 4*zt(2)*Li(2, - y)*a1*a2*d1 - 16*zt(2)*Li(2, - y)*a1*c1*d1*[1-y]^-1 + 8*zt(2)*Li(2 , - y)*a1*c1*d1 + 16*zt(2)*Li(2, - y)*a1*c1*f1*[1-y]^-1 - 8*zt(2)*Li( 2, - y)*a1*c1*f1 + 8*zt(2)*Li(2, - y)*a1*c1^2*[1-y]^-1 - 4*zt(2)*Li(2 , - y)*a1*c1^2 - 16*zt(2)*Li(2, - y)*a1*d1*f1*[1-y]^-1 + 8*zt(2)*Li(2 , - y)*a1*d1*f1 + 8*zt(2)*Li(2, - y)*a1*d1^2*[1-y]^-1 - 4*zt(2)*Li(2, - y)*a1*d1^2 - 4*zt(2)*Li(2, - y)*a1^2*c1*[1-y]^-1 + 2*zt(2)*Li(2, - y)*a1^2*c1 + 4*zt(2)*Li(2, - y)*a1^2*d1*[1-y]^-1 - 2*zt(2)*Li(2, - y)*a1^2*d1 - 16*zt(2)*Li(2, - y)*a2*c1*d1*[1-y]^-1 + 8*zt(2)*Li(2, - y)*a2*c1*d1 + 16*zt(2)*Li(2, - y)*a2*c1*f1*[1-y]^-1 - 8*zt(2)*Li(2 , - y)*a2*c1*f1 + 8*zt(2)*Li(2, - y)*a2*c1^2*[1-y]^-1 - 4*zt(2)*Li(2, - y)*a2*c1^2 - 16*zt(2)*Li(2, - y)*a2*d1*f1*[1-y]^-1 + 8*zt(2)*Li(2, - y)*a2*d1*f1 + 8*zt(2)*Li(2, - y)*a2*d1^2*[1-y]^-1 - 4*zt(2)*Li(2, - y)*a2*d1^2 - 4*zt(2)*Li(2, - y)*a2^2*c1*[1-y]^-1 + 2*zt(2)*Li(2, - y)*a2^2*c1 + 4*zt(2)*Li(2, - y)*a2^2*d1*[1-y]^-1 - 2*zt(2)*Li(2, - y)*a2^2*d1 + 32*zt(2)*Li(2, - y)*c1*d1*f1*[1-y]^-1 - 16*zt(2)*Li(2 , - y)*c1*d1*f1 - 12*zt(2)*Li(2, - y)*c1*d1^2*[1-y]^-1 + 6*zt(2)*Li(2 , - y)*c1*d1^2 - 12*zt(2)*Li(2, - y)*c1*f1^2*[1-y]^-1 + 6*zt(2)*Li(2, - y)*c1*f1^2 + 12*zt(2)*Li(2, - y)*c1^2*d1*[1-y]^-1 - 6*zt(2)*Li(2, - y)*c1^2*d1 - 16*zt(2)*Li(2, - y)*c1^2*f1*[1-y]^-1 + 8*zt(2)*Li(2, - y)*c1^2*f1 - 4*zt(2)*Li(2, - y)*c1^3*[1-y]^-1 + 2*zt(2)*Li(2, - y) *c1^3 + 12*zt(2)*Li(2, - y)*d1*f1^2*[1-y]^-1 - 6*zt(2)*Li(2, - y)*d1* f1^2 - 16*zt(2)*Li(2, - y)*d1^2*f1*[1-y]^-1 + 8*zt(2)*Li(2, - y)*d1^2 *f1 + 4*zt(2)*Li(2, - y)*d1^3*[1-y]^-1 - 2*zt(2)*Li(2, - y)*d1^3 - 4* zt(2)*Li(2,y)*a1*c1*d1*[1-y]^-1 + 2*zt(2)*Li(2,y)*a1*c1*d1 - 4*zt(2)* Li(2,y)*a1*c1*f1*[1-y]^-1 + 2*zt(2)*Li(2,y)*a1*c1*f1 + 4*zt(2)*Li(2,y )*a1*c1^2*[1-y]^-1 - 2*zt(2)*Li(2,y)*a1*c1^2 + 4*zt(2)*Li(2,y)*a1*d1* f1*[1-y]^-1 - 2*zt(2)*Li(2,y)*a1*d1*f1 - 4*zt(2)*Li(2,y)*a2*c1*d1* [1-y]^-1 + 2*zt(2)*Li(2,y)*a2*c1*d1 - 4*zt(2)*Li(2,y)*a2*c1*f1* [1-y]^-1 + 2*zt(2)*Li(2,y)*a2*c1*f1 + 4*zt(2)*Li(2,y)*a2*c1^2* [1-y]^-1 - 2*zt(2)*Li(2,y)*a2*c1^2 + 4*zt(2)*Li(2,y)*a2*d1*f1* [1-y]^-1 - 2*zt(2)*Li(2,y)*a2*d1*f1 + 4*zt(2)*Li(2,y)*c1*d1*f1* [1-y]^-1 - 2*zt(2)*Li(2,y)*c1*d1*f1 - 4*zt(2)*Li(2,y)*c1*d1^2* [1-y]^-1 + 2*zt(2)*Li(2,y)*c1*d1^2 + 12*zt(2)*Li(2,y)*c1*f1^2* [1-y]^-1 - 6*zt(2)*Li(2,y)*c1*f1^2 + 8*zt(2)*Li(2,y)*c1^2*d1*[1-y]^-1 - 4*zt(2)*Li(2,y)*c1^2*d1 - 8*zt(2)*Li(2,y)*c1^2*f1*[1-y]^-1 + 4*zt( 2)*Li(2,y)*c1^2*f1 - 4*zt(2)*Li(2,y)*c1^3*[1-y]^-1 + 2*zt(2)*Li(2,y)* c1^3 - 12*zt(2)*Li(2,y)*d1*f1^2*[1-y]^-1 + 6*zt(2)*Li(2,y)*d1*f1^2 + 4*zt(2)*Li(2,y)*d1^2*f1*[1-y]^-1 - 2*zt(2)*Li(2,y)*d1^2*f1 - 5*zt(4)* a1*a2*c1*[1-y]^-1 + zt(4)*a1*a2*c1 + 5*zt(4)*a1*a2*d1*[1-y]^-1 - zt(4 )*a1*a2*d1 - 15*zt(4)*a1*a2*f1 + 1/4*zt(4)*a1*a2^2 - 2*zt(4)*a1*c1*d1 *[1-y]^-1 - 1/2*zt(4)*a1*c1*d1 + 15*zt(4)*a1*c1*f1*[1-y]^-1 + 15/2* zt(4)*a1*c1*f1 - 3*zt(4)*a1*c1^2*[1-y]^-1 + 11/4*zt(4)*a1*c1^2 - 15* zt(4)*a1*d1*f1*[1-y]^-1 - 15/2*zt(4)*a1*d1*f1 + 5*zt(4)*a1*d1^2* [1-y]^-1 - 9/4*zt(4)*a1*d1^2 + 45/4*zt(4)*a1*f1^2 + 1/4*zt(4)*a1^2*a2 - 5/2*zt(4)*a1^2*c1*[1-y]^-1 + zt(4)*a1^2*c1 + 5/2*zt(4)*a1^2*d1* [1-y]^-1 - zt(4)*a1^2*d1 - 9/4*zt(4)*a1^3 - 2*zt(4)*a2*c1*d1*[1-y]^-1 - 1/2*zt(4)*a2*c1*d1 + 15*zt(4)*a2*c1*f1*[1-y]^-1 + 15/2*zt(4)*a2*c1 *f1 - 3*zt(4)*a2*c1^2*[1-y]^-1 + 11/4*zt(4)*a2*c1^2 - 15*zt(4)*a2*d1* f1*[1-y]^-1 - 15/2*zt(4)*a2*d1*f1 + 5*zt(4)*a2*d1^2*[1-y]^-1 - 9/4* zt(4)*a2*d1^2 + 45/4*zt(4)*a2*f1^2 - 5/2*zt(4)*a2^2*c1*[1-y]^-1 + zt( 4)*a2^2*c1 + 5/2*zt(4)*a2^2*d1*[1-y]^-1 - zt(4)*a2^2*d1 - 9/4*zt(4)* a2^3 + 15*zt(4)*c1*d1*f1 + 1/2*zt(4)*c1*d1^2*[1-y]^-1 - 1/2*zt(4)*c1* d1^2 - 45/2*zt(4)*c1*f1^2*[1-y]^-1 - 5/2*zt(4)*c1^2*d1*[1-y]^-1 + 5/2 *zt(4)*c1^2*d1 + 15*zt(4)*c1^2*f1*[1-y]^-1 - 45/2*zt(4)*c1^2*f1 - 1/2 *zt(4)*c1^3*[1-y]^-1 + 3/2*zt(4)*c1^3 + 45/2*zt(4)*d1*f1^2*[1-y]^-1 - 15*zt(4)*d1^2*f1*[1-y]^-1 + 15/2*zt(4)*d1^2*f1 + 5/2*zt(4)*d1^3* [1-y]^-1 - 7/2*zt(4)*d1^3 - 8*Li(2, - y)*Li(2, - y)*a1*a2*c1*[1-y]^-1 + 4*Li(2, - y)*Li(2, - y)*a1*a2*c1 + 8*Li(2, - y)*Li(2, - y)*a1*a2* d1*[1-y]^-1 - 4*Li(2, - y)*Li(2, - y)*a1*a2*d1 - 16*Li(2, - y)*Li(2, - y)*a1*c1*d1*[1-y]^-1 + 8*Li(2, - y)*Li(2, - y)*a1*c1*d1 + 8*Li(2, - y)*Li(2, - y)*a1*c1*f1*[1-y]^-1 - 4*Li(2, - y)*Li(2, - y)*a1*c1*f1 + 8*Li(2, - y)*Li(2, - y)*a1*c1^2*[1-y]^-1 - 4*Li(2, - y)*Li(2, - y) *a1*c1^2 - 8*Li(2, - y)*Li(2, - y)*a1*d1*f1*[1-y]^-1 + 4*Li(2, - y)* Li(2, - y)*a1*d1*f1 + 8*Li(2, - y)*Li(2, - y)*a1*d1^2*[1-y]^-1 - 4* Li(2, - y)*Li(2, - y)*a1*d1^2 - 4*Li(2, - y)*Li(2, - y)*a1^2*c1* [1-y]^-1 + 2*Li(2, - y)*Li(2, - y)*a1^2*c1 + 4*Li(2, - y)*Li(2, - y)* a1^2*d1*[1-y]^-1 - 2*Li(2, - y)*Li(2, - y)*a1^2*d1 - 16*Li(2, - y)* Li(2, - y)*a2*c1*d1*[1-y]^-1 + 8*Li(2, - y)*Li(2, - y)*a2*c1*d1 + 8* Li(2, - y)*Li(2, - y)*a2*c1*f1*[1-y]^-1 - 4*Li(2, - y)*Li(2, - y)*a2* c1*f1 + 8*Li(2, - y)*Li(2, - y)*a2*c1^2*[1-y]^-1 - 4*Li(2, - y)*Li(2, - y)*a2*c1^2 - 8*Li(2, - y)*Li(2, - y)*a2*d1*f1*[1-y]^-1 + 4*Li(2, - y)*Li(2, - y)*a2*d1*f1 + 8*Li(2, - y)*Li(2, - y)*a2*d1^2*[1-y]^-1 - 4*Li(2, - y)*Li(2, - y)*a2*d1^2 - 4*Li(2, - y)*Li(2, - y)*a2^2*c1* [1-y]^-1 + 2*Li(2, - y)*Li(2, - y)*a2^2*c1 + 4*Li(2, - y)*Li(2, - y)* a2^2*d1*[1-y]^-1 - 2*Li(2, - y)*Li(2, - y)*a2^2*d1 + 16*Li(2, - y)* Li(2, - y)*c1*d1*f1*[1-y]^-1 - 8*Li(2, - y)*Li(2, - y)*c1*d1*f1 - 12* Li(2, - y)*Li(2, - y)*c1*d1^2*[1-y]^-1 + 6*Li(2, - y)*Li(2, - y)*c1* d1^2 - 4*Li(2, - y)*Li(2, - y)*c1*f1^2*[1-y]^-1 + 2*Li(2, - y)*Li(2, - y)*c1*f1^2 + 12*Li(2, - y)*Li(2, - y)*c1^2*d1*[1-y]^-1 - 6*Li(2, - y)*Li(2, - y)*c1^2*d1 - 8*Li(2, - y)*Li(2, - y)*c1^2*f1*[1-y]^-1 + 4*Li(2, - y)*Li(2, - y)*c1^2*f1 - 4*Li(2, - y)*Li(2, - y)*c1^3* [1-y]^-1 + 2*Li(2, - y)*Li(2, - y)*c1^3 + 4*Li(2, - y)*Li(2, - y)*d1* f1^2*[1-y]^-1 - 2*Li(2, - y)*Li(2, - y)*d1*f1^2 - 8*Li(2, - y)*Li(2, - y)*d1^2*f1*[1-y]^-1 + 4*Li(2, - y)*Li(2, - y)*d1^2*f1 + 4*Li(2, - y)*Li(2, - y)*d1^3*[1-y]^-1 - 2*Li(2, - y)*Li(2, - y)*d1^3 - 8*Li(2,y )*Li(2, - y)*a1*c1*f1*[1-y]^-1 + 4*Li(2,y)*Li(2, - y)*a1*c1*f1 + 8* Li(2,y)*Li(2, - y)*a1*d1*f1*[1-y]^-1 - 4*Li(2,y)*Li(2, - y)*a1*d1*f1 - 8*Li(2,y)*Li(2, - y)*a2*c1*f1*[1-y]^-1 + 4*Li(2,y)*Li(2, - y)*a2* c1*f1 + 8*Li(2,y)*Li(2, - y)*a2*d1*f1*[1-y]^-1 - 4*Li(2,y)*Li(2, - y) *a2*d1*f1 - 16*Li(2,y)*Li(2, - y)*c1*d1*f1*[1-y]^-1 + 8*Li(2,y)*Li(2, - y)*c1*d1*f1 + 8*Li(2,y)*Li(2, - y)*c1*f1^2*[1-y]^-1 - 4*Li(2,y)* Li(2, - y)*c1*f1^2 + 8*Li(2,y)*Li(2, - y)*c1^2*f1*[1-y]^-1 - 4*Li(2,y )*Li(2, - y)*c1^2*f1 - 8*Li(2,y)*Li(2, - y)*d1*f1^2*[1-y]^-1 + 4*Li(2 ,y)*Li(2, - y)*d1*f1^2 + 8*Li(2,y)*Li(2, - y)*d1^2*f1*[1-y]^-1 - 4* Li(2,y)*Li(2, - y)*d1^2*f1 - 4*Li(2,y)*Li(2,y)*c1*d1*f1*[1-y]^-1 + 2* Li(2,y)*Li(2,y)*c1*d1*f1 - 4*Li(2,y)*Li(2,y)*c1*f1^2*[1-y]^-1 + 2*Li( 2,y)*Li(2,y)*c1*f1^2 + 4*Li(2,y)*Li(2,y)*c1^2*f1*[1-y]^-1 - 2*Li(2,y) *Li(2,y)*c1^2*f1 + 4*Li(2,y)*Li(2,y)*d1*f1^2*[1-y]^-1 - 2*Li(2,y)*Li( 2,y)*d1*f1^2 - 12*Li(4, - y)*a1*a2*c1 + 12*Li(4, - y)*a1*a2*d1 - 4* Li(4, - y)*a1*a2*f1 + 2*Li(4, - y)*a1*a2^2 + 12*Li(4, - y)*a1*c1*d1* [1-y]^-1 - 18*Li(4, - y)*a1*c1*d1 + 4*Li(4, - y)*a1*c1*f1 - 12*Li(4, - y)*a1*c1^2*[1-y]^-1 + 16*Li(4, - y)*a1*c1^2 - 4*Li(4, - y)*a1*d1* f1 + 2*Li(4, - y)*a1*d1^2 + 2*Li(4, - y)*a1^2*a2 - 2*Li(4, - y)*a1^2* c1 + 2*Li(4, - y)*a1^2*d1 + 2*Li(4, - y)*a1^2*f1 - 2*Li(4, - y)*a1^3 + 12*Li(4, - y)*a2*c1*d1*[1-y]^-1 - 18*Li(4, - y)*a2*c1*d1 + 4*Li(4, - y)*a2*c1*f1 - 12*Li(4, - y)*a2*c1^2*[1-y]^-1 + 16*Li(4, - y)*a2* c1^2 - 4*Li(4, - y)*a2*d1*f1 + 2*Li(4, - y)*a2*d1^2 - 2*Li(4, - y)* a2^2*c1 + 2*Li(4, - y)*a2^2*d1 + 2*Li(4, - y)*a2^2*f1 - 2*Li(4, - y)* a2^3 - 12*Li(4, - y)*c1*d1*f1*[1-y]^-1 + 10*Li(4, - y)*c1*d1*f1 + 12* Li(4, - y)*c1*d1^2*[1-y]^-1 - 8*Li(4, - y)*c1*d1^2 - 24*Li(4, - y)* c1^2*d1*[1-y]^-1 + 22*Li(4, - y)*c1^2*d1 + 12*Li(4, - y)*c1^2*f1* [1-y]^-1 - 12*Li(4, - y)*c1^2*f1 + 12*Li(4, - y)*c1^3*[1-y]^-1 - 12* Li(4, - y)*c1^3 + 2*Li(4, - y)*d1^2*f1 - 2*Li(4, - y)*d1^3 - 8*Li(4, - [1-y]*[1+y]^-1)*a1*a2*f1 + 8*Li(4, - [1-y]*[1+y]^-1)*a1*c1*f1 - 8* Li(4, - [1-y]*[1+y]^-1)*a1*d1*f1 + 12*Li(4, - [1-y]*[1+y]^-1)*a1*f1^2 - 4*Li(4, - [1-y]*[1+y]^-1)*a1^2*f1 + 8*Li(4, - [1-y]*[1+y]^-1)*a2* c1*f1 - 8*Li(4, - [1-y]*[1+y]^-1)*a2*d1*f1 + 12*Li(4, - [1-y]* [1+y]^-1)*a2*f1^2 - 4*Li(4, - [1-y]*[1+y]^-1)*a2^2*f1 + 8*Li(4, - [1-y]*[1+y]^-1)*c1*d1*f1 - 12*Li(4, - [1-y]*[1+y]^-1)*c1*f1^2 - 4*Li( 4, - [1-y]*[1+y]^-1)*c1^2*f1 + 12*Li(4, - [1-y]*[1+y]^-1)*d1*f1^2 - 4 *Li(4, - [1-y]*[1+y]^-1)*d1^2*f1 - 8*Li(4, - [1-y]*[1+y]^-1)*f1^3 + 8 *Li(4,[1-y]*[1+y]^-1)*a1*a2*f1 - 8*Li(4,[1-y]*[1+y]^-1)*a1*c1*f1 + 8* Li(4,[1-y]*[1+y]^-1)*a1*d1*f1 - 12*Li(4,[1-y]*[1+y]^-1)*a1*f1^2 + 4* Li(4,[1-y]*[1+y]^-1)*a1^2*f1 - 8*Li(4,[1-y]*[1+y]^-1)*a2*c1*f1 + 8* Li(4,[1-y]*[1+y]^-1)*a2*d1*f1 - 12*Li(4,[1-y]*[1+y]^-1)*a2*f1^2 + 4* Li(4,[1-y]*[1+y]^-1)*a2^2*f1 - 8*Li(4,[1-y]*[1+y]^-1)*c1*d1*f1 + 12* Li(4,[1-y]*[1+y]^-1)*c1*f1^2 + 4*Li(4,[1-y]*[1+y]^-1)*c1^2*f1 - 12* Li(4,[1-y]*[1+y]^-1)*d1*f1^2 + 4*Li(4,[1-y]*[1+y]^-1)*d1^2*f1 + 8*Li( 4,[1-y]*[1+y]^-1)*f1^3 - 6*Li(4,y)*a1*a2*c1 + 6*Li(4,y)*a1*a2*d1 + 4* Li(4,y)*a1*a2*f1 + 4*Li(4,y)*a1*c1*d1*[1-y]^-1 - 8*Li(4,y)*a1*c1*d1 - 4*Li(4,y)*a1*c1*f1 - 4*Li(4,y)*a1*c1^2*[1-y]^-1 + 8*Li(4,y)*a1* c1^2 + 4*Li(4,y)*a1*d1*f1 - 2*Li(4,y)*a1^2*f1 + 4*Li(4,y)*a2*c1*d1* [1-y]^-1 - 8*Li(4,y)*a2*c1*d1 - 4*Li(4,y)*a2*c1*f1 - 4*Li(4,y)*a2* c1^2*[1-y]^-1 + 8*Li(4,y)*a2*c1^2 + 4*Li(4,y)*a2*d1*f1 - 2*Li(4,y)* a2^2*f1 + 12*Li(4,y)*c1*d1*f1*[1-y]^-1 - 10*Li(4,y)*c1*d1*f1 + 4*Li(4 ,y)*c1*d1^2*[1-y]^-1 - 2*Li(4,y)*c1*d1^2 - 16*Li(4,y)*c1^2*d1* [1-y]^-1 + 14*Li(4,y)*c1^2*d1 - 12*Li(4,y)*c1^2*f1*[1-y]^-1 + 12*Li(4 ,y)*c1^2*f1 + 12*Li(4,y)*c1^3*[1-y]^-1 - 12*Li(4,y)*c1^3 - 2*Li(4,y)* d1^2*f1 ) + ep^4*F32C(a1,a2,d1,f1,c1,z,ep^4) ) ;