sys21 = {-2*n - omd - Omd + S1/C - G*m*M*(C11*((e*R*Cos[i])/(a^2*S1) - (3*e^3*R*Cos[i])/(4*a^2*S1) - (e*R*(S1 - S3)*Cos[i])/(a^2*S1^2) + (3*e^3*R*(S1 - S3)*Cos[i])/(4*a^2*S1^2) + (e*R*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i])/(2*a^2*Sqrt[1 - (S1 - S3)^2/S1^2]) - (3*e^3*R*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i])/(8*a^2*Sqrt[1 - (S1 - S3)^2/S1^2])) + C22*((51*e^2*R^2*Cos[i])/(4*a^3*S1) - (115*e^4*R^2*Cos[i])/(4*a^3*S1) - (51*e^2*R^2*(S1 - S3)*Cos[i])/(4*a^3*S1^2) + (115*e^4*R^2*(S1 - S3)*Cos[i])/(4*a^3*S1^2) + (51*e^2*R^2*(S1 - S3)*Cos[i]^2)/(4*a^3*S1^2) - (115*e^4*R^2*(S1 - S3)*Cos[i]^2)/(4*a^3*S1^2) - (51*e^2*R^2*(S1 - S3)^2*Cos[i]^2)/(4*a^3*S1^3) + (115*e^4*R^2*(S1 - S3)^2*Cos[i]^2)/(4*a^3*S1^3) + (51*e^2*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]) - (115*e^4*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]) + (51*e^2*R^2*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(4*a^3*S1) - (115*e^4*R^2*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(4*a^3*S1) + (51*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]) - (115*e^4*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]) - (51*e^2*R^2*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i])/(4*a^3*S1^2) + (115*e^4*R^2*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i])/(4*a^3*S1^2) + (51*e^2*R^2*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i]^2)/(8*a^3) - (115*e^4*R^2*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i]^2)/(8*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)) + C33*((7335*e^3*R^3*Cos[i])/(32*a^4*S1) - (7335*e^3*R^3*(S1 - S3)*Cos[i])/(32*a^4*S1^2) + (7335*e^3*R^3*(S1 - S3)*Cos[i]^2)/(16*a^4*S1^2) - (7335*e^3*R^3*(S1 - S3)^2*Cos[i]^2)/(16*a^4*S1^3) + (7335*e^3*R^3*(S1 - S3)^2*Cos[i]^3)/(32*a^4*S1^3) - (7335*e^3*R^3*(S1 - S3)^3*Cos[i]^3)/(32*a^4*S1^4) + (7335*e^3*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i])/(64*a^4*Sqrt[1 - (S1 - S3)^2/S1^2]) + (7335*e^3*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(16*a^4*S1) + (7335*e^3*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i])/(32*a^4*S1*Sqrt[1 - (S1 - S3)^2/S1^2]) - (7335*e^3*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i])/(16*a^4*S1^2) + (7335*e^3*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^2*Sin[i])/(16*a^4*S1^2) + (7335*e^3*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^2*Sin[i])/(64*a^4*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) - (7335*e^3*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^2*Sin[i])/(16*a^4*S1^3) + (7335*e^3*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i]^2)/(32*a^4) + (7335*e^3*R^3*(1 - (S1 - S3)^2/S1^2)*Cos[i]*Sin[i]^2)/(32*a^4*S1) + (7335*e^3*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i]^2)/(32*a^4*S1) - (7335*e^3*R^3*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]*Sin[i]^2)/(32*a^4*S1^2) + (7335*e^3*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sqrt[1 - (S1 - S3)^2/S1^2]*Sin[i]^3)/(64*a^4)) + C31*((-9*e*R^3*Cos[i])/(4*a^4*S1) - (33*e^3*R^3*Cos[i])/(16*a^4*S1) + (9*e*R^3*(S1 - S3)*Cos[i])/(4*a^4*S1^2) + (33*e^3*R^3*(S1 - S3)*Cos[i])/(16*a^4*S1^2) + (45*e*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Cos[i]^2)/(16*a^4) + (165*e^3*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Cos[i]^2)/(64*a^4) + (135*e*R^3*(1 - (S1 - S3)^2/S1^2)*Cos[i]^3)/(16*a^4*S1) + (495*e^3*R^3*(1 - (S1 - S3)^2/S1^2)*Cos[i]^3)/(64*a^4*S1) + (135*e*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]^3)/(16*a^4*S1) + (495*e^3*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]^3)/(64*a^4*S1) - (135*e*R^3*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^3)/(16*a^4*S1^2) - (495*e^3*R^3*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^3)/(64*a^4*S1^2) - (9*e*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i])/(8*a^4*Sqrt[1 - (S1 - S3)^2/S1^2]) - (33*e^3*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i])/(32*a^4*Sqrt[1 - (S1 - S3)^2/S1^2]) - (45*e*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(8*a^4*S1) - (165*e^3*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(32*a^4*S1) - (45*e*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i])/(16*a^4*S1*Sqrt[1 - (S1 - S3)^2/S1^2]) - (165*e^3*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i])/(64*a^4*S1*Sqrt[1 - (S1 - S3)^2/S1^2]) + (45*e*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i])/(8*a^4*S1^2) + (165*e^3*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i])/(32*a^4*S1^2) + (405*e*R^3*((-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^4) + (1485*e^3*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]^2*Sin[i])/(128*a^4) - (135*e*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^2*Sin[i])/(4*a^4*S1^2) - (495*e^3*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^2*Sin[i])/(16*a^4*S1^2) - (135*e*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^2*Sin[i])/(16*a^4*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) - (495*e^3*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^2*Sin[i])/(64*a^4*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) + (135*e*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^2*Sin[i])/(4*a^4*S1^3) + (495*e^3*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^2*Sin[i])/(16*a^4*S1^3) + (45*e*R^3*(S1 - S3)*Sin[i]^2)/(8*a^4*S1^2) + (165*e^3*R^3*(S1 - S3)*Sin[i]^2)/(32*a^4*S1^2) - (45*e*R^3*(S1 - S3)^2*Sin[i]^2)/(8*a^4*S1^3) - (165*e^3*R^3*(S1 - S3)^2*Sin[i]^2)/(32*a^4*S1^3) - (135*e*R^3*(1 - (S1 - S3)^2/S1^2)*Cos[i]*Sin[i]^2)/(8*a^4*S1) - (495*e^3*R^3*(1 - (S1 - S3)^2/S1^2)*Cos[i]*Sin[i]^2)/(32*a^4*S1) - (135*e*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i]^2)/(8*a^4*S1) - (495*e^3*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i]^2)/(32*a^4*S1) + (135*e*R^3*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]*Sin[i]^2)/(8*a^4*S1^2) + (495*e^3*R^3*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]*Sin[i]^2)/(32*a^4*S1^2) + (405*e*R^3*(S1 - S3)^2*Cos[i]*Sin[i]^2)/(16*a^4*S1^3) + (1485*e^3*R^3*(S1 - S3)^2*Cos[i]*Sin[i]^2)/(64*a^4*S1^3) - (405*e*R^3*(S1 - S3)^3*Cos[i]*Sin[i]^2)/(16*a^4*S1^4) - (1485*e^3*R^3*(S1 - S3)^3*Cos[i]*Sin[i]^2)/(64*a^4*S1^4) + (135*e*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Sin[i]^3)/(8*a^4*S1^2) + (495*e^3*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Sin[i]^3)/(32*a^4*S1^2) + (135*e*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Sin[i]^3)/(32*a^4*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) + (495*e^3*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Sin[i]^3)/(128*a^4*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) - (135*e*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Sin[i]^3)/(8*a^4*S1^3) - (495*e^3*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Sin[i]^3)/(32*a^4*S1^3)) + C44*((5565*e^2*R^4*Cos[i])/(4*a^5*S1) + (147945*e^4*R^4*Cos[i])/(32*a^5*S1) - (5565*e^2*R^4*(S1 - S3)*Cos[i])/(4*a^5*S1^2) - (147945*e^4*R^4*(S1 - S3)*Cos[i])/(32*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) + (1470*e^2*R^4*(S1 - S3)*Cos[i]^2)/(a^5*S1^2) + (358645*e^4*R^4*(S1 - S3)*Cos[i]^2)/(32*a^5*S1^2) + (525*R^4*(S1 - S3)^2*Cos[i]^2)/(2*a^5*S1^3) - (1470*e^2*R^4*(S1 - S3)^2*Cos[i]^2)/(a^5*S1^3) - (358645*e^4*R^4*(S1 - S3)^2*Cos[i]^2)/(32*a^5*S1^3) - (5565*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]^3)/(4*a^5*S1) + (6265*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]^3)/(8*a^5*S1) - (5565*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]^3)/(4*a^5*S1) + (6265*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]^3)/(8*a^5*S1) + (5565*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^3)/(4*a^5*S1^2) - (6265*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^3)/(8*a^5*S1^2) - (16695*e^2*R^4*(S1 - S3)^2*Cos[i]^3)/(4*a^5*S1^3) + (594195*e^4*R^4*(S1 - S3)^2*Cos[i]^3)/(32*a^5*S1^3) + (16695*e^2*R^4*(S1 - S3)^3*Cos[i]^3)/(4*a^5*S1^4) - (594195*e^4*R^4*(S1 - S3)^3*Cos[i]^3)/(32*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) - (2415*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^4)/(2*a^5*S1^2) + (45605*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) - (2415*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^4)/(4*a^5*S1^2) + (45605*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) + (2415*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Cos[i]^4)/(2*a^5*S1^3) - (45605*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) - (6405*e^2*R^4*(S1 - S3)^3*Cos[i]^4)/(2*a^5*S1^4) + (471695*e^4*R^4*(S1 - S3)^3*Cos[i]^4)/(32*a^5*S1^4) - (945*R^4*(S1 - S3)^4*Cos[i]^4)/(2*a^5*S1^5) + (6405*e^2*R^4*(S1 - S3)^4*Cos[i]^4)/(2*a^5*S1^5) - (471695*e^4*R^4*(S1 - S3)^4*Cos[i]^4)/(32*a^5*S1^5) + (5565*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i])/(8*a^5*Sqrt[1 - (S1 - S3)^2/S1^2]) + (147945*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i])/(64*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) + (1995*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(a^5*S1) + (402745*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(32*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]) + (1995*e^2*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]) + (402745*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]) + (315*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i])/(2*a^5*S1^2) - (1995*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i])/(a^5*S1^2) - (402745*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i])/(32*a^5*S1^2) - (16695*e^2*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])/(8*a^5) + (18795*e^4*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])/(16*a^5) - (5565*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^2*Sin[i])/(2*a^5*S1^2) + (544075*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^2*Sin[i])/(16*a^5*S1^2) - (5565*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^2*Sin[i])/(8*a^5*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) + (544075*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^2*Sin[i])/(64*a^5*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) + (5565*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^2*Sin[i])/(2*a^5*S1^3) - (544075*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^2*Sin[i])/(16*a^5*S1^3) + (315*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(2*a^5*S1) - (1995*e^2*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(a^5*S1) + (58135*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) - (5985*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])/(2*a^5*S1) + (174405*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) + (1995*e^2*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*(S1 - S3)*Cos[i]^3*Sin[i])/(a^5*S1^2) - (58135*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) - (5985*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^3*Sin[i])/(a^5*S1^3) + (867825*e^4*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*((-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]) - (1995*e^2*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]) + (289275*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^3*Cos[i]^3*Sin[i])/(64*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) + (5985*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^3*Cos[i]^3*Sin[i])/(a^5*S1^4) - (867825*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^3*Cos[i]^3*Sin[i])/(32*a^5*S1^4) - (525*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i]^2)/(4*a^5) + (735*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i]^2)/a^5 + (358645*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i]^2)/(64*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) - (5565*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]*Sin[i]^2)/(4*a^5*S1) + (544075*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]*Sin[i]^2)/(32*a^5*S1) - (5565*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i]^2)/(4*a^5*S1) + (544075*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i]^2)/(32*a^5*S1) + (5565*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]*Sin[i]^2)/(4*a^5*S1^2) - (544075*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]*Sin[i]^2)/(32*a^5*S1^2) - (16695*e^2*R^4*(S1 - S3)^2*Cos[i]*Sin[i]^2)/(4*a^5*S1^3) + (18795*e^4*R^4*(S1 - S3)^2*Cos[i]*Sin[i]^2)/(8*a^5*S1^3) + (16695*e^2*R^4*(S1 - S3)^3*Cos[i]*Sin[i]^2)/(4*a^5*S1^4) - (18795*e^4*R^4*(S1 - S3)^3*Cos[i]*Sin[i]^2)/(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]^2*Sin[i]^2)/a^5 - (2415*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)/(2*a^5) + (45605*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) - (2415*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^2*Sin[i]^2)/(a^5*S1^2) + (883855*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)/(a^5*S1^2) - (2415*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^2*Sin[i]^2)/(2*a^5*S1^2) + (883855*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) - (630*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Cos[i]^2*Sin[i]^2)/(a^5*S1^3) + (2415*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Cos[i]^2*Sin[i]^2)/(a^5*S1^3) - (883855*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Cos[i]^2*Sin[i]^2)/(32*a^5*S1^3) + (630*R^4*(S1 - S3)^3*Cos[i]^2*Sin[i]^2)/(a^5*S1^4) - (2415*e^2*R^4*(S1 - S3)^3*Cos[i]^2*Sin[i]^2)/(a^5*S1^4) + (45605*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) + (2415*e^2*R^4*(S1 - S3)^4*Cos[i]^2*Sin[i]^2)/(a^5*S1^5) - (45605*e^4*R^4*(S1 - S3)^4*Cos[i]^2*Sin[i]^2)/(4*a^5*S1^5) - (16695*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)/(8*a^5) + (594195*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)/(64*a^5) - (5565*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Sin[i]^3)/(2*a^5*S1^2) + (6265*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Sin[i]^3)/(4*a^5*S1^2) - (5565*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Sin[i]^3)/(8*a^5*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) + (6265*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Sin[i]^3)/(16*a^5*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) + (5565*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Sin[i]^3)/(2*a^5*S1^3) - (6265*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Sin[i]^3)/(4*a^5*S1^3) + (315*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]*Sin[i]^3)/(2*a^5*S1) - (1995*e^2*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]*Sin[i]^3)/(a^5*S1) + (289275*e^4*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]*Sin[i]^3)/(32*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) - (5985*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)/(2*a^5*S1) + (867825*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)/(64*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) + (1995*e^2*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*(S1 - S3)*Cos[i]*Sin[i]^3)/(a^5*S1^2) - (289275*e^4*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*(S1 - S3)*Cos[i]*Sin[i]^3)/(32*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) - (5985*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]*Sin[i]^3)/(a^5*S1^3) + (174405*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]) - (1995*e^2*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]) + (58135*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) + (5985*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^3*Cos[i]*Sin[i]^3)/(a^5*S1^4) - (174405*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) - (6405*e^2*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) + (471695*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(1 - (S1 - S3)^2/S1^2)*Sin[i]^4)/(64*a^5) + (315*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Sin[i]^4)/(a^5*S1^2) - (2415*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Sin[i]^4)/(2*a^5*S1^2) + (45605*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) - (2415*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Sin[i]^4)/(4*a^5*S1^2) + (45605*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) + (2415*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Sin[i]^4)/(2*a^5*S1^3) - (45605*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*((-795*e^2*R^4*Cos[i])/(16*a^5*S1) + (895*e^4*R^4*Cos[i])/(32*a^5*S1) + (795*e^2*R^4*(S1 - S3)*Cos[i])/(16*a^5*S1^2) - (895*e^4*R^4*(S1 - S3)*Cos[i])/(32*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) - (195*e^2*R^4*(S1 - S3)*Cos[i]^2)/(16*a^5*S1^2) + (4045*e^4*R^4*(S1 - S3)*Cos[i]^2)/(32*a^5*S1^2) - (15*R^4*(S1 - S3)^2*Cos[i]^2)/(2*a^5*S1^3) + (195*e^2*R^4*(S1 - S3)^2*Cos[i]^2)/(16*a^5*S1^3) - (4045*e^4*R^4*(S1 - S3)^2*Cos[i]^2)/(32*a^5*S1^3) + (5565*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]^3)/(32*a^5*S1) - (6265*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]^3)/(64*a^5*S1) + (5565*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]^3)/(32*a^5*S1) - (6265*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]^3)/(64*a^5*S1) - (5565*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^3)/(32*a^5*S1^2) + (6265*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^3)/(64*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) + (2415*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^4)/(16*a^5*S1^2) - (45605*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) + (2415*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^4)/(32*a^5*S1^2) - (45605*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) - (2415*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Cos[i]^4)/(16*a^5*S1^3) + (45605*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Cos[i]^4)/(64*a^5*S1^3) - (795*e^2*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]) + (895*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i])/(64*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) - (3945*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(16*a^5*S1) - (31285*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]) - (3945*e^2*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]) - (31285*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) + (3945*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i])/(16*a^5*S1^2) + (31285*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i])/(64*a^5*S1^2) + (16695*e^2*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])/(64*a^5) - (18795*e^4*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])/(128*a^5) - (5565*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^2*Sin[i])/(8*a^5*S1^2) + (6265*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^2*Sin[i])/(16*a^5*S1^2) - (5565*e^2*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]) + (6265*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^2*Sin[i])/(64*a^5*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) + (5565*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^2*Sin[i])/(8*a^5*S1^3) - (6265*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^2*Sin[i])/(16*a^5*S1^3) + (315*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(8*a^5*S1) + (8715*e^2*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(16*a^5*S1) + (20545*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) + (26145*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) + (61635*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) - (8715*e^2*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*(S1 - S3)*Cos[i]^3*Sin[i])/(16*a^5*S1^2) - (20545*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) - (7245*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^3*Sin[i])/(16*a^5*S1^3) + (136815*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]) - (2415*e^2*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]) + (45605*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) + (7245*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^3*Cos[i]^3*Sin[i])/(16*a^5*S1^4) - (136815*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) - (195*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i]^2)/(32*a^5) + (4045*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i]^2)/(64*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) - (5565*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]*Sin[i]^2)/(16*a^5*S1) + (6265*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]*Sin[i]^2)/(32*a^5*S1) - (5565*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i]^2)/(16*a^5*S1) + (6265*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i]^2)/(32*a^5*S1) + (5565*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]*Sin[i]^2)/(16*a^5*S1^2) - (6265*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]*Sin[i]^2)/(32*a^5*S1^2) + (16695*e^2*R^4*(S1 - S3)^2*Cos[i]*Sin[i]^2)/(32*a^5*S1^3) - (18795*e^4*R^4*(S1 - S3)^2*Cos[i]*Sin[i]^2)/(64*a^5*S1^3) - (16695*e^2*R^4*(S1 - S3)^3*Cos[i]*Sin[i]^2)/(32*a^5*S1^4) + (18795*e^4*R^4*(S1 - S3)^3*Cos[i]*Sin[i]^2)/(64*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) + (2415*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)/(16*a^5) - (45605*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) - (1785*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^2*Sin[i]^2)/(a^5*S1^2) - (8015*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) - (1785*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^2*Sin[i]^2)/(2*a^5*S1^2) - (8015*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) + (1785*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Cos[i]^2*Sin[i]^2)/(a^5*S1^3) + (8015*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) + (2415*e^2*R^4*(S1 - S3)^3*Cos[i]^2*Sin[i]^2)/(8*a^5*S1^4) - (45605*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) - (2415*e^2*R^4*(S1 - S3)^4*Cos[i]^2*Sin[i]^2)/(8*a^5*S1^5) + (45605*e^4*R^4*(S1 - S3)^4*Cos[i]^2*Sin[i]^2)/(32*a^5*S1^5) + (5565*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Sin[i]^3)/(16*a^5*S1^2) - (6265*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Sin[i]^3)/(32*a^5*S1^2) + (5565*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Sin[i]^3)/(64*a^5*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) - (6265*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Sin[i]^3)/(128*a^5*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) - (5565*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Sin[i]^3)/(16*a^5*S1^3) + (6265*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Sin[i]^3)/(32*a^5*S1^3) + (315*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]*Sin[i]^3)/(8*a^5*S1) - (2415*e^2*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]*Sin[i]^3)/(16*a^5*S1) + (45605*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) - (7245*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)/(32*a^5*S1) + (136815*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) + (2415*e^2*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*(S1 - S3)*Cos[i]*Sin[i]^3)/(16*a^5*S1^2) - (45605*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) + (26145*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]*Sin[i]^3)/(16*a^5*S1^3) + (61635*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]) + (8715*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]) + (20545*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) - (26145*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^3*Cos[i]*Sin[i]^3)/(16*a^5*S1^4) - (61635*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) + (2415*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Sin[i]^4)/(16*a^5*S1^2) - (45605*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) + (2415*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Sin[i]^4)/(32*a^5*S1^2) - (45605*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) - (2415*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Sin[i]^4)/(16*a^5*S1^3) + (45605*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*(C11*(-((e*R*Cos[i])/(a^2*S1)) + (3*e^3*R*Cos[i])/(4*a^2*S1) + (e*R*(S1 - S3)*Sin[i])/(a^2*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) - (3*e^3*R*(S1 - S3)*Sin[i])/(4*a^2*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2])) + C22*((-51*e^2*R^2*Cos[i])/(4*a^3*S1) + (115*e^4*R^2*Cos[i])/(4*a^3*S1) - (51*e^2*R^2*(S1 - S3)*Cos[i]^2)/(4*a^3*S1^2) + (115*e^4*R^2*(S1 - S3)*Cos[i]^2)/(4*a^3*S1^2) + (51*e^2*R^2*(S1 - S3)*Sin[i])/(4*a^3*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) - (115*e^4*R^2*(S1 - S3)*Sin[i])/(4*a^3*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) - (51*e^2*R^2*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(4*a^3*S1) + (115*e^4*R^2*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(4*a^3*S1) + (51*e^2*R^2*(S1 - S3)^2*Cos[i]*Sin[i])/(4*a^3*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) - (115*e^4*R^2*(S1 - S3)^2*Cos[i]*Sin[i])/(4*a^3*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) + (51*e^2*R^2*(S1 - S3)*Sin[i]^2)/(4*a^3*S1^2) - (115*e^4*R^2*(S1 - S3)*Sin[i]^2)/(4*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)) + C33*((-7335*e^3*R^3*Cos[i])/(32*a^4*S1) - (7335*e^3*R^3*(S1 - S3)*Cos[i]^2)/(16*a^4*S1^2) - (7335*e^3*R^3*(S1 - S3)^2*Cos[i]^3)/(32*a^4*S1^3) + (7335*e^3*R^3*(S1 - S3)*Sin[i])/(32*a^4*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) - (7335*e^3*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(16*a^4*S1) + (7335*e^3*R^3*(S1 - S3)^2*Cos[i]*Sin[i])/(16*a^4*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) - (7335*e^3*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^2*Sin[i])/(16*a^4*S1^2) + (7335*e^3*R^3*(S1 - S3)^3*Cos[i]^2*Sin[i])/(32*a^4*S1^4*Sqrt[1 - (S1 - S3)^2/S1^2]) + (7335*e^3*R^3*(S1 - S3)*Sin[i]^2)/(16*a^4*S1^2) - (7335*e^3*R^3*(1 - (S1 - S3)^2/S1^2)*Cos[i]*Sin[i]^2)/(32*a^4*S1) + (7335*e^3*R^3*(S1 - S3)^2*Cos[i]*Sin[i]^2)/(16*a^4*S1^3) + (7335*e^3*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Sin[i]^3)/(32*a^4*S1^2)) + C31*((9*e*R^3*Cos[i])/(4*a^4*S1) + (33*e^3*R^3*Cos[i])/(16*a^4*S1) + (45*e*R^3*(S1 - S3)*Cos[i]^2)/(8*a^4*S1^2) + (165*e^3*R^3*(S1 - S3)*Cos[i]^2)/(32*a^4*S1^2) - (135*e*R^3*(1 - (S1 - S3)^2/S1^2)*Cos[i]^3)/(16*a^4*S1) - (495*e^3*R^3*(1 - (S1 - S3)^2/S1^2)*Cos[i]^3)/(64*a^4*S1) + (135*e*R^3*(S1 - S3)^2*Cos[i]^3)/(8*a^4*S1^3) + (495*e^3*R^3*(S1 - S3)^2*Cos[i]^3)/(32*a^4*S1^3) - (9*e*R^3*(S1 - S3)*Sin[i])/(4*a^4*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) - (33*e^3*R^3*(S1 - S3)*Sin[i])/(16*a^4*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) + (45*e*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(8*a^4*S1) + (165*e^3*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(32*a^4*S1) - (45*e*R^3*(S1 - S3)^2*Cos[i]*Sin[i])/(8*a^4*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) - (165*e^3*R^3*(S1 - S3)^2*Cos[i]*Sin[i])/(32*a^4*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) + (945*e*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^2*Sin[i])/(16*a^4*S1^2) + (3465*e^3*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^2*Sin[i])/(64*a^4*S1^2) - (135*e*R^3*(S1 - S3)^3*Cos[i]^2*Sin[i])/(8*a^4*S1^4*Sqrt[1 - (S1 - S3)^2/S1^2]) - (495*e^3*R^3*(S1 - S3)^3*Cos[i]^2*Sin[i])/(32*a^4*S1^4*Sqrt[1 - (S1 - S3)^2/S1^2]) - (45*e*R^3*(S1 - S3)*Sin[i]^2)/(8*a^4*S1^2) - (165*e^3*R^3*(S1 - S3)*Sin[i]^2)/(32*a^4*S1^2) + (135*e*R^3*(1 - (S1 - S3)^2/S1^2)*Cos[i]*Sin[i]^2)/(8*a^4*S1) + (495*e^3*R^3*(1 - (S1 - S3)^2/S1^2)*Cos[i]*Sin[i]^2)/(32*a^4*S1) - (945*e*R^3*(S1 - S3)^2*Cos[i]*Sin[i]^2)/(16*a^4*S1^3) - (3465*e^3*R^3*(S1 - S3)^2*Cos[i]*Sin[i]^2)/(64*a^4*S1^3) - (135*e*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Sin[i]^3)/(8*a^4*S1^2) - (495*e^3*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Sin[i]^3)/(32*a^4*S1^2) + (135*e*R^3*(S1 - S3)^3*Sin[i]^3)/(16*a^4*S1^4*Sqrt[1 - (S1 - S3)^2/S1^2]) + (495*e^3*R^3*(S1 - S3)^3*Sin[i]^3)/(64*a^4*S1^4*Sqrt[1 - (S1 - S3)^2/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*((795*e^2*R^4*Cos[i])/(16*a^5*S1) - (895*e^4*R^4*Cos[i])/(32*a^5*S1) + (315*R^4*(S1 - S3)*Cos[i]^2)/(8*a^5*S1^2) + (3945*e^2*R^4*(S1 - S3)*Cos[i]^2)/(16*a^5*S1^2) + (31285*e^4*R^4*(S1 - S3)*Cos[i]^2)/(64*a^5*S1^2) - (5565*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]^3)/(32*a^5*S1) + (6265*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]^3)/(64*a^5*S1) + (5565*e^2*R^4*(S1 - S3)^2*Cos[i]^3)/(16*a^5*S1^3) - (6265*e^4*R^4*(S1 - S3)^2*Cos[i]^3)/(32*a^5*S1^3) - (315*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^4)/(8*a^5*S1^2) - (8715*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^4)/(16*a^5*S1^2) - (20545*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) + (2415*e^2*R^4*(S1 - S3)^3*Cos[i]^4)/(16*a^5*S1^4) - (45605*e^4*R^4*(S1 - S3)^3*Cos[i]^4)/(64*a^5*S1^4) - (795*e^2*R^4*(S1 - S3)*Sin[i])/(16*a^5*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) + (895*e^4*R^4*(S1 - S3)*Sin[i])/(32*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) + (3945*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(16*a^5*S1) + (31285*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]) - (3945*e^2*R^4*(S1 - S3)^2*Cos[i]*Sin[i])/(16*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) - (31285*e^4*R^4*(S1 - S3)^2*Cos[i]*Sin[i])/(64*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) + (38955*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^2*Sin[i])/(32*a^5*S1^2) - (43855*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^2*Sin[i])/(64*a^5*S1^2) - (5565*e^2*R^4*(S1 - S3)^3*Cos[i]^2*Sin[i])/(16*a^5*S1^4*Sqrt[1 - (S1 - S3)^2/S1^2]) + (6265*e^4*R^4*(S1 - S3)^3*Cos[i]^2*Sin[i])/(32*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) - (8715*e^2*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(16*a^5*S1) - (20545*e^4*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(64*a^5*S1) + (16695*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^3*Sin[i])/(8*a^5*S1^3) - (18795*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*(S1 - S3)^4*Cos[i]^3*Sin[i])/(8*a^5*S1^5*Sqrt[1 - (S1 - S3)^2/S1^2]) - (2415*e^2*R^4*(S1 - S3)^4*Cos[i]^3*Sin[i])/(16*a^5*S1^5*Sqrt[1 - (S1 - S3)^2/S1^2]) + (45605*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) - (3945*e^2*R^4*(S1 - S3)*Sin[i]^2)/(16*a^5*S1^2) - (31285*e^4*R^4*(S1 - S3)*Sin[i]^2)/(64*a^5*S1^2) + (5565*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]*Sin[i]^2)/(16*a^5*S1) - (6265*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]*Sin[i]^2)/(32*a^5*S1) - (38955*e^2*R^4*(S1 - S3)^2*Cos[i]*Sin[i]^2)/(32*a^5*S1^3) + (43855*e^4*R^4*(S1 - S3)^2*Cos[i]*Sin[i]^2)/(64*a^5*S1^3) + (16695*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^2*Sin[i]^2)/(8*a^5*S1^2) - (18795*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^2*Sin[i]^2)/(16*a^5*S1^2) - (16695*e^2*R^4*(S1 - S3)^3*Cos[i]^2*Sin[i]^2)/(8*a^5*S1^4) + (18795*e^4*R^4*(S1 - S3)^3*Cos[i]^2*Sin[i]^2)/(16*a^5*S1^4) - (5565*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Sin[i]^3)/(16*a^5*S1^2) + (6265*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Sin[i]^3)/(32*a^5*S1^2) + (5565*e^2*R^4*(S1 - S3)^3*Sin[i]^3)/(32*a^5*S1^4*Sqrt[1 - (S1 - S3)^2/S1^2]) - (6265*e^4*R^4*(S1 - S3)^3*Sin[i]^3)/(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]*Sin[i]^3)/(8*a^5*S1) + (2415*e^2*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]*Sin[i]^3)/(16*a^5*S1) - (45605*e^4*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]*Sin[i]^3)/(64*a^5*S1) - (16695*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]*Sin[i]^3)/(8*a^5*S1^3) + (18795*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*(S1 - S3)^4*Cos[i]*Sin[i]^3)/(8*a^5*S1^5*Sqrt[1 - (S1 - S3)^2/S1^2]) + (8715*e^2*R^4*(S1 - S3)^4*Cos[i]*Sin[i]^3)/(16*a^5*S1^5*Sqrt[1 - (S1 - S3)^2/S1^2]) + (20545*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) - (2415*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Sin[i]^4)/(16*a^5*S1^2) + (45605*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) + (8715*e^2*R^4*(S1 - S3)^3*Sin[i]^4)/(16*a^5*S1^4) + (20545*e^4*R^4*(S1 - S3)^3*Sin[i]^4)/(64*a^5*S1^4)) + C44*((-5565*e^2*R^4*Cos[i])/(4*a^5*S1) - (147945*e^4*R^4*Cos[i])/(32*a^5*S1) + (315*R^4*(S1 - S3)*Cos[i]^2)/(2*a^5*S1^2) - (1995*e^2*R^4*(S1 - S3)*Cos[i]^2)/(a^5*S1^2) - (402745*e^4*R^4*(S1 - S3)*Cos[i]^2)/(32*a^5*S1^2) + (5565*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]^3)/(4*a^5*S1) - (6265*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]^3)/(8*a^5*S1) + (5565*e^2*R^4*(S1 - S3)^2*Cos[i]^3)/(4*a^5*S1^3) - (544075*e^4*R^4*(S1 - S3)^2*Cos[i]^3)/(32*a^5*S1^3) - (315*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^4)/(2*a^5*S1^2) + (1995*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^4)/(a^5*S1^2) - (58135*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) + (1995*e^2*R^4*(S1 - S3)^3*Cos[i]^4)/(a^5*S1^4) - (289275*e^4*R^4*(S1 - S3)^3*Cos[i]^4)/(32*a^5*S1^4) + (5565*e^2*R^4*(S1 - S3)*Sin[i])/(4*a^5*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) + (147945*e^4*R^4*(S1 - S3)*Sin[i])/(32*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) - (1995*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(a^5*S1) - (402745*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(32*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]) + (1995*e^2*R^4*(S1 - S3)^2*Cos[i]*Sin[i])/(a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) + (402745*e^4*R^4*(S1 - S3)^2*Cos[i]*Sin[i])/(32*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) - (5565*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^2*Sin[i])/(4*a^5*S1^2) - (506485*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^2*Sin[i])/(16*a^5*S1^2) - (5565*e^2*R^4*(S1 - S3)^3*Cos[i]^2*Sin[i])/(4*a^5*S1^4*Sqrt[1 - (S1 - S3)^2/S1^2]) + (544075*e^4*R^4*(S1 - S3)^3*Cos[i]^2*Sin[i])/(32*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) + (1995*e^2*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(a^5*S1) - (58135*e^4*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(16*a^5*S1) - (519015*e^4*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])/(2*a^5*S1^5*Sqrt[1 - (S1 - S3)^2/S1^2]) - (1995*e^2*R^4*(S1 - S3)^4*Cos[i]^3*Sin[i])/(a^5*S1^5*Sqrt[1 - (S1 - S3)^2/S1^2]) + (289275*e^4*R^4*(S1 - S3)^4*Cos[i]^3*Sin[i])/(32*a^5*S1^5*Sqrt[1 - (S1 - S3)^2/S1^2]) - (315*R^4*(S1 - S3)*Sin[i]^2)/(2*a^5*S1^2) + (1995*e^2*R^4*(S1 - S3)*Sin[i]^2)/(a^5*S1^2) + (402745*e^4*R^4*(S1 - S3)*Sin[i]^2)/(32*a^5*S1^2) + (5565*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]*Sin[i]^2)/(4*a^5*S1) - (544075*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]*Sin[i]^2)/(32*a^5*S1) + (5565*e^2*R^4*(S1 - S3)^2*Cos[i]*Sin[i]^2)/(4*a^5*S1^3) + (506485*e^4*R^4*(S1 - S3)^2*Cos[i]*Sin[i]^2)/(16*a^5*S1^3) - (519015*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^2*Sin[i]^2)/(32*a^5*S1^2) + (519015*e^4*R^4*(S1 - S3)^3*Cos[i]^2*Sin[i]^2)/(32*a^5*S1^4) - (5565*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Sin[i]^3)/(4*a^5*S1^2) + (544075*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Sin[i]^3)/(32*a^5*S1^2) - (5565*e^2*R^4*(S1 - S3)^3*Sin[i]^3)/(4*a^5*S1^4*Sqrt[1 - (S1 - S3)^2/S1^2]) + (6265*e^4*R^4*(S1 - S3)^3*Sin[i]^3)/(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]*Sin[i]^3)/(2*a^5*S1) + (1995*e^2*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]*Sin[i]^3)/(a^5*S1) - (289275*e^4*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]*Sin[i]^3)/(32*a^5*S1) + (519015*e^4*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)/(2*a^5*S1^5*Sqrt[1 - (S1 - S3)^2/S1^2]) - (1995*e^2*R^4*(S1 - S3)^4*Cos[i]*Sin[i]^3)/(a^5*S1^5*Sqrt[1 - (S1 - S3)^2/S1^2]) + (58135*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) - (1995*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Sin[i]^4)/(a^5*S1^2) + (289275*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Sin[i]^4)/(32*a^5*S1^2) + (315*R^4*(S1 - S3)^3*Sin[i]^4)/(2*a^5*S1^4) - (1995*e^2*R^4*(S1 - S3)^3*Sin[i]^4)/(a^5*S1^4) + (58135*e^4*R^4*(S1 - S3)^3*Sin[i]^4)/(16*a^5*S1^4))) == 0}