sys31 = {-3*n - omd - Omd + S1/C - G*m*M*(C11*((27*e^2*R*Cos[i])/(16*a^2*S1) - (27*e^4*R*Cos[i])/(16*a^2*S1) - (27*e^2*R*(S1 - S3)*Cos[i])/(16*a^2*S1^2) + (27*e^4*R*(S1 - S3)*Cos[i])/(16*a^2*S1^2) + (27*e^2*R*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i])/(32*a^2*Sqrt[1 - (S1 - S3)^2/S1^2]) - (27*e^4*R*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i])/(32*a^2*Sqrt[1 - (S1 - S3)^2/S1^2])) + C22*((1599*e^4*R^2*Cos[i])/(32*a^3*S1) - (1599*e^4*R^2*(S1 - S3)*Cos[i])/(32*a^3*S1^2) + (1599*e^4*R^2*(S1 - S3)*Cos[i]^2)/(32*a^3*S1^2) - (1599*e^4*R^2*(S1 - S3)^2*Cos[i]^2)/(32*a^3*S1^3) + (1599*e^4*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]) + (1599*e^4*R^2*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(32*a^3*S1) + (1599*e^4*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]) - (1599*e^4*R^2*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i])/(32*a^3*S1^2) + (1599*e^4*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)) + C31*((-159*e^2*R^3*Cos[i])/(32*a^4*S1) - (117*e^4*R^3*Cos[i])/(64*a^4*S1) + (159*e^2*R^3*(S1 - S3)*Cos[i])/(32*a^4*S1^2) + (117*e^4*R^3*(S1 - S3)*Cos[i])/(64*a^4*S1^2) + (795*e^2*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Cos[i]^2)/(128*a^4) + (585*e^4*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Cos[i]^2)/(256*a^4) + (2385*e^2*R^3*(1 - (S1 - S3)^2/S1^2)*Cos[i]^3)/(128*a^4*S1) + (1755*e^4*R^3*(1 - (S1 - S3)^2/S1^2)*Cos[i]^3)/(256*a^4*S1) + (2385*e^2*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]^3)/(128*a^4*S1) + (1755*e^4*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]^3)/(256*a^4*S1) - (2385*e^2*R^3*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^3)/(128*a^4*S1^2) - (1755*e^4*R^3*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^3)/(256*a^4*S1^2) - (159*e^2*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]) - (117*e^4*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i])/(128*a^4*Sqrt[1 - (S1 - S3)^2/S1^2]) - (795*e^2*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(64*a^4*S1) - (585*e^4*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(128*a^4*S1) - (795*e^2*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i])/(128*a^4*S1*Sqrt[1 - (S1 - S3)^2/S1^2]) - (585*e^4*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i])/(256*a^4*S1*Sqrt[1 - (S1 - S3)^2/S1^2]) + (795*e^2*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i])/(64*a^4*S1^2) + (585*e^4*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i])/(128*a^4*S1^2) + (7155*e^2*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])/(256*a^4) + (5265*e^4*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])/(512*a^4) - (2385*e^2*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^2*Sin[i])/(32*a^4*S1^2) - (1755*e^4*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^2*Sin[i])/(64*a^4*S1^2) - (2385*e^2*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^2*Sin[i])/(128*a^4*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) - (1755*e^4*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^2*Sin[i])/(256*a^4*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) + (2385*e^2*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^2*Sin[i])/(32*a^4*S1^3) + (1755*e^4*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^2*Sin[i])/(64*a^4*S1^3) + (795*e^2*R^3*(S1 - S3)*Sin[i]^2)/(64*a^4*S1^2) + (585*e^4*R^3*(S1 - S3)*Sin[i]^2)/(128*a^4*S1^2) - (795*e^2*R^3*(S1 - S3)^2*Sin[i]^2)/(64*a^4*S1^3) - (585*e^4*R^3*(S1 - S3)^2*Sin[i]^2)/(128*a^4*S1^3) - (2385*e^2*R^3*(1 - (S1 - S3)^2/S1^2)*Cos[i]*Sin[i]^2)/(64*a^4*S1) - (1755*e^4*R^3*(1 - (S1 - S3)^2/S1^2)*Cos[i]*Sin[i]^2)/(128*a^4*S1) - (2385*e^2*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i]^2)/(64*a^4*S1) - (1755*e^4*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i]^2)/(128*a^4*S1) + (2385*e^2*R^3*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]*Sin[i]^2)/(64*a^4*S1^2) + (1755*e^4*R^3*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]*Sin[i]^2)/(128*a^4*S1^2) + (7155*e^2*R^3*(S1 - S3)^2*Cos[i]*Sin[i]^2)/(128*a^4*S1^3) + (5265*e^4*R^3*(S1 - S3)^2*Cos[i]*Sin[i]^2)/(256*a^4*S1^3) - (7155*e^2*R^3*(S1 - S3)^3*Cos[i]*Sin[i]^2)/(128*a^4*S1^4) - (5265*e^4*R^3*(S1 - S3)^3*Cos[i]*Sin[i]^2)/(256*a^4*S1^4) + (2385*e^2*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Sin[i]^3)/(64*a^4*S1^2) + (1755*e^4*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Sin[i]^3)/(128*a^4*S1^2) + (2385*e^2*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Sin[i]^3)/(256*a^4*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) + (1755*e^4*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Sin[i]^3)/(512*a^4*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) - (2385*e^2*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Sin[i]^3)/(64*a^4*S1^3) - (1755*e^4*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Sin[i]^3)/(128*a^4*S1^3)) + C44*((58275*e^4*R^4*Cos[i])/(8*a^5*S1) - (58275*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) - (2625*e^2*R^4*(S1 - S3)*Cos[i]^2)/(2*a^5*S1^2) + (177975*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) + (2625*e^2*R^4*(S1 - S3)^2*Cos[i]^2)/(2*a^5*S1^3) - (177975*e^4*R^4*(S1 - S3)^2*Cos[i]^2)/(16*a^5*S1^3) - (58275*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]^3)/(8*a^5*S1) - (58275*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]^3)/(8*a^5*S1) + (58275*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^3)/(8*a^5*S1^2) - (174825*e^4*R^4*(S1 - S3)^2*Cos[i]^3)/(8*a^5*S1^3) + (174825*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) + (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^4)/(a^5*S1^2) - (83475*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) + (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) - (83475*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) - (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Cos[i]^4)/(a^5*S1^3) + (83475*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) + (4725*e^2*R^4*(S1 - S3)^3*Cos[i]^4)/(2*a^5*S1^4) - (366975*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) - (4725*e^2*R^4*(S1 - S3)^4*Cos[i]^4)/(2*a^5*S1^5) + (366975*e^4*R^4*(S1 - S3)^4*Cos[i]^4)/(16*a^5*S1^5) + (58275*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) - (1575*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(2*a^5*S1) + (200025*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]) - (1575*e^2*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]) + (200025*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) + (1575*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i])/(2*a^5*S1^2) - (200025*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i])/(16*a^5*S1^2) - (174825*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) - (58275*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^2*Sin[i])/(4*a^5*S1^2) - (58275*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]) + (58275*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) + (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(2*a^5*S1) - (200025*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) + (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) - (600075*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) - (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) + (200025*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) + (4725*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^3*Sin[i])/(2*a^5*S1^3) - (600075*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]) + (1575*e^2*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]) - (200025*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) - (4725*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^3*Cos[i]^3*Sin[i])/(2*a^5*S1^4) + (600075*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) - (2625*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i]^2)/(4*a^5) + (177975*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) - (58275*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]*Sin[i]^2)/(8*a^5*S1) - (58275*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) + (58275*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]*Sin[i]^2)/(8*a^5*S1^2) - (174825*e^4*R^4*(S1 - S3)^2*Cos[i]*Sin[i]^2)/(8*a^5*S1^3) + (174825*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 + (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 - (83475*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) + (3150*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^2*Sin[i]^2)/(a^5*S1^2) - (83475*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^2*Sin[i]^2)/(4*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) + (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)/(a^5*S1^2) - (83475*e^4*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) - (630*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Cos[i]^2*Sin[i]^2)/(a^5*S1^3) - (3150*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Cos[i]^2*Sin[i]^2)/(a^5*S1^3) + (83475*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Cos[i]^2*Sin[i]^2)/(4*a^5*S1^3) + (630*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) - (83475*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) - (3150*e^2*R^4*(S1 - S3)^4*Cos[i]^2*Sin[i]^2)/(a^5*S1^5) + (83475*e^4*R^4*(S1 - S3)^4*Cos[i]^2*Sin[i]^2)/(4*a^5*S1^5) - (174825*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) - (58275*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Sin[i]^3)/(4*a^5*S1^2) - (58275*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]) + (58275*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) + (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]*Sin[i]^3)/(2*a^5*S1) - (200025*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) + (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)/(4*a^5*S1) - (600075*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) - (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*(S1 - S3)*Cos[i]*Sin[i]^3)/(2*a^5*S1^2) + (200025*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) + (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) - (600075*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]) + (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]) - (200025*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) - (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) + (600075*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) + (4725*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) - (366975*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) + (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Sin[i]^4)/(a^5*S1^2) - (83475*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) + (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) - (83475*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) - (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Sin[i]^4)/(a^5*S1^3) + (83475*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*((-8325*e^4*R^4*Cos[i])/(32*a^5*S1) + (8325*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) + (75*e^2*R^4*(S1 - S3)*Cos[i]^2)/(2*a^5*S1^2) - (5175*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) - (75*e^2*R^4*(S1 - S3)^2*Cos[i]^2)/(2*a^5*S1^3) + (5175*e^4*R^4*(S1 - S3)^2*Cos[i]^2)/(32*a^5*S1^3) + (58275*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]^3)/(64*a^5*S1) + (58275*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]^3)/(64*a^5*S1) - (58275*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) - (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^4)/(8*a^5*S1^2) + (83475*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) - (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) + (83475*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) + (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Cos[i]^4)/(8*a^5*S1^3) - (83475*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Cos[i]^4)/(64*a^5*S1^3) - (8325*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) - (1575*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(8*a^5*S1) - (49725*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]) - (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]) - (49725*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) + (1575*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i])/(8*a^5*S1^2) + (49725*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i])/(64*a^5*S1^2) + (174825*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) - (58275*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^2*Sin[i])/(16*a^5*S1^2) - (58275*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]) + (58275*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) + (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(8*a^5*S1) + (149625*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) + (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) + (448875*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) - (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) - (149625*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) + (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) - (250425*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]) + (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]) - (83475*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) - (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) + (250425*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) + (75*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i]^2)/(4*a^5) - (5175*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) - (58275*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]*Sin[i]^2)/(32*a^5*S1) - (58275*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) + (58275*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]*Sin[i]^2)/(32*a^5*S1^2) + (174825*e^4*R^4*(S1 - S3)^2*Cos[i]*Sin[i]^2)/(64*a^5*S1^3) - (174825*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) - (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) + (83475*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) - (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^2*Sin[i]^2)/(4*a^5*S1^2) - (266175*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) - (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) - (266175*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) + (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) + (266175*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) - (1575*e^2*R^4*(S1 - S3)^3*Cos[i]^2*Sin[i]^2)/(4*a^5*S1^4) + (83475*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) + (1575*e^2*R^4*(S1 - S3)^4*Cos[i]^2*Sin[i]^2)/(4*a^5*S1^5) - (83475*e^4*R^4*(S1 - S3)^4*Cos[i]^2*Sin[i]^2)/(32*a^5*S1^5) + (58275*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Sin[i]^3)/(32*a^5*S1^2) + (58275*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]) - (58275*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) + (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]*Sin[i]^3)/(8*a^5*S1) - (83475*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) + (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) - (250425*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) - (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) + (83475*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) + (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) + (448875*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]) + (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]) + (149625*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) - (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) - (448875*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) - (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Sin[i]^4)/(8*a^5*S1^2) + (83475*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) - (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) + (83475*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) + (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Sin[i]^4)/(8*a^5*S1^3) - (83475*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*((-27*e^2*R*Cos[i])/(16*a^2*S1) + (27*e^4*R*Cos[i])/(16*a^2*S1) + (27*e^2*R*(S1 - S3)*Sin[i])/(16*a^2*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) - (27*e^4*R*(S1 - S3)*Sin[i])/(16*a^2*S1^2*Sqrt[1 - (S1 - S3)^2/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)) + C22*((-1599*e^4*R^2*Cos[i])/(32*a^3*S1) - (1599*e^4*R^2*(S1 - S3)*Cos[i]^2)/(32*a^3*S1^2) + (1599*e^4*R^2*(S1 - S3)*Sin[i])/(32*a^3*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) - (1599*e^4*R^2*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(32*a^3*S1) + (1599*e^4*R^2*(S1 - S3)^2*Cos[i]*Sin[i])/(32*a^3*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) + (1599*e^4*R^2*(S1 - S3)*Sin[i]^2)/(32*a^3*S1^2)) + C31*((159*e^2*R^3*Cos[i])/(32*a^4*S1) + (117*e^4*R^3*Cos[i])/(64*a^4*S1) + (795*e^2*R^3*(S1 - S3)*Cos[i]^2)/(64*a^4*S1^2) + (585*e^4*R^3*(S1 - S3)*Cos[i]^2)/(128*a^4*S1^2) - (2385*e^2*R^3*(1 - (S1 - S3)^2/S1^2)*Cos[i]^3)/(128*a^4*S1) - (1755*e^4*R^3*(1 - (S1 - S3)^2/S1^2)*Cos[i]^3)/(256*a^4*S1) + (2385*e^2*R^3*(S1 - S3)^2*Cos[i]^3)/(64*a^4*S1^3) + (1755*e^4*R^3*(S1 - S3)^2*Cos[i]^3)/(128*a^4*S1^3) - (159*e^2*R^3*(S1 - S3)*Sin[i])/(32*a^4*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) - (117*e^4*R^3*(S1 - S3)*Sin[i])/(64*a^4*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) + (795*e^2*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(64*a^4*S1) + (585*e^4*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(128*a^4*S1) - (795*e^2*R^3*(S1 - S3)^2*Cos[i]*Sin[i])/(64*a^4*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) - (585*e^4*R^3*(S1 - S3)^2*Cos[i]*Sin[i])/(128*a^4*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) + (16695*e^2*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^2*Sin[i])/(128*a^4*S1^2) + (12285*e^4*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^2*Sin[i])/(256*a^4*S1^2) - (2385*e^2*R^3*(S1 - S3)^3*Cos[i]^2*Sin[i])/(64*a^4*S1^4*Sqrt[1 - (S1 - S3)^2/S1^2]) - (1755*e^4*R^3*(S1 - S3)^3*Cos[i]^2*Sin[i])/(128*a^4*S1^4*Sqrt[1 - (S1 - S3)^2/S1^2]) - (795*e^2*R^3*(S1 - S3)*Sin[i]^2)/(64*a^4*S1^2) - (585*e^4*R^3*(S1 - S3)*Sin[i]^2)/(128*a^4*S1^2) + (2385*e^2*R^3*(1 - (S1 - S3)^2/S1^2)*Cos[i]*Sin[i]^2)/(64*a^4*S1) + (1755*e^4*R^3*(1 - (S1 - S3)^2/S1^2)*Cos[i]*Sin[i]^2)/(128*a^4*S1) - (16695*e^2*R^3*(S1 - S3)^2*Cos[i]*Sin[i]^2)/(128*a^4*S1^3) - (12285*e^4*R^3*(S1 - S3)^2*Cos[i]*Sin[i]^2)/(256*a^4*S1^3) - (2385*e^2*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Sin[i]^3)/(64*a^4*S1^2) - (1755*e^4*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Sin[i]^3)/(128*a^4*S1^2) + (2385*e^2*R^3*(S1 - S3)^3*Sin[i]^3)/(128*a^4*S1^4*Sqrt[1 - (S1 - S3)^2/S1^2]) + (1755*e^4*R^3*(S1 - S3)^3*Sin[i]^3)/(256*a^4*S1^4*Sqrt[1 - (S1 - S3)^2/S1^2])) + C44*((-58275*e^4*R^4*Cos[i])/(8*a^5*S1) + (315*R^4*(S1 - S3)*Cos[i]^2)/(2*a^5*S1^2) + (1575*e^2*R^4*(S1 - S3)*Cos[i]^2)/(2*a^5*S1^2) - (200025*e^4*R^4*(S1 - S3)*Cos[i]^2)/(16*a^5*S1^2) + (58275*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]^3)/(8*a^5*S1) + (58275*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) - (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^4)/(2*a^5*S1^2) + (200025*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) - (1575*e^2*R^4*(S1 - S3)^3*Cos[i]^4)/(2*a^5*S1^4) + (200025*e^4*R^4*(S1 - S3)^3*Cos[i]^4)/(16*a^5*S1^4) + (58275*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) + (1575*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(2*a^5*S1) - (200025*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]) - (1575*e^2*R^4*(S1 - S3)^2*Cos[i]*Sin[i])/(2*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) + (200025*e^4*R^4*(S1 - S3)^2*Cos[i]*Sin[i])/(16*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) - (58275*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^2*Sin[i])/(8*a^5*S1^2) - (58275*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) - (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(2*a^5*S1) + (200025*e^4*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(16*a^5*S1) + (315*R^4*(S1 - S3)^4*Cos[i]^3*Sin[i])/(2*a^5*S1^5*Sqrt[1 - (S1 - S3)^2/S1^2]) + (1575*e^2*R^4*(S1 - S3)^4*Cos[i]^3*Sin[i])/(2*a^5*S1^5*Sqrt[1 - (S1 - S3)^2/S1^2]) - (200025*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) - (1575*e^2*R^4*(S1 - S3)*Sin[i]^2)/(2*a^5*S1^2) + (200025*e^4*R^4*(S1 - S3)*Sin[i]^2)/(16*a^5*S1^2) + (58275*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]*Sin[i]^2)/(8*a^5*S1) + (58275*e^4*R^4*(S1 - S3)^2*Cos[i]*Sin[i]^2)/(8*a^5*S1^3) - (58275*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Sin[i]^3)/(8*a^5*S1^2) - (58275*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) - (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]*Sin[i]^3)/(2*a^5*S1) + (200025*e^4*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]*Sin[i]^3)/(16*a^5*S1) + (315*R^4*(S1 - S3)^4*Cos[i]*Sin[i]^3)/(2*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]) - (200025*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) + (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Sin[i]^4)/(2*a^5*S1^2) - (200025*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) + (1575*e^2*R^4*(S1 - S3)^3*Sin[i]^4)/(2*a^5*S1^4) - (200025*e^4*R^4*(S1 - S3)^3*Sin[i]^4)/(16*a^5*S1^4)) + 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*((8325*e^4*R^4*Cos[i])/(32*a^5*S1) + (315*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) + (49725*e^4*R^4*(S1 - S3)*Cos[i]^2)/(64*a^5*S1^2) - (58275*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]^3)/(64*a^5*S1) + (58275*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) - (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^4)/(8*a^5*S1^2) - (149625*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) - (1575*e^2*R^4*(S1 - S3)^3*Cos[i]^4)/(8*a^5*S1^4) + (83475*e^4*R^4*(S1 - S3)^3*Cos[i]^4)/(64*a^5*S1^4) - (8325*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) + (1575*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(8*a^5*S1) + (49725*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]) - (1575*e^2*R^4*(S1 - S3)^2*Cos[i]*Sin[i])/(8*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) - (49725*e^4*R^4*(S1 - S3)^2*Cos[i]*Sin[i])/(64*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) + (407925*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^2*Sin[i])/(64*a^5*S1^2) - (58275*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) - (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(8*a^5*S1) - (149625*e^4*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(64*a^5*S1) + (174825*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]) + (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]) - (83475*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) - (1575*e^2*R^4*(S1 - S3)*Sin[i]^2)/(8*a^5*S1^2) - (49725*e^4*R^4*(S1 - S3)*Sin[i]^2)/(64*a^5*S1^2) + (58275*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]*Sin[i]^2)/(32*a^5*S1) - (407925*e^4*R^4*(S1 - S3)^2*Cos[i]*Sin[i]^2)/(64*a^5*S1^3) + (174825*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^2*Sin[i]^2)/(16*a^5*S1^2) - (174825*e^4*R^4*(S1 - S3)^3*Cos[i]^2*Sin[i]^2)/(16*a^5*S1^4) - (58275*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Sin[i]^3)/(32*a^5*S1^2) + (58275*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) - (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]*Sin[i]^3)/(8*a^5*S1) + (83475*e^4*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]*Sin[i]^3)/(64*a^5*S1) - (174825*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]) + (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]) + (149625*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) + (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Sin[i]^4)/(8*a^5*S1^2) - (83475*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) + (1575*e^2*R^4*(S1 - S3)^3*Sin[i]^4)/(8*a^5*S1^4) + (149625*e^4*R^4*(S1 - S3)^3*Sin[i]^4)/(64*a^5*S1^4))) == 0}