sys32 = {(-3*n)/2 - omd - Omd + S1/C - G*m*M*(C22*((21*e*R^2*Cos[i])/(4*a^3*S1) - (369*e^3*R^2*Cos[i])/(32*a^3*S1) - (21*e*R^2*(S1 - S3)*Cos[i])/(4*a^3*S1^2) + (369*e^3*R^2*(S1 - S3)*Cos[i])/(32*a^3*S1^2) + (21*e*R^2*(S1 - S3)*Cos[i]^2)/(4*a^3*S1^2) - (369*e^3*R^2*(S1 - S3)*Cos[i]^2)/(32*a^3*S1^2) - (21*e*R^2*(S1 - S3)^2*Cos[i]^2)/(4*a^3*S1^3) + (369*e^3*R^2*(S1 - S3)^2*Cos[i]^2)/(32*a^3*S1^3) + (21*e*R^2*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i])/(8*a^3*Sqrt[1 - (S1 - S3)^2/S1^2]) - (369*e^3*R^2*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i])/(64*a^3*Sqrt[1 - (S1 - S3)^2/S1^2]) + (21*e*R^2*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(4*a^3*S1) - (369*e^3*R^2*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(32*a^3*S1) + (21*e*R^2*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i])/(8*a^3*S1*Sqrt[1 - (S1 - S3)^2/S1^2]) - (369*e^3*R^2*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i])/(64*a^3*S1*Sqrt[1 - (S1 - S3)^2/S1^2]) - (21*e*R^2*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i])/(4*a^3*S1^2) + (369*e^3*R^2*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i])/(32*a^3*S1^2) + (21*e*R^2*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i]^2)/(8*a^3) - (369*e^3*R^2*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i]^2)/(64*a^3)) + C20*((3*R^2*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Cos[i]^2)/(4*a^3) + (9*e^2*R^2*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Cos[i]^2)/(8*a^3) + (45*e^4*R^2*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Cos[i]^2)/(32*a^3) - (3*R^2*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(2*a^3*S1) - (9*e^2*R^2*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(4*a^3*S1) - (45*e^4*R^2*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(16*a^3*S1) - (3*R^2*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i])/(4*a^3*S1*Sqrt[1 - (S1 - S3)^2/S1^2]) - (9*e^2*R^2*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i])/(8*a^3*S1*Sqrt[1 - (S1 - S3)^2/S1^2]) - (45*e^4*R^2*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i])/(32*a^3*S1*Sqrt[1 - (S1 - S3)^2/S1^2]) + (3*R^2*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i])/(2*a^3*S1^2) + (9*e^2*R^2*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i])/(4*a^3*S1^2) + (45*e^4*R^2*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i])/(16*a^3*S1^2) + (3*R^2*(S1 - S3)*Sin[i]^2)/(2*a^3*S1^2) + (9*e^2*R^2*(S1 - S3)*Sin[i]^2)/(4*a^3*S1^2) + (45*e^4*R^2*(S1 - S3)*Sin[i]^2)/(16*a^3*S1^2) - (3*R^2*(S1 - S3)^2*Sin[i]^2)/(2*a^3*S1^3) - (9*e^2*R^2*(S1 - S3)^2*Sin[i]^2)/(4*a^3*S1^3) - (45*e^4*R^2*(S1 - S3)^2*Sin[i]^2)/(16*a^3*S1^3)) + C44*((945*e*R^4*Cos[i])/(2*a^5*S1) + (5355*e^2*R^4*Cos[i])/(8*a^5*S1) - (315*e^3*R^4*Cos[i])/(16*a^5*S1) - (33705*e^4*R^4*Cos[i])/(8*a^5*S1) - (945*e*R^4*(S1 - S3)*Cos[i])/(2*a^5*S1^2) - (5355*e^2*R^4*(S1 - S3)*Cos[i])/(8*a^5*S1^2) + (315*e^3*R^4*(S1 - S3)*Cos[i])/(16*a^5*S1^2) + (33705*e^4*R^4*(S1 - S3)*Cos[i])/(8*a^5*S1^2) - (105*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Cos[i]^2)/(2*a^5) - (525*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Cos[i]^2)/(2*a^5) - (11025*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Cos[i]^2)/(16*a^5) - (525*R^4*(S1 - S3)*Cos[i]^2)/(2*a^5*S1^2) + (945*e*R^4*(S1 - S3)*Cos[i]^2)/(a^5*S1^2) + (5565*e^2*R^4*(S1 - S3)*Cos[i]^2)/(8*a^5*S1^2) - (315*e^3*R^4*(S1 - S3)*Cos[i]^2)/(8*a^5*S1^2) - (257355*e^4*R^4*(S1 - S3)*Cos[i]^2)/(16*a^5*S1^2) + (525*R^4*(S1 - S3)^2*Cos[i]^2)/(2*a^5*S1^3) - (945*e*R^4*(S1 - S3)^2*Cos[i]^2)/(a^5*S1^3) - (5565*e^2*R^4*(S1 - S3)^2*Cos[i]^2)/(8*a^5*S1^3) + (315*e^3*R^4*(S1 - S3)^2*Cos[i]^2)/(8*a^5*S1^3) + (257355*e^4*R^4*(S1 - S3)^2*Cos[i]^2)/(16*a^5*S1^3) - (945*e*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]^3)/(2*a^5*S1) + (315*e^3*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]^3)/(16*a^5*S1) - (945*e*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]^3)/(2*a^5*S1) + (315*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]^3)/(16*a^5*S1) + (945*e*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^3)/(2*a^5*S1^2) - (315*e^3*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^3)/(16*a^5*S1^2) - (2835*e*R^4*(S1 - S3)^2*Cos[i]^3)/(2*a^5*S1^3) + (16065*e^2*R^4*(S1 - S3)^2*Cos[i]^3)/(8*a^5*S1^3) + (945*e^3*R^4*(S1 - S3)^2*Cos[i]^3)/(16*a^5*S1^3) - (101115*e^4*R^4*(S1 - S3)^2*Cos[i]^3)/(8*a^5*S1^3) + (2835*e*R^4*(S1 - S3)^3*Cos[i]^3)/(2*a^5*S1^4) - (16065*e^2*R^4*(S1 - S3)^3*Cos[i]^3)/(8*a^5*S1^4) - (945*e^3*R^4*(S1 - S3)^3*Cos[i]^3)/(16*a^5*S1^4) + (101115*e^4*R^4*(S1 - S3)^3*Cos[i]^3)/(8*a^5*S1^4) + (315*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(1 - (S1 - S3)^2/S1^2)*Cos[i]^4)/(4*a^5) + (1575*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(1 - (S1 - S3)^2/S1^2)*Cos[i]^4)/(4*a^5) + (33075*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(1 - (S1 - S3)^2/S1^2)*Cos[i]^4)/(32*a^5) + (315*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^4)/(a^5*S1^2) - (945*e*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^4)/(a^5*S1^2) + (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^4)/(a^5*S1^2) + (315*e^3*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^4)/(8*a^5*S1^2) + (33075*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^4)/(8*a^5*S1^2) + (315*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^4)/(2*a^5*S1^2) - (945*e*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^4)/(2*a^5*S1^2) + (1575*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^4)/(2*a^5*S1^2) + (315*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^4)/(16*a^5*S1^2) + (33075*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^4)/(16*a^5*S1^2) - (315*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Cos[i]^4)/(a^5*S1^3) + (945*e*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Cos[i]^4)/(a^5*S1^3) - (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Cos[i]^4)/(a^5*S1^3) - (315*e^3*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Cos[i]^4)/(8*a^5*S1^3) - (33075*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Cos[i]^4)/(8*a^5*S1^3) + (945*R^4*(S1 - S3)^3*Cos[i]^4)/(2*a^5*S1^4) - (1890*e*R^4*(S1 - S3)^3*Cos[i]^4)/(a^5*S1^4) + (24255*e^2*R^4*(S1 - S3)^3*Cos[i]^4)/(8*a^5*S1^4) + (315*e^3*R^4*(S1 - S3)^3*Cos[i]^4)/(4*a^5*S1^4) + (31815*e^4*R^4*(S1 - S3)^3*Cos[i]^4)/(16*a^5*S1^4) - (945*R^4*(S1 - S3)^4*Cos[i]^4)/(2*a^5*S1^5) + (1890*e*R^4*(S1 - S3)^4*Cos[i]^4)/(a^5*S1^5) - (24255*e^2*R^4*(S1 - S3)^4*Cos[i]^4)/(8*a^5*S1^5) - (315*e^3*R^4*(S1 - S3)^4*Cos[i]^4)/(4*a^5*S1^5) - (31815*e^4*R^4*(S1 - S3)^4*Cos[i]^4)/(16*a^5*S1^5) + (945*e*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i])/(4*a^5*Sqrt[1 - (S1 - S3)^2/S1^2]) + (5355*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i])/(16*a^5*Sqrt[1 - (S1 - S3)^2/S1^2]) - (315*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i])/(32*a^5*Sqrt[1 - (S1 - S3)^2/S1^2]) - (33705*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i])/(16*a^5*Sqrt[1 - (S1 - S3)^2/S1^2]) - (315*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(2*a^5*S1) + (945*e*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(a^5*S1) + (9765*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(8*a^5*S1) - (315*e^3*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(8*a^5*S1) - (235305*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(16*a^5*S1) - (315*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i])/(4*a^5*S1*Sqrt[1 - (S1 - S3)^2/S1^2]) + (945*e*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i])/(2*a^5*S1*Sqrt[1 - (S1 - S3)^2/S1^2]) + (9765*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i])/(16*a^5*S1*Sqrt[1 - (S1 - S3)^2/S1^2]) - (315*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i])/(16*a^5*S1*Sqrt[1 - (S1 - S3)^2/S1^2]) - (235305*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i])/(32*a^5*S1*Sqrt[1 - (S1 - S3)^2/S1^2]) + (315*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i])/(2*a^5*S1^2) - (945*e*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i])/(a^5*S1^2) - (9765*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i])/(8*a^5*S1^2) + (315*e^3*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i])/(8*a^5*S1^2) + (235305*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i])/(16*a^5*S1^2) - (2835*e*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]^2*Sin[i])/(4*a^5) + (945*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]^2*Sin[i])/(32*a^5) - (945*e*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^2*Sin[i])/(a^5*S1^2) + (16065*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^2*Sin[i])/(4*a^5*S1^2) + (315*e^3*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^2*Sin[i])/(8*a^5*S1^2) - (101115*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^2*Sin[i])/(4*a^5*S1^2) - (945*e*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^2*Sin[i])/(4*a^5*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) + (16065*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^2*Sin[i])/(16*a^5*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) + (315*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^2*Sin[i])/(32*a^5*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) - (101115*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^2*Sin[i])/(16*a^5*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) + (945*e*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^2*Sin[i])/(a^5*S1^3) - (16065*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^2*Sin[i])/(4*a^5*S1^3) - (315*e^3*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^2*Sin[i])/(8*a^5*S1^3) + (101115*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^2*Sin[i])/(4*a^5*S1^3) + (315*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(2*a^5*S1) - (945*e*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(a^5*S1) + (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(2*a^5*S1) + (315*e^3*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(8*a^5*S1) + (33075*e^4*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(16*a^5*S1) + (945*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^3*Sin[i])/(4*a^5*S1) - (2835*e*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^3*Sin[i])/(2*a^5*S1) + (4725*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^3*Sin[i])/(4*a^5*S1) + (945*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^3*Sin[i])/(16*a^5*S1) + (99225*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^3*Sin[i])/(32*a^5*S1) - (315*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*(S1 - S3)*Cos[i]^3*Sin[i])/(2*a^5*S1^2) + (945*e*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*(S1 - S3)*Cos[i]^3*Sin[i])/(a^5*S1^2) - (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*(S1 - S3)*Cos[i]^3*Sin[i])/(2*a^5*S1^2) - (315*e^3*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*(S1 - S3)*Cos[i]^3*Sin[i])/(8*a^5*S1^2) - (33075*e^4*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*(S1 - S3)*Cos[i]^3*Sin[i])/(16*a^5*S1^2) + (945*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^3*Sin[i])/(2*a^5*S1^3) - (2835*e*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^3*Sin[i])/(a^5*S1^3) + (34965*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^3*Sin[i])/(8*a^5*S1^3) + (945*e^3*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^3*Sin[i])/(8*a^5*S1^3) - (103005*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^3*Sin[i])/(16*a^5*S1^3) + (315*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^3*Cos[i]^3*Sin[i])/(4*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) - (945*e*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^3*Cos[i]^3*Sin[i])/(2*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) + (11655*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^3*Cos[i]^3*Sin[i])/(16*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) + (315*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^3*Cos[i]^3*Sin[i])/(16*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) - (34335*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^3*Cos[i]^3*Sin[i])/(32*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) - (945*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^3*Cos[i]^3*Sin[i])/(2*a^5*S1^4) + (2835*e*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^3*Cos[i]^3*Sin[i])/(a^5*S1^4) - (34965*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^3*Cos[i]^3*Sin[i])/(8*a^5*S1^4) - (945*e^3*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^3*Cos[i]^3*Sin[i])/(8*a^5*S1^4) + (103005*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^3*Cos[i]^3*Sin[i])/(16*a^5*S1^4) - (525*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i]^2)/(4*a^5) + (945*e*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i]^2)/(2*a^5) + (5565*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i]^2)/(16*a^5) - (315*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i]^2)/(16*a^5) - (257355*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i]^2)/(32*a^5) - (105*R^4*(S1 - S3)*Sin[i]^2)/(a^5*S1^2) - (525*e^2*R^4*(S1 - S3)*Sin[i]^2)/(a^5*S1^2) - (11025*e^4*R^4*(S1 - S3)*Sin[i]^2)/(8*a^5*S1^2) + (105*R^4*(S1 - S3)^2*Sin[i]^2)/(a^5*S1^3) + (525*e^2*R^4*(S1 - S3)^2*Sin[i]^2)/(a^5*S1^3) + (11025*e^4*R^4*(S1 - S3)^2*Sin[i]^2)/(8*a^5*S1^3) - (945*e*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]*Sin[i]^2)/(2*a^5*S1) + (16065*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]*Sin[i]^2)/(8*a^5*S1) + (315*e^3*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]*Sin[i]^2)/(16*a^5*S1) - (101115*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]*Sin[i]^2)/(8*a^5*S1) - (945*e*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i]^2)/(2*a^5*S1) + (16065*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i]^2)/(8*a^5*S1) + (315*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i]^2)/(16*a^5*S1) - (101115*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i]^2)/(8*a^5*S1) + (945*e*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]*Sin[i]^2)/(2*a^5*S1^2) - (16065*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]*Sin[i]^2)/(8*a^5*S1^2) - (315*e^3*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]*Sin[i]^2)/(16*a^5*S1^2) + (101115*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]*Sin[i]^2)/(8*a^5*S1^2) - (2835*e*R^4*(S1 - S3)^2*Cos[i]*Sin[i]^2)/(2*a^5*S1^3) + (945*e^3*R^4*(S1 - S3)^2*Cos[i]*Sin[i]^2)/(16*a^5*S1^3) + (2835*e*R^4*(S1 - S3)^3*Cos[i]*Sin[i]^2)/(2*a^5*S1^4) - (945*e^3*R^4*(S1 - S3)^3*Cos[i]*Sin[i]^2)/(16*a^5*S1^4) + (315*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(1 - (S1 - S3)^2/S1^2)*Cos[i]^2*Sin[i]^2)/a^5 - (945*e*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(1 - (S1 - S3)^2/S1^2)*Cos[i]^2*Sin[i]^2)/a^5 + (1575*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(1 - (S1 - S3)^2/S1^2)*Cos[i]^2*Sin[i]^2)/a^5 + (315*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(1 - (S1 - S3)^2/S1^2)*Cos[i]^2*Sin[i]^2)/(8*a^5) + (33075*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(1 - (S1 - S3)^2/S1^2)*Cos[i]^2*Sin[i]^2)/(8*a^5) + (630*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^2*Sin[i]^2)/(a^5*S1^2) - (1890*e*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^2*Sin[i]^2)/(a^5*S1^2) + (41265*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^2*Sin[i]^2)/(8*a^5*S1^2) + (315*e^3*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^2*Sin[i]^2)/(4*a^5*S1^2) - (34965*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^2*Sin[i]^2)/(8*a^5*S1^2) + (315*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^2*Sin[i]^2)/(a^5*S1^2) - (945*e*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^2*Sin[i]^2)/(a^5*S1^2) + (41265*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^2*Sin[i]^2)/(16*a^5*S1^2) + (315*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^2*Sin[i]^2)/(8*a^5*S1^2) - (34965*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^2*Sin[i]^2)/(16*a^5*S1^2) - (630*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Cos[i]^2*Sin[i]^2)/(a^5*S1^3) + (1890*e*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Cos[i]^2*Sin[i]^2)/(a^5*S1^3) - (41265*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Cos[i]^2*Sin[i]^2)/(8*a^5*S1^3) - (315*e^3*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Cos[i]^2*Sin[i]^2)/(4*a^5*S1^3) + (34965*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Cos[i]^2*Sin[i]^2)/(8*a^5*S1^3) + (630*R^4*(S1 - S3)^3*Cos[i]^2*Sin[i]^2)/(a^5*S1^4) - (1890*e*R^4*(S1 - S3)^3*Cos[i]^2*Sin[i]^2)/(a^5*S1^4) + (3150*e^2*R^4*(S1 - S3)^3*Cos[i]^2*Sin[i]^2)/(a^5*S1^4) + (315*e^3*R^4*(S1 - S3)^3*Cos[i]^2*Sin[i]^2)/(4*a^5*S1^4) + (33075*e^4*R^4*(S1 - S3)^3*Cos[i]^2*Sin[i]^2)/(4*a^5*S1^4) - (630*R^4*(S1 - S3)^4*Cos[i]^2*Sin[i]^2)/(a^5*S1^5) + (1890*e*R^4*(S1 - S3)^4*Cos[i]^2*Sin[i]^2)/(a^5*S1^5) - (3150*e^2*R^4*(S1 - S3)^4*Cos[i]^2*Sin[i]^2)/(a^5*S1^5) - (315*e^3*R^4*(S1 - S3)^4*Cos[i]^2*Sin[i]^2)/(4*a^5*S1^5) - (33075*e^4*R^4*(S1 - S3)^4*Cos[i]^2*Sin[i]^2)/(4*a^5*S1^5) - (2835*e*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sqrt[1 - (S1 - S3)^2/S1^2]*Sin[i]^3)/(4*a^5) + (16065*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sqrt[1 - (S1 - S3)^2/S1^2]*Sin[i]^3)/(16*a^5) + (945*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sqrt[1 - (S1 - S3)^2/S1^2]*Sin[i]^3)/(32*a^5) - (101115*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sqrt[1 - (S1 - S3)^2/S1^2]*Sin[i]^3)/(16*a^5) - (945*e*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Sin[i]^3)/(a^5*S1^2) + (315*e^3*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Sin[i]^3)/(8*a^5*S1^2) - (945*e*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Sin[i]^3)/(4*a^5*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) + (315*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Sin[i]^3)/(32*a^5*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) + (945*e*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Sin[i]^3)/(a^5*S1^3) - (315*e^3*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Sin[i]^3)/(8*a^5*S1^3) + (315*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]*Sin[i]^3)/(2*a^5*S1) - (945*e*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]*Sin[i]^3)/(a^5*S1) + (11655*e^2*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]*Sin[i]^3)/(8*a^5*S1) + (315*e^3*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]*Sin[i]^3)/(8*a^5*S1) - (34335*e^4*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]*Sin[i]^3)/(16*a^5*S1) + (945*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i]^3)/(4*a^5*S1) - (2835*e*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i]^3)/(2*a^5*S1) + (34965*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i]^3)/(16*a^5*S1) + (945*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i]^3)/(16*a^5*S1) - (103005*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i]^3)/(32*a^5*S1) - (315*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*(S1 - S3)*Cos[i]*Sin[i]^3)/(2*a^5*S1^2) + (945*e*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*(S1 - S3)*Cos[i]*Sin[i]^3)/(a^5*S1^2) - (11655*e^2*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*(S1 - S3)*Cos[i]*Sin[i]^3)/(8*a^5*S1^2) - (315*e^3*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*(S1 - S3)*Cos[i]*Sin[i]^3)/(8*a^5*S1^2) + (34335*e^4*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*(S1 - S3)*Cos[i]*Sin[i]^3)/(16*a^5*S1^2) + (945*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]*Sin[i]^3)/(2*a^5*S1^3) - (2835*e*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]*Sin[i]^3)/(a^5*S1^3) + (4725*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]*Sin[i]^3)/(2*a^5*S1^3) + (945*e^3*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]*Sin[i]^3)/(8*a^5*S1^3) + (99225*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]*Sin[i]^3)/(16*a^5*S1^3) + (315*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^3*Cos[i]*Sin[i]^3)/(4*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) - (945*e*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^3*Cos[i]*Sin[i]^3)/(2*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) + (1575*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^3*Cos[i]*Sin[i]^3)/(4*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) + (315*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^3*Cos[i]*Sin[i]^3)/(16*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) + (33075*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^3*Cos[i]*Sin[i]^3)/(32*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) - (945*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^3*Cos[i]*Sin[i]^3)/(2*a^5*S1^4) + (2835*e*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^3*Cos[i]*Sin[i]^3)/(a^5*S1^4) - (4725*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^3*Cos[i]*Sin[i]^3)/(2*a^5*S1^4) - (945*e^3*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^3*Cos[i]*Sin[i]^3)/(8*a^5*S1^4) - (99225*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^3*Cos[i]*Sin[i]^3)/(16*a^5*S1^4) + (945*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(1 - (S1 - S3)^2/S1^2)*Sin[i]^4)/(4*a^5) - (945*e*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(1 - (S1 - S3)^2/S1^2)*Sin[i]^4)/a^5 + (24255*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(1 - (S1 - S3)^2/S1^2)*Sin[i]^4)/(16*a^5) + (315*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(1 - (S1 - S3)^2/S1^2)*Sin[i]^4)/(8*a^5) + (31815*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(1 - (S1 - S3)^2/S1^2)*Sin[i]^4)/(32*a^5) + (315*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Sin[i]^4)/(a^5*S1^2) - (945*e*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Sin[i]^4)/(a^5*S1^2) + (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Sin[i]^4)/(a^5*S1^2) + (315*e^3*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Sin[i]^4)/(8*a^5*S1^2) + (33075*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Sin[i]^4)/(8*a^5*S1^2) + (315*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Sin[i]^4)/(2*a^5*S1^2) - (945*e*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Sin[i]^4)/(2*a^5*S1^2) + (1575*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Sin[i]^4)/(2*a^5*S1^2) + (315*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Sin[i]^4)/(16*a^5*S1^2) + (33075*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Sin[i]^4)/(16*a^5*S1^2) - (315*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Sin[i]^4)/(a^5*S1^3) + (945*e*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Sin[i]^4)/(a^5*S1^3) - (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Sin[i]^4)/(a^5*S1^3) - (315*e^3*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Sin[i]^4)/(8*a^5*S1^3) - (33075*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Sin[i]^4)/(8*a^5*S1^3) + (315*R^4*(S1 - S3)^3*Sin[i]^4)/(2*a^5*S1^4) + (1575*e^2*R^4*(S1 - S3)^3*Sin[i]^4)/(2*a^5*S1^4) + (33075*e^4*R^4*(S1 - S3)^3*Sin[i]^4)/(16*a^5*S1^4) - (315*R^4*(S1 - S3)^4*Sin[i]^4)/(2*a^5*S1^5) - (1575*e^2*R^4*(S1 - S3)^4*Sin[i]^4)/(2*a^5*S1^5) - (33075*e^4*R^4*(S1 - S3)^4*Sin[i]^4)/(16*a^5*S1^5)) + C40*((-15*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Cos[i]^2)/(8*a^5) - (75*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Cos[i]^2)/(8*a^5) - (1575*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Cos[i]^2)/(64*a^5) + (105*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(1 - (S1 - S3)^2/S1^2)*Cos[i]^4)/(32*a^5) + (525*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(1 - (S1 - S3)^2/S1^2)*Cos[i]^4)/(32*a^5) + (11025*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(1 - (S1 - S3)^2/S1^2)*Cos[i]^4)/(256*a^5) + (15*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(4*a^5*S1) + (75*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(4*a^5*S1) + (1575*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(32*a^5*S1) + (15*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i])/(8*a^5*S1*Sqrt[1 - (S1 - S3)^2/S1^2]) + (75*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i])/(8*a^5*S1*Sqrt[1 - (S1 - S3)^2/S1^2]) + (1575*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i])/(64*a^5*S1*Sqrt[1 - (S1 - S3)^2/S1^2]) - (15*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i])/(4*a^5*S1^2) - (75*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i])/(4*a^5*S1^2) - (1575*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i])/(32*a^5*S1^2) - (105*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(16*a^5*S1) - (525*e^2*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(16*a^5*S1) - (11025*e^4*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(128*a^5*S1) - (315*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^3*Sin[i])/(32*a^5*S1) - (1575*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^3*Sin[i])/(32*a^5*S1) - (33075*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^3*Sin[i])/(256*a^5*S1) + (105*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*(S1 - S3)*Cos[i]^3*Sin[i])/(16*a^5*S1^2) + (525*e^2*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*(S1 - S3)*Cos[i]^3*Sin[i])/(16*a^5*S1^2) + (11025*e^4*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*(S1 - S3)*Cos[i]^3*Sin[i])/(128*a^5*S1^2) - (15*R^4*(S1 - S3)*Sin[i]^2)/(4*a^5*S1^2) - (75*e^2*R^4*(S1 - S3)*Sin[i]^2)/(4*a^5*S1^2) - (1575*e^4*R^4*(S1 - S3)*Sin[i]^2)/(32*a^5*S1^2) + (15*R^4*(S1 - S3)^2*Sin[i]^2)/(4*a^5*S1^3) + (75*e^2*R^4*(S1 - S3)^2*Sin[i]^2)/(4*a^5*S1^3) + (1575*e^4*R^4*(S1 - S3)^2*Sin[i]^2)/(32*a^5*S1^3) + (315*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^2*Sin[i]^2)/(16*a^5*S1^2) + (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^2*Sin[i]^2)/(16*a^5*S1^2) + (33075*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^2*Sin[i]^2)/(128*a^5*S1^2) + (315*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^2*Sin[i]^2)/(32*a^5*S1^2) + (1575*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^2*Sin[i]^2)/(32*a^5*S1^2) + (33075*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^2*Sin[i]^2)/(256*a^5*S1^2) - (315*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Cos[i]^2*Sin[i]^2)/(16*a^5*S1^3) - (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Cos[i]^2*Sin[i]^2)/(16*a^5*S1^3) - (33075*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Cos[i]^2*Sin[i]^2)/(128*a^5*S1^3) - (315*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]*Sin[i]^3)/(16*a^5*S1^3) - (1575*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]*Sin[i]^3)/(16*a^5*S1^3) - (33075*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]*Sin[i]^3)/(128*a^5*S1^3) - (105*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^3*Cos[i]*Sin[i]^3)/(32*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) - (525*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^3*Cos[i]*Sin[i]^3)/(32*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) - (11025*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^3*Cos[i]*Sin[i]^3)/(256*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) + (315*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^3*Cos[i]*Sin[i]^3)/(16*a^5*S1^4) + (1575*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^3*Cos[i]*Sin[i]^3)/(16*a^5*S1^4) + (33075*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^3*Cos[i]*Sin[i]^3)/(128*a^5*S1^4) + (105*R^4*(S1 - S3)^3*Sin[i]^4)/(16*a^5*S1^4) + (525*e^2*R^4*(S1 - S3)^3*Sin[i]^4)/(16*a^5*S1^4) + (11025*e^4*R^4*(S1 - S3)^3*Sin[i]^4)/(128*a^5*S1^4) - (105*R^4*(S1 - S3)^4*Sin[i]^4)/(16*a^5*S1^5) - (525*e^2*R^4*(S1 - S3)^4*Sin[i]^4)/(16*a^5*S1^5) - (11025*e^4*R^4*(S1 - S3)^4*Sin[i]^4)/(128*a^5*S1^5)) + C42*((-135*e*R^4*Cos[i])/(8*a^5*S1) + (45*e^3*R^4*Cos[i])/(64*a^5*S1) + (135*e*R^4*(S1 - S3)*Cos[i])/(8*a^5*S1^2) - (45*e^3*R^4*(S1 - S3)*Cos[i])/(64*a^5*S1^2) + (375*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Cos[i]^2)/(16*a^5) + (1875*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Cos[i]^2)/(16*a^5) + (39375*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Cos[i]^2)/(128*a^5) + (15*R^4*(S1 - S3)*Cos[i]^2)/(2*a^5*S1^2) - (135*e*R^4*(S1 - S3)*Cos[i]^2)/(8*a^5*S1^2) + (75*e^2*R^4*(S1 - S3)*Cos[i]^2)/(2*a^5*S1^2) + (45*e^3*R^4*(S1 - S3)*Cos[i]^2)/(64*a^5*S1^2) + (1575*e^4*R^4*(S1 - S3)*Cos[i]^2)/(16*a^5*S1^2) - (15*R^4*(S1 - S3)^2*Cos[i]^2)/(2*a^5*S1^3) + (135*e*R^4*(S1 - S3)^2*Cos[i]^2)/(8*a^5*S1^3) - (75*e^2*R^4*(S1 - S3)^2*Cos[i]^2)/(2*a^5*S1^3) - (45*e^3*R^4*(S1 - S3)^2*Cos[i]^2)/(64*a^5*S1^3) - (1575*e^4*R^4*(S1 - S3)^2*Cos[i]^2)/(16*a^5*S1^3) + (945*e*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]^3)/(16*a^5*S1) - (315*e^3*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]^3)/(128*a^5*S1) + (945*e*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]^3)/(16*a^5*S1) - (315*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]^3)/(128*a^5*S1) - (945*e*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^3)/(16*a^5*S1^2) + (315*e^3*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^3)/(128*a^5*S1^2) - (315*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(1 - (S1 - S3)^2/S1^2)*Cos[i]^4)/(8*a^5) - (1575*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(1 - (S1 - S3)^2/S1^2)*Cos[i]^4)/(8*a^5) - (33075*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(1 - (S1 - S3)^2/S1^2)*Cos[i]^4)/(64*a^5) - (315*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^4)/(8*a^5*S1^2) + (945*e*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^4)/(8*a^5*S1^2) - (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^4)/(8*a^5*S1^2) - (315*e^3*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^4)/(64*a^5*S1^2) - (33075*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^4)/(64*a^5*S1^2) - (315*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^4)/(16*a^5*S1^2) + (945*e*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^4)/(16*a^5*S1^2) - (1575*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^4)/(16*a^5*S1^2) - (315*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^4)/(128*a^5*S1^2) - (33075*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^4)/(128*a^5*S1^2) + (315*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Cos[i]^4)/(8*a^5*S1^3) - (945*e*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Cos[i]^4)/(8*a^5*S1^3) + (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Cos[i]^4)/(8*a^5*S1^3) + (315*e^3*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Cos[i]^4)/(64*a^5*S1^3) + (33075*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Cos[i]^4)/(64*a^5*S1^3) - (135*e*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i])/(16*a^5*Sqrt[1 - (S1 - S3)^2/S1^2]) + (45*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i])/(128*a^5*Sqrt[1 - (S1 - S3)^2/S1^2]) - (315*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(8*a^5*S1) - (135*e*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(8*a^5*S1) - (1575*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(8*a^5*S1) + (45*e^3*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(64*a^5*S1) - (33075*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(64*a^5*S1) - (315*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i])/(16*a^5*S1*Sqrt[1 - (S1 - S3)^2/S1^2]) - (135*e*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i])/(16*a^5*S1*Sqrt[1 - (S1 - S3)^2/S1^2]) - (1575*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i])/(16*a^5*S1*Sqrt[1 - (S1 - S3)^2/S1^2]) + (45*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i])/(128*a^5*S1*Sqrt[1 - (S1 - S3)^2/S1^2]) - (33075*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i])/(128*a^5*S1*Sqrt[1 - (S1 - S3)^2/S1^2]) + (315*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i])/(8*a^5*S1^2) + (135*e*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i])/(8*a^5*S1^2) + (1575*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i])/(8*a^5*S1^2) - (45*e^3*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i])/(64*a^5*S1^2) + (33075*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i])/(64*a^5*S1^2) + (2835*e*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]^2*Sin[i])/(32*a^5) - (945*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]^2*Sin[i])/(256*a^5) - (945*e*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^2*Sin[i])/(4*a^5*S1^2) + (315*e^3*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^2*Sin[i])/(32*a^5*S1^2) - (945*e*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^2*Sin[i])/(16*a^5*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) + (315*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^2*Sin[i])/(128*a^5*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) + (945*e*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^2*Sin[i])/(4*a^5*S1^3) - (315*e^3*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^2*Sin[i])/(32*a^5*S1^3) + (315*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(8*a^5*S1) + (945*e*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(8*a^5*S1) + (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(8*a^5*S1) - (315*e^3*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(64*a^5*S1) + (33075*e^4*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(64*a^5*S1) + (945*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^3*Sin[i])/(16*a^5*S1) + (2835*e*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^3*Sin[i])/(16*a^5*S1) + (4725*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^3*Sin[i])/(16*a^5*S1) - (945*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^3*Sin[i])/(128*a^5*S1) + (99225*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^3*Sin[i])/(128*a^5*S1) - (315*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*(S1 - S3)*Cos[i]^3*Sin[i])/(8*a^5*S1^2) - (945*e*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*(S1 - S3)*Cos[i]^3*Sin[i])/(8*a^5*S1^2) - (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*(S1 - S3)*Cos[i]^3*Sin[i])/(8*a^5*S1^2) + (315*e^3*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*(S1 - S3)*Cos[i]^3*Sin[i])/(64*a^5*S1^2) - (33075*e^4*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*(S1 - S3)*Cos[i]^3*Sin[i])/(64*a^5*S1^2) + (945*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^3*Sin[i])/(8*a^5*S1^3) - (2835*e*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^3*Sin[i])/(8*a^5*S1^3) + (4725*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^3*Sin[i])/(8*a^5*S1^3) + (945*e^3*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^3*Sin[i])/(64*a^5*S1^3) + (99225*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^3*Sin[i])/(64*a^5*S1^3) + (315*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^3*Cos[i]^3*Sin[i])/(16*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) - (945*e*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^3*Cos[i]^3*Sin[i])/(16*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) + (1575*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^3*Cos[i]^3*Sin[i])/(16*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) + (315*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^3*Cos[i]^3*Sin[i])/(128*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) + (33075*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^3*Cos[i]^3*Sin[i])/(128*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) - (945*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^3*Cos[i]^3*Sin[i])/(8*a^5*S1^4) + (2835*e*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^3*Cos[i]^3*Sin[i])/(8*a^5*S1^4) - (4725*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^3*Cos[i]^3*Sin[i])/(8*a^5*S1^4) - (945*e^3*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^3*Cos[i]^3*Sin[i])/(64*a^5*S1^4) - (99225*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^3*Cos[i]^3*Sin[i])/(64*a^5*S1^4) + (15*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i]^2)/(4*a^5) - (135*e*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i]^2)/(16*a^5) + (75*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i]^2)/(4*a^5) + (45*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i]^2)/(128*a^5) + (1575*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i]^2)/(32*a^5) + (375*R^4*(S1 - S3)*Sin[i]^2)/(8*a^5*S1^2) + (1875*e^2*R^4*(S1 - S3)*Sin[i]^2)/(8*a^5*S1^2) + (39375*e^4*R^4*(S1 - S3)*Sin[i]^2)/(64*a^5*S1^2) - (375*R^4*(S1 - S3)^2*Sin[i]^2)/(8*a^5*S1^3) - (1875*e^2*R^4*(S1 - S3)^2*Sin[i]^2)/(8*a^5*S1^3) - (39375*e^4*R^4*(S1 - S3)^2*Sin[i]^2)/(64*a^5*S1^3) - (945*e*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]*Sin[i]^2)/(8*a^5*S1) + (315*e^3*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]*Sin[i]^2)/(64*a^5*S1) - (945*e*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i]^2)/(8*a^5*S1) + (315*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i]^2)/(64*a^5*S1) + (945*e*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]*Sin[i]^2)/(8*a^5*S1^2) - (315*e^3*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]*Sin[i]^2)/(64*a^5*S1^2) + (2835*e*R^4*(S1 - S3)^2*Cos[i]*Sin[i]^2)/(16*a^5*S1^3) - (945*e^3*R^4*(S1 - S3)^2*Cos[i]*Sin[i]^2)/(128*a^5*S1^3) - (2835*e*R^4*(S1 - S3)^3*Cos[i]*Sin[i]^2)/(16*a^5*S1^4) + (945*e^3*R^4*(S1 - S3)^3*Cos[i]*Sin[i]^2)/(128*a^5*S1^4) - (315*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(1 - (S1 - S3)^2/S1^2)*Cos[i]^2*Sin[i]^2)/(8*a^5) + (945*e*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(1 - (S1 - S3)^2/S1^2)*Cos[i]^2*Sin[i]^2)/(8*a^5) - (1575*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(1 - (S1 - S3)^2/S1^2)*Cos[i]^2*Sin[i]^2)/(8*a^5) - (315*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(1 - (S1 - S3)^2/S1^2)*Cos[i]^2*Sin[i]^2)/(64*a^5) - (33075*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(1 - (S1 - S3)^2/S1^2)*Cos[i]^2*Sin[i]^2)/(64*a^5) - (315*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^2*Sin[i]^2)/(4*a^5*S1^2) - (945*e*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^2*Sin[i]^2)/(2*a^5*S1^2) - (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^2*Sin[i]^2)/(4*a^5*S1^2) + (315*e^3*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^2*Sin[i]^2)/(16*a^5*S1^2) - (33075*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^2*Sin[i]^2)/(32*a^5*S1^2) - (315*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^2*Sin[i]^2)/(8*a^5*S1^2) - (945*e*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^2*Sin[i]^2)/(4*a^5*S1^2) - (1575*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^2*Sin[i]^2)/(8*a^5*S1^2) + (315*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^2*Sin[i]^2)/(32*a^5*S1^2) - (33075*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^2*Sin[i]^2)/(64*a^5*S1^2) + (315*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Cos[i]^2*Sin[i]^2)/(4*a^5*S1^3) + (945*e*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Cos[i]^2*Sin[i]^2)/(2*a^5*S1^3) + (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Cos[i]^2*Sin[i]^2)/(4*a^5*S1^3) - (315*e^3*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Cos[i]^2*Sin[i]^2)/(16*a^5*S1^3) + (33075*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Cos[i]^2*Sin[i]^2)/(32*a^5*S1^3) - (315*R^4*(S1 - S3)^3*Cos[i]^2*Sin[i]^2)/(4*a^5*S1^4) + (945*e*R^4*(S1 - S3)^3*Cos[i]^2*Sin[i]^2)/(4*a^5*S1^4) - (1575*e^2*R^4*(S1 - S3)^3*Cos[i]^2*Sin[i]^2)/(4*a^5*S1^4) - (315*e^3*R^4*(S1 - S3)^3*Cos[i]^2*Sin[i]^2)/(32*a^5*S1^4) - (33075*e^4*R^4*(S1 - S3)^3*Cos[i]^2*Sin[i]^2)/(32*a^5*S1^4) + (315*R^4*(S1 - S3)^4*Cos[i]^2*Sin[i]^2)/(4*a^5*S1^5) - (945*e*R^4*(S1 - S3)^4*Cos[i]^2*Sin[i]^2)/(4*a^5*S1^5) + (1575*e^2*R^4*(S1 - S3)^4*Cos[i]^2*Sin[i]^2)/(4*a^5*S1^5) + (315*e^3*R^4*(S1 - S3)^4*Cos[i]^2*Sin[i]^2)/(32*a^5*S1^5) + (33075*e^4*R^4*(S1 - S3)^4*Cos[i]^2*Sin[i]^2)/(32*a^5*S1^5) + (945*e*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Sin[i]^3)/(8*a^5*S1^2) - (315*e^3*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Sin[i]^3)/(64*a^5*S1^2) + (945*e*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Sin[i]^3)/(32*a^5*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) - (315*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Sin[i]^3)/(256*a^5*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) - (945*e*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Sin[i]^3)/(8*a^5*S1^3) + (315*e^3*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Sin[i]^3)/(64*a^5*S1^3) + (315*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]*Sin[i]^3)/(8*a^5*S1) - (945*e*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]*Sin[i]^3)/(8*a^5*S1) + (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]*Sin[i]^3)/(8*a^5*S1) + (315*e^3*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]*Sin[i]^3)/(64*a^5*S1) + (33075*e^4*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]*Sin[i]^3)/(64*a^5*S1) + (945*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i]^3)/(16*a^5*S1) - (2835*e*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i]^3)/(16*a^5*S1) + (4725*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i]^3)/(16*a^5*S1) + (945*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i]^3)/(128*a^5*S1) + (99225*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i]^3)/(128*a^5*S1) - (315*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*(S1 - S3)*Cos[i]*Sin[i]^3)/(8*a^5*S1^2) + (945*e*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*(S1 - S3)*Cos[i]*Sin[i]^3)/(8*a^5*S1^2) - (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*(S1 - S3)*Cos[i]*Sin[i]^3)/(8*a^5*S1^2) - (315*e^3*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*(S1 - S3)*Cos[i]*Sin[i]^3)/(64*a^5*S1^2) - (33075*e^4*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*(S1 - S3)*Cos[i]*Sin[i]^3)/(64*a^5*S1^2) + (945*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]*Sin[i]^3)/(8*a^5*S1^3) + (2835*e*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]*Sin[i]^3)/(8*a^5*S1^3) + (4725*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]*Sin[i]^3)/(8*a^5*S1^3) - (945*e^3*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]*Sin[i]^3)/(64*a^5*S1^3) + (99225*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]*Sin[i]^3)/(64*a^5*S1^3) + (315*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^3*Cos[i]*Sin[i]^3)/(16*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) + (945*e*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^3*Cos[i]*Sin[i]^3)/(16*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) + (1575*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^3*Cos[i]*Sin[i]^3)/(16*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) - (315*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^3*Cos[i]*Sin[i]^3)/(128*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) + (33075*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^3*Cos[i]*Sin[i]^3)/(128*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) - (945*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^3*Cos[i]*Sin[i]^3)/(8*a^5*S1^4) - (2835*e*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^3*Cos[i]*Sin[i]^3)/(8*a^5*S1^4) - (4725*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^3*Cos[i]*Sin[i]^3)/(8*a^5*S1^4) + (945*e^3*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^3*Cos[i]*Sin[i]^3)/(64*a^5*S1^4) - (99225*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^3*Cos[i]*Sin[i]^3)/(64*a^5*S1^4) - (315*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Sin[i]^4)/(8*a^5*S1^2) + (945*e*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Sin[i]^4)/(8*a^5*S1^2) - (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Sin[i]^4)/(8*a^5*S1^2) - (315*e^3*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Sin[i]^4)/(64*a^5*S1^2) - (33075*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Sin[i]^4)/(64*a^5*S1^2) - (315*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Sin[i]^4)/(16*a^5*S1^2) + (945*e*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Sin[i]^4)/(16*a^5*S1^2) - (1575*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Sin[i]^4)/(16*a^5*S1^2) - (315*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Sin[i]^4)/(128*a^5*S1^2) - (33075*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Sin[i]^4)/(128*a^5*S1^2) + (315*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Sin[i]^4)/(8*a^5*S1^3) - (945*e*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Sin[i]^4)/(8*a^5*S1^3) + (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Sin[i]^4)/(8*a^5*S1^3) + (315*e^3*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Sin[i]^4)/(64*a^5*S1^3) + (33075*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Sin[i]^4)/(64*a^5*S1^3) - (315*R^4*(S1 - S3)^3*Sin[i]^4)/(4*a^5*S1^4) - (1575*e^2*R^4*(S1 - S3)^3*Sin[i]^4)/(4*a^5*S1^4) - (33075*e^4*R^4*(S1 - S3)^3*Sin[i]^4)/(32*a^5*S1^4) + (315*R^4*(S1 - S3)^4*Sin[i]^4)/(4*a^5*S1^5) + (1575*e^2*R^4*(S1 - S3)^4*Sin[i]^4)/(4*a^5*S1^5) + (33075*e^4*R^4*(S1 - S3)^4*Sin[i]^4)/(32*a^5*S1^5))) == 0, Omd - G*m*M*(C22*((-21*e*R^2*Cos[i])/(4*a^3*S1) + (369*e^3*R^2*Cos[i])/(32*a^3*S1) - (21*e*R^2*(S1 - S3)*Cos[i]^2)/(4*a^3*S1^2) + (369*e^3*R^2*(S1 - S3)*Cos[i]^2)/(32*a^3*S1^2) + (21*e*R^2*(S1 - S3)*Sin[i])/(4*a^3*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) - (369*e^3*R^2*(S1 - S3)*Sin[i])/(32*a^3*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) - (21*e*R^2*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(4*a^3*S1) + (369*e^3*R^2*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(32*a^3*S1) + (21*e*R^2*(S1 - S3)^2*Cos[i]*Sin[i])/(4*a^3*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) - (369*e^3*R^2*(S1 - S3)^2*Cos[i]*Sin[i])/(32*a^3*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) + (21*e*R^2*(S1 - S3)*Sin[i]^2)/(4*a^3*S1^2) - (369*e^3*R^2*(S1 - S3)*Sin[i]^2)/(32*a^3*S1^2)) + C20*((3*R^2*(S1 - S3)*Cos[i]^2)/(2*a^3*S1^2) + (9*e^2*R^2*(S1 - S3)*Cos[i]^2)/(4*a^3*S1^2) + (45*e^4*R^2*(S1 - S3)*Cos[i]^2)/(16*a^3*S1^2) + (3*R^2*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(2*a^3*S1) + (9*e^2*R^2*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(4*a^3*S1) + (45*e^4*R^2*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(16*a^3*S1) - (3*R^2*(S1 - S3)^2*Cos[i]*Sin[i])/(2*a^3*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) - (9*e^2*R^2*(S1 - S3)^2*Cos[i]*Sin[i])/(4*a^3*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) - (45*e^4*R^2*(S1 - S3)^2*Cos[i]*Sin[i])/(16*a^3*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) - (3*R^2*(S1 - S3)*Sin[i]^2)/(2*a^3*S1^2) - (9*e^2*R^2*(S1 - S3)*Sin[i]^2)/(4*a^3*S1^2) - (45*e^4*R^2*(S1 - S3)*Sin[i]^2)/(16*a^3*S1^2)) + C40*((-15*R^4*(S1 - S3)*Cos[i]^2)/(4*a^5*S1^2) - (75*e^2*R^4*(S1 - S3)*Cos[i]^2)/(4*a^5*S1^2) - (1575*e^4*R^4*(S1 - S3)*Cos[i]^2)/(32*a^5*S1^2) + (105*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^4)/(16*a^5*S1^2) + (525*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^4)/(16*a^5*S1^2) + (11025*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^4)/(128*a^5*S1^2) - (15*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(4*a^5*S1) - (75*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(4*a^5*S1) - (1575*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(32*a^5*S1) + (15*R^4*(S1 - S3)^2*Cos[i]*Sin[i])/(4*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) + (75*e^2*R^4*(S1 - S3)^2*Cos[i]*Sin[i])/(4*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) + (1575*e^4*R^4*(S1 - S3)^2*Cos[i]*Sin[i])/(32*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) + (105*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(16*a^5*S1) + (525*e^2*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(16*a^5*S1) + (11025*e^4*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(128*a^5*S1) - (315*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^3*Sin[i])/(16*a^5*S1^3) - (1575*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^3*Sin[i])/(16*a^5*S1^3) - (33075*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^3*Sin[i])/(128*a^5*S1^3) + (15*R^4*(S1 - S3)*Sin[i]^2)/(4*a^5*S1^2) + (75*e^2*R^4*(S1 - S3)*Sin[i]^2)/(4*a^5*S1^2) + (1575*e^4*R^4*(S1 - S3)*Sin[i]^2)/(32*a^5*S1^2) - (315*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^2*Sin[i]^2)/(16*a^5*S1^2) - (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^2*Sin[i]^2)/(16*a^5*S1^2) - (33075*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^2*Sin[i]^2)/(128*a^5*S1^2) + (315*R^4*(S1 - S3)^3*Cos[i]^2*Sin[i]^2)/(16*a^5*S1^4) + (1575*e^2*R^4*(S1 - S3)^3*Cos[i]^2*Sin[i]^2)/(16*a^5*S1^4) + (33075*e^4*R^4*(S1 - S3)^3*Cos[i]^2*Sin[i]^2)/(128*a^5*S1^4) + (315*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]*Sin[i]^3)/(16*a^5*S1^3) + (1575*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]*Sin[i]^3)/(16*a^5*S1^3) + (33075*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]*Sin[i]^3)/(128*a^5*S1^3) - (105*R^4*(S1 - S3)^4*Cos[i]*Sin[i]^3)/(16*a^5*S1^5*Sqrt[1 - (S1 - S3)^2/S1^2]) - (525*e^2*R^4*(S1 - S3)^4*Cos[i]*Sin[i]^3)/(16*a^5*S1^5*Sqrt[1 - (S1 - S3)^2/S1^2]) - (11025*e^4*R^4*(S1 - S3)^4*Cos[i]*Sin[i]^3)/(128*a^5*S1^5*Sqrt[1 - (S1 - S3)^2/S1^2]) - (105*R^4*(S1 - S3)^3*Sin[i]^4)/(16*a^5*S1^4) - (525*e^2*R^4*(S1 - S3)^3*Sin[i]^4)/(16*a^5*S1^4) - (11025*e^4*R^4*(S1 - S3)^3*Sin[i]^4)/(128*a^5*S1^4)) + C42*((135*e*R^4*Cos[i])/(8*a^5*S1) - (45*e^3*R^4*Cos[i])/(64*a^5*S1) + (315*R^4*(S1 - S3)*Cos[i]^2)/(8*a^5*S1^2) + (135*e*R^4*(S1 - S3)*Cos[i]^2)/(8*a^5*S1^2) + (1575*e^2*R^4*(S1 - S3)*Cos[i]^2)/(8*a^5*S1^2) - (45*e^3*R^4*(S1 - S3)*Cos[i]^2)/(64*a^5*S1^2) + (33075*e^4*R^4*(S1 - S3)*Cos[i]^2)/(64*a^5*S1^2) - (945*e*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]^3)/(16*a^5*S1) + (315*e^3*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]^3)/(128*a^5*S1) + (945*e*R^4*(S1 - S3)^2*Cos[i]^3)/(8*a^5*S1^3) - (315*e^3*R^4*(S1 - S3)^2*Cos[i]^3)/(64*a^5*S1^3) - (315*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^4)/(8*a^5*S1^2) - (945*e*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^4)/(8*a^5*S1^2) - (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^4)/(8*a^5*S1^2) + (315*e^3*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^4)/(64*a^5*S1^2) - (33075*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^4)/(64*a^5*S1^2) - (315*R^4*(S1 - S3)^3*Cos[i]^4)/(8*a^5*S1^4) + (945*e*R^4*(S1 - S3)^3*Cos[i]^4)/(8*a^5*S1^4) - (1575*e^2*R^4*(S1 - S3)^3*Cos[i]^4)/(8*a^5*S1^4) - (315*e^3*R^4*(S1 - S3)^3*Cos[i]^4)/(64*a^5*S1^4) - (33075*e^4*R^4*(S1 - S3)^3*Cos[i]^4)/(64*a^5*S1^4) - (135*e*R^4*(S1 - S3)*Sin[i])/(8*a^5*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) + (45*e^3*R^4*(S1 - S3)*Sin[i])/(64*a^5*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) + (315*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(8*a^5*S1) + (135*e*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(8*a^5*S1) + (1575*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(8*a^5*S1) - (45*e^3*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(64*a^5*S1) + (33075*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(64*a^5*S1) - (315*R^4*(S1 - S3)^2*Cos[i]*Sin[i])/(8*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) - (135*e*R^4*(S1 - S3)^2*Cos[i]*Sin[i])/(8*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) - (1575*e^2*R^4*(S1 - S3)^2*Cos[i]*Sin[i])/(8*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) + (45*e^3*R^4*(S1 - S3)^2*Cos[i]*Sin[i])/(64*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) - (33075*e^4*R^4*(S1 - S3)^2*Cos[i]*Sin[i])/(64*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) + (6615*e*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^2*Sin[i])/(16*a^5*S1^2) - (2205*e^3*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^2*Sin[i])/(128*a^5*S1^2) - (945*e*R^4*(S1 - S3)^3*Cos[i]^2*Sin[i])/(8*a^5*S1^4*Sqrt[1 - (S1 - S3)^2/S1^2]) + (315*e^3*R^4*(S1 - S3)^3*Cos[i]^2*Sin[i])/(64*a^5*S1^4*Sqrt[1 - (S1 - S3)^2/S1^2]) - (315*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(8*a^5*S1) - (945*e*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(8*a^5*S1) - (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(8*a^5*S1) + (315*e^3*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(64*a^5*S1) - (33075*e^4*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(64*a^5*S1) + (2835*e*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^3*Sin[i])/(4*a^5*S1^3) - (945*e^3*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^3*Sin[i])/(32*a^5*S1^3) + (315*R^4*(S1 - S3)^4*Cos[i]^3*Sin[i])/(8*a^5*S1^5*Sqrt[1 - (S1 - S3)^2/S1^2]) - (945*e*R^4*(S1 - S3)^4*Cos[i]^3*Sin[i])/(8*a^5*S1^5*Sqrt[1 - (S1 - S3)^2/S1^2]) + (1575*e^2*R^4*(S1 - S3)^4*Cos[i]^3*Sin[i])/(8*a^5*S1^5*Sqrt[1 - (S1 - S3)^2/S1^2]) + (315*e^3*R^4*(S1 - S3)^4*Cos[i]^3*Sin[i])/(64*a^5*S1^5*Sqrt[1 - (S1 - S3)^2/S1^2]) + (33075*e^4*R^4*(S1 - S3)^4*Cos[i]^3*Sin[i])/(64*a^5*S1^5*Sqrt[1 - (S1 - S3)^2/S1^2]) - (315*R^4*(S1 - S3)*Sin[i]^2)/(8*a^5*S1^2) - (135*e*R^4*(S1 - S3)*Sin[i]^2)/(8*a^5*S1^2) - (1575*e^2*R^4*(S1 - S3)*Sin[i]^2)/(8*a^5*S1^2) + (45*e^3*R^4*(S1 - S3)*Sin[i]^2)/(64*a^5*S1^2) - (33075*e^4*R^4*(S1 - S3)*Sin[i]^2)/(64*a^5*S1^2) + (945*e*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]*Sin[i]^2)/(8*a^5*S1) - (315*e^3*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]*Sin[i]^2)/(64*a^5*S1) - (6615*e*R^4*(S1 - S3)^2*Cos[i]*Sin[i]^2)/(16*a^5*S1^3) + (2205*e^3*R^4*(S1 - S3)^2*Cos[i]*Sin[i]^2)/(128*a^5*S1^3) + (2835*e*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^2*Sin[i]^2)/(4*a^5*S1^2) - (945*e^3*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^2*Sin[i]^2)/(32*a^5*S1^2) - (2835*e*R^4*(S1 - S3)^3*Cos[i]^2*Sin[i]^2)/(4*a^5*S1^4) + (945*e^3*R^4*(S1 - S3)^3*Cos[i]^2*Sin[i]^2)/(32*a^5*S1^4) - (945*e*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Sin[i]^3)/(8*a^5*S1^2) + (315*e^3*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Sin[i]^3)/(64*a^5*S1^2) + (945*e*R^4*(S1 - S3)^3*Sin[i]^3)/(16*a^5*S1^4*Sqrt[1 - (S1 - S3)^2/S1^2]) - (315*e^3*R^4*(S1 - S3)^3*Sin[i]^3)/(128*a^5*S1^4*Sqrt[1 - (S1 - S3)^2/S1^2]) - (315*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]*Sin[i]^3)/(8*a^5*S1) + (945*e*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]*Sin[i]^3)/(8*a^5*S1) - (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]*Sin[i]^3)/(8*a^5*S1) - (315*e^3*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]*Sin[i]^3)/(64*a^5*S1) - (33075*e^4*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]*Sin[i]^3)/(64*a^5*S1) - (2835*e*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]*Sin[i]^3)/(4*a^5*S1^3) + (945*e^3*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]*Sin[i]^3)/(32*a^5*S1^3) + (315*R^4*(S1 - S3)^4*Cos[i]*Sin[i]^3)/(8*a^5*S1^5*Sqrt[1 - (S1 - S3)^2/S1^2]) + (945*e*R^4*(S1 - S3)^4*Cos[i]*Sin[i]^3)/(8*a^5*S1^5*Sqrt[1 - (S1 - S3)^2/S1^2]) + (1575*e^2*R^4*(S1 - S3)^4*Cos[i]*Sin[i]^3)/(8*a^5*S1^5*Sqrt[1 - (S1 - S3)^2/S1^2]) - (315*e^3*R^4*(S1 - S3)^4*Cos[i]*Sin[i]^3)/(64*a^5*S1^5*Sqrt[1 - (S1 - S3)^2/S1^2]) + (33075*e^4*R^4*(S1 - S3)^4*Cos[i]*Sin[i]^3)/(64*a^5*S1^5*Sqrt[1 - (S1 - S3)^2/S1^2]) + (315*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Sin[i]^4)/(8*a^5*S1^2) - (945*e*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Sin[i]^4)/(8*a^5*S1^2) + (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Sin[i]^4)/(8*a^5*S1^2) + (315*e^3*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Sin[i]^4)/(64*a^5*S1^2) + (33075*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Sin[i]^4)/(64*a^5*S1^2) + (315*R^4*(S1 - S3)^3*Sin[i]^4)/(8*a^5*S1^4) + (945*e*R^4*(S1 - S3)^3*Sin[i]^4)/(8*a^5*S1^4) + (1575*e^2*R^4*(S1 - S3)^3*Sin[i]^4)/(8*a^5*S1^4) - (315*e^3*R^4*(S1 - S3)^3*Sin[i]^4)/(64*a^5*S1^4) + (33075*e^4*R^4*(S1 - S3)^3*Sin[i]^4)/(64*a^5*S1^4)) + C44*((-945*e*R^4*Cos[i])/(2*a^5*S1) - (5355*e^2*R^4*Cos[i])/(8*a^5*S1) + (315*e^3*R^4*Cos[i])/(16*a^5*S1) + (33705*e^4*R^4*Cos[i])/(8*a^5*S1) + (315*R^4*(S1 - S3)*Cos[i]^2)/(2*a^5*S1^2) - (945*e*R^4*(S1 - S3)*Cos[i]^2)/(a^5*S1^2) - (9765*e^2*R^4*(S1 - S3)*Cos[i]^2)/(8*a^5*S1^2) + (315*e^3*R^4*(S1 - S3)*Cos[i]^2)/(8*a^5*S1^2) + (235305*e^4*R^4*(S1 - S3)*Cos[i]^2)/(16*a^5*S1^2) + (945*e*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]^3)/(2*a^5*S1) - (315*e^3*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]^3)/(16*a^5*S1) + (945*e*R^4*(S1 - S3)^2*Cos[i]^3)/(2*a^5*S1^3) - (16065*e^2*R^4*(S1 - S3)^2*Cos[i]^3)/(8*a^5*S1^3) - (315*e^3*R^4*(S1 - S3)^2*Cos[i]^3)/(16*a^5*S1^3) + (101115*e^4*R^4*(S1 - S3)^2*Cos[i]^3)/(8*a^5*S1^3) - (315*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^4)/(2*a^5*S1^2) + (945*e*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^4)/(a^5*S1^2) - (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^4)/(2*a^5*S1^2) - (315*e^3*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^4)/(8*a^5*S1^2) - (33075*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^4)/(16*a^5*S1^2) - (315*R^4*(S1 - S3)^3*Cos[i]^4)/(2*a^5*S1^4) + (945*e*R^4*(S1 - S3)^3*Cos[i]^4)/(a^5*S1^4) - (11655*e^2*R^4*(S1 - S3)^3*Cos[i]^4)/(8*a^5*S1^4) - (315*e^3*R^4*(S1 - S3)^3*Cos[i]^4)/(8*a^5*S1^4) + (34335*e^4*R^4*(S1 - S3)^3*Cos[i]^4)/(16*a^5*S1^4) + (945*e*R^4*(S1 - S3)*Sin[i])/(2*a^5*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) + (5355*e^2*R^4*(S1 - S3)*Sin[i])/(8*a^5*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) - (315*e^3*R^4*(S1 - S3)*Sin[i])/(16*a^5*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) - (33705*e^4*R^4*(S1 - S3)*Sin[i])/(8*a^5*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) + (315*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(2*a^5*S1) - (945*e*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(a^5*S1) - (9765*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(8*a^5*S1) + (315*e^3*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(8*a^5*S1) + (235305*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(16*a^5*S1) - (315*R^4*(S1 - S3)^2*Cos[i]*Sin[i])/(2*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) + (945*e*R^4*(S1 - S3)^2*Cos[i]*Sin[i])/(a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) + (9765*e^2*R^4*(S1 - S3)^2*Cos[i]*Sin[i])/(8*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) - (315*e^3*R^4*(S1 - S3)^2*Cos[i]*Sin[i])/(8*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) - (235305*e^4*R^4*(S1 - S3)^2*Cos[i]*Sin[i])/(16*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) - (945*e*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^2*Sin[i])/(2*a^5*S1^2) - (16065*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^2*Sin[i])/(4*a^5*S1^2) + (315*e^3*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^2*Sin[i])/(16*a^5*S1^2) + (101115*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^2*Sin[i])/(4*a^5*S1^2) - (945*e*R^4*(S1 - S3)^3*Cos[i]^2*Sin[i])/(2*a^5*S1^4*Sqrt[1 - (S1 - S3)^2/S1^2]) + (16065*e^2*R^4*(S1 - S3)^3*Cos[i]^2*Sin[i])/(8*a^5*S1^4*Sqrt[1 - (S1 - S3)^2/S1^2]) + (315*e^3*R^4*(S1 - S3)^3*Cos[i]^2*Sin[i])/(16*a^5*S1^4*Sqrt[1 - (S1 - S3)^2/S1^2]) - (101115*e^4*R^4*(S1 - S3)^3*Cos[i]^2*Sin[i])/(8*a^5*S1^4*Sqrt[1 - (S1 - S3)^2/S1^2]) - (315*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(2*a^5*S1) + (945*e*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(a^5*S1) - (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(2*a^5*S1) - (315*e^3*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(8*a^5*S1) - (33075*e^4*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(16*a^5*S1) - (16065*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^3*Sin[i])/(8*a^5*S1^3) + (101115*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^3*Sin[i])/(8*a^5*S1^3) + (315*R^4*(S1 - S3)^4*Cos[i]^3*Sin[i])/(2*a^5*S1^5*Sqrt[1 - (S1 - S3)^2/S1^2]) - (945*e*R^4*(S1 - S3)^4*Cos[i]^3*Sin[i])/(a^5*S1^5*Sqrt[1 - (S1 - S3)^2/S1^2]) + (11655*e^2*R^4*(S1 - S3)^4*Cos[i]^3*Sin[i])/(8*a^5*S1^5*Sqrt[1 - (S1 - S3)^2/S1^2]) + (315*e^3*R^4*(S1 - S3)^4*Cos[i]^3*Sin[i])/(8*a^5*S1^5*Sqrt[1 - (S1 - S3)^2/S1^2]) - (34335*e^4*R^4*(S1 - S3)^4*Cos[i]^3*Sin[i])/(16*a^5*S1^5*Sqrt[1 - (S1 - S3)^2/S1^2]) - (315*R^4*(S1 - S3)*Sin[i]^2)/(2*a^5*S1^2) + (945*e*R^4*(S1 - S3)*Sin[i]^2)/(a^5*S1^2) + (9765*e^2*R^4*(S1 - S3)*Sin[i]^2)/(8*a^5*S1^2) - (315*e^3*R^4*(S1 - S3)*Sin[i]^2)/(8*a^5*S1^2) - (235305*e^4*R^4*(S1 - S3)*Sin[i]^2)/(16*a^5*S1^2) + (945*e*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]*Sin[i]^2)/(2*a^5*S1) - (16065*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]*Sin[i]^2)/(8*a^5*S1) - (315*e^3*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]*Sin[i]^2)/(16*a^5*S1) + (101115*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]*Sin[i]^2)/(8*a^5*S1) + (945*e*R^4*(S1 - S3)^2*Cos[i]*Sin[i]^2)/(2*a^5*S1^3) + (16065*e^2*R^4*(S1 - S3)^2*Cos[i]*Sin[i]^2)/(4*a^5*S1^3) - (315*e^3*R^4*(S1 - S3)^2*Cos[i]*Sin[i]^2)/(16*a^5*S1^3) - (101115*e^4*R^4*(S1 - S3)^2*Cos[i]*Sin[i]^2)/(4*a^5*S1^3) - (16065*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^2*Sin[i]^2)/(8*a^5*S1^2) + (101115*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^2*Sin[i]^2)/(8*a^5*S1^2) + (16065*e^2*R^4*(S1 - S3)^3*Cos[i]^2*Sin[i]^2)/(8*a^5*S1^4) - (101115*e^4*R^4*(S1 - S3)^3*Cos[i]^2*Sin[i]^2)/(8*a^5*S1^4) - (945*e*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Sin[i]^3)/(2*a^5*S1^2) + (16065*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Sin[i]^3)/(8*a^5*S1^2) + (315*e^3*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Sin[i]^3)/(16*a^5*S1^2) - (101115*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Sin[i]^3)/(8*a^5*S1^2) - (945*e*R^4*(S1 - S3)^3*Sin[i]^3)/(2*a^5*S1^4*Sqrt[1 - (S1 - S3)^2/S1^2]) + (315*e^3*R^4*(S1 - S3)^3*Sin[i]^3)/(16*a^5*S1^4*Sqrt[1 - (S1 - S3)^2/S1^2]) - (315*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]*Sin[i]^3)/(2*a^5*S1) + (945*e*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]*Sin[i]^3)/(a^5*S1) - (11655*e^2*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]*Sin[i]^3)/(8*a^5*S1) - (315*e^3*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]*Sin[i]^3)/(8*a^5*S1) + (34335*e^4*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]*Sin[i]^3)/(16*a^5*S1) + (16065*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]*Sin[i]^3)/(8*a^5*S1^3) - (101115*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]*Sin[i]^3)/(8*a^5*S1^3) + (315*R^4*(S1 - S3)^4*Cos[i]*Sin[i]^3)/(2*a^5*S1^5*Sqrt[1 - (S1 - S3)^2/S1^2]) - (945*e*R^4*(S1 - S3)^4*Cos[i]*Sin[i]^3)/(a^5*S1^5*Sqrt[1 - (S1 - S3)^2/S1^2]) + (1575*e^2*R^4*(S1 - S3)^4*Cos[i]*Sin[i]^3)/(2*a^5*S1^5*Sqrt[1 - (S1 - S3)^2/S1^2]) + (315*e^3*R^4*(S1 - S3)^4*Cos[i]*Sin[i]^3)/(8*a^5*S1^5*Sqrt[1 - (S1 - S3)^2/S1^2]) + (33075*e^4*R^4*(S1 - S3)^4*Cos[i]*Sin[i]^3)/(16*a^5*S1^5*Sqrt[1 - (S1 - S3)^2/S1^2]) + (315*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Sin[i]^4)/(2*a^5*S1^2) - (945*e*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Sin[i]^4)/(a^5*S1^2) + (11655*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Sin[i]^4)/(8*a^5*S1^2) + (315*e^3*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Sin[i]^4)/(8*a^5*S1^2) - (34335*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Sin[i]^4)/(16*a^5*S1^2) + (315*R^4*(S1 - S3)^3*Sin[i]^4)/(2*a^5*S1^4) - (945*e*R^4*(S1 - S3)^3*Sin[i]^4)/(a^5*S1^4) + (1575*e^2*R^4*(S1 - S3)^3*Sin[i]^4)/(2*a^5*S1^4) + (315*e^3*R^4*(S1 - S3)^3*Sin[i]^4)/(8*a^5*S1^4) + (33075*e^4*R^4*(S1 - S3)^3*Sin[i]^4)/(16*a^5*S1^4))) == 0}