sys52 = {(-5*n)/2 - omd - Omd + S1/C - G*m*M*(C22*((845*e^3*R^2*Cos[i])/(32*a^3*S1) - (845*e^3*R^2*(S1 - S3)*Cos[i])/(32*a^3*S1^2) + (845*e^3*R^2*(S1 - S3)*Cos[i]^2)/(32*a^3*S1^2) - (845*e^3*R^2*(S1 - S3)^2*Cos[i]^2)/(32*a^3*S1^3) + (845*e^3*R^2*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i])/(64*a^3*Sqrt[1 - (S1 - S3)^2/S1^2]) + (845*e^3*R^2*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(32*a^3*S1) + (845*e^3*R^2*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i])/(64*a^3*S1*Sqrt[1 - (S1 - S3)^2/S1^2]) - (845*e^3*R^2*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i])/(32*a^3*S1^2) + (845*e^3*R^2*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i]^2)/(64*a^3)) + C20*((3*R^2*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Cos[i]^2)/(4*a^3) + (9*e^2*R^2*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Cos[i]^2)/(8*a^3) + (45*e^4*R^2*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Cos[i]^2)/(32*a^3) - (3*R^2*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(2*a^3*S1) - (9*e^2*R^2*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(4*a^3*S1) - (45*e^4*R^2*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(16*a^3*S1) - (3*R^2*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i])/(4*a^3*S1*Sqrt[1 - (S1 - S3)^2/S1^2]) - (9*e^2*R^2*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i])/(8*a^3*S1*Sqrt[1 - (S1 - S3)^2/S1^2]) - (45*e^4*R^2*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i])/(32*a^3*S1*Sqrt[1 - (S1 - S3)^2/S1^2]) + (3*R^2*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i])/(2*a^3*S1^2) + (9*e^2*R^2*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i])/(4*a^3*S1^2) + (45*e^4*R^2*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i])/(16*a^3*S1^2) + (3*R^2*(S1 - S3)*Sin[i]^2)/(2*a^3*S1^2) + (9*e^2*R^2*(S1 - S3)*Sin[i]^2)/(4*a^3*S1^2) + (45*e^4*R^2*(S1 - S3)*Sin[i]^2)/(16*a^3*S1^2) - (3*R^2*(S1 - S3)^2*Sin[i]^2)/(2*a^3*S1^3) - (9*e^2*R^2*(S1 - S3)^2*Sin[i]^2)/(4*a^3*S1^3) - (45*e^4*R^2*(S1 - S3)^2*Sin[i]^2)/(16*a^3*S1^3)) + C44*((53935*e^3*R^4*Cos[i])/(16*a^5*S1) - (53935*e^3*R^4*(S1 - S3)*Cos[i])/(16*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) + (53935*e^3*R^4*(S1 - S3)*Cos[i]^2)/(8*a^5*S1^2) - (55125*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) - (53935*e^3*R^4*(S1 - S3)^2*Cos[i]^2)/(8*a^5*S1^3) + (55125*e^4*R^4*(S1 - S3)^2*Cos[i]^2)/(16*a^5*S1^3) - (53935*e^3*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]^3)/(16*a^5*S1) - (53935*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]^3)/(16*a^5*S1) + (53935*e^3*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^3)/(16*a^5*S1^2) - (161805*e^3*R^4*(S1 - S3)^2*Cos[i]^3)/(16*a^5*S1^3) + (161805*e^3*R^4*(S1 - S3)^3*Cos[i]^3)/(16*a^5*S1^4) + (315*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(1 - (S1 - S3)^2/S1^2)*Cos[i]^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) - (53935*e^3*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^4)/(8*a^5*S1^2) + (33075*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^4)/(8*a^5*S1^2) + (315*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^4)/(2*a^5*S1^2) + (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) - (53935*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^4)/(16*a^5*S1^2) + (33075*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^4)/(16*a^5*S1^2) - (315*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Cos[i]^4)/(a^5*S1^3) - (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Cos[i]^4)/(a^5*S1^3) + (53935*e^3*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Cos[i]^4)/(8*a^5*S1^3) - (33075*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Cos[i]^4)/(8*a^5*S1^3) + (945*R^4*(S1 - S3)^3*Cos[i]^4)/(2*a^5*S1^4) + (4725*e^2*R^4*(S1 - S3)^3*Cos[i]^4)/(2*a^5*S1^4) - (53935*e^3*R^4*(S1 - S3)^3*Cos[i]^4)/(4*a^5*S1^4) + (99225*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) + (53935*e^3*R^4*(S1 - S3)^4*Cos[i]^4)/(4*a^5*S1^5) - (99225*e^4*R^4*(S1 - S3)^4*Cos[i]^4)/(16*a^5*S1^5) + (53935*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i])/(32*a^5*Sqrt[1 - (S1 - S3)^2/S1^2]) - (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) + (53935*e^3*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(8*a^5*S1) - (33075*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]) + (53935*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i])/(16*a^5*S1*Sqrt[1 - (S1 - S3)^2/S1^2]) - (33075*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) - (53935*e^3*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i])/(8*a^5*S1^2) + (33075*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i])/(16*a^5*S1^2) - (161805*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]^2*Sin[i])/(32*a^5) - (53935*e^3*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^2*Sin[i])/(8*a^5*S1^2) - (53935*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^2*Sin[i])/(32*a^5*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) + (53935*e^3*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^2*Sin[i])/(8*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) - (53935*e^3*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(8*a^5*S1) + (33075*e^4*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(16*a^5*S1) + (945*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^3*Sin[i])/(4*a^5*S1) + (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) - (161805*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^3*Sin[i])/(16*a^5*S1) + (99225*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^3*Sin[i])/(32*a^5*S1) - (315*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*(S1 - S3)*Cos[i]^3*Sin[i])/(2*a^5*S1^2) - (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) + (53935*e^3*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*(S1 - S3)*Cos[i]^3*Sin[i])/(8*a^5*S1^2) - (33075*e^4*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*(S1 - S3)*Cos[i]^3*Sin[i])/(16*a^5*S1^2) + (945*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^3*Sin[i])/(2*a^5*S1^3) + (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) - (161805*e^3*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^3*Sin[i])/(8*a^5*S1^3) + (99225*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]) - (53935*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^3*Cos[i]^3*Sin[i])/(16*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) + (33075*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) + (161805*e^3*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^3*Cos[i]^3*Sin[i])/(8*a^5*S1^4) - (99225*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) + (53935*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i]^2)/(16*a^5) - (55125*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) - (53935*e^3*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]*Sin[i]^2)/(16*a^5*S1) - (53935*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i]^2)/(16*a^5*S1) + (53935*e^3*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]*Sin[i]^2)/(16*a^5*S1^2) - (161805*e^3*R^4*(S1 - S3)^2*Cos[i]*Sin[i]^2)/(16*a^5*S1^3) + (161805*e^3*R^4*(S1 - S3)^3*Cos[i]*Sin[i]^2)/(16*a^5*S1^4) + (315*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(1 - (S1 - S3)^2/S1^2)*Cos[i]^2*Sin[i]^2)/a^5 + (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 - (53935*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(1 - (S1 - S3)^2/S1^2)*Cos[i]^2*Sin[i]^2)/(8*a^5) + (33075*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(1 - (S1 - S3)^2/S1^2)*Cos[i]^2*Sin[i]^2)/(8*a^5) + (630*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^2*Sin[i]^2)/(a^5*S1^2) + (3150*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^2*Sin[i]^2)/(a^5*S1^2) - (53935*e^3*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^2*Sin[i]^2)/(4*a^5*S1^2) + (33075*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) - (53935*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^2*Sin[i]^2)/(8*a^5*S1^2) + (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)/(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) + (53935*e^3*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Cos[i]^2*Sin[i]^2)/(4*a^5*S1^3) - (33075*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) - (53935*e^3*R^4*(S1 - S3)^3*Cos[i]^2*Sin[i]^2)/(4*a^5*S1^4) + (33075*e^4*R^4*(S1 - S3)^3*Cos[i]^2*Sin[i]^2)/(4*a^5*S1^4) - (630*R^4*(S1 - S3)^4*Cos[i]^2*Sin[i]^2)/(a^5*S1^5) - (3150*e^2*R^4*(S1 - S3)^4*Cos[i]^2*Sin[i]^2)/(a^5*S1^5) + (53935*e^3*R^4*(S1 - S3)^4*Cos[i]^2*Sin[i]^2)/(4*a^5*S1^5) - (33075*e^4*R^4*(S1 - S3)^4*Cos[i]^2*Sin[i]^2)/(4*a^5*S1^5) - (161805*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sqrt[1 - (S1 - S3)^2/S1^2]*Sin[i]^3)/(32*a^5) - (53935*e^3*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Sin[i]^3)/(8*a^5*S1^2) - (53935*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Sin[i]^3)/(32*a^5*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) + (53935*e^3*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Sin[i]^3)/(8*a^5*S1^3) + (315*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]*Sin[i]^3)/(2*a^5*S1) + (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]*Sin[i]^3)/(2*a^5*S1) - (53935*e^3*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]*Sin[i]^3)/(8*a^5*S1) + (33075*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) - (161805*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i]^3)/(16*a^5*S1) + (99225*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i]^3)/(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) + (53935*e^3*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*(S1 - S3)*Cos[i]*Sin[i]^3)/(8*a^5*S1^2) - (33075*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) - (161805*e^3*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]*Sin[i]^3)/(8*a^5*S1^3) + (99225*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]*Sin[i]^3)/(16*a^5*S1^3) + (315*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^3*Cos[i]*Sin[i]^3)/(4*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) + (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]) - (53935*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^3*Cos[i]*Sin[i]^3)/(16*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) + (33075*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^3*Cos[i]*Sin[i]^3)/(32*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) - (945*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^3*Cos[i]*Sin[i]^3)/(2*a^5*S1^4) - (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) + (161805*e^3*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^3*Cos[i]*Sin[i]^3)/(8*a^5*S1^4) - (99225*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^3*Cos[i]*Sin[i]^3)/(16*a^5*S1^4) + (945*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(1 - (S1 - S3)^2/S1^2)*Sin[i]^4)/(4*a^5) + (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) - (53935*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(1 - (S1 - S3)^2/S1^2)*Sin[i]^4)/(8*a^5) + (99225*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) - (53935*e^3*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Sin[i]^4)/(8*a^5*S1^2) + (33075*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Sin[i]^4)/(8*a^5*S1^2) + (315*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Sin[i]^4)/(2*a^5*S1^2) + (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) - (53935*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Sin[i]^4)/(16*a^5*S1^2) + (33075*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Sin[i]^4)/(16*a^5*S1^2) - (315*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Sin[i]^4)/(a^5*S1^3) - (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Sin[i]^4)/(a^5*S1^3) + (53935*e^3*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Sin[i]^4)/(8*a^5*S1^3) - (33075*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Sin[i]^4)/(8*a^5*S1^3) + (315*R^4*(S1 - S3)^3*Sin[i]^4)/(2*a^5*S1^4) + (1575*e^2*R^4*(S1 - S3)^3*Sin[i]^4)/(2*a^5*S1^4) + (33075*e^4*R^4*(S1 - S3)^3*Sin[i]^4)/(16*a^5*S1^4) - (315*R^4*(S1 - S3)^4*Sin[i]^4)/(2*a^5*S1^5) - (1575*e^2*R^4*(S1 - S3)^4*Sin[i]^4)/(2*a^5*S1^5) - (33075*e^4*R^4*(S1 - S3)^4*Sin[i]^4)/(16*a^5*S1^5)) + C40*((-15*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Cos[i]^2)/(8*a^5) - (75*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Cos[i]^2)/(8*a^5) - (1575*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Cos[i]^2)/(64*a^5) + (105*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(1 - (S1 - S3)^2/S1^2)*Cos[i]^4)/(32*a^5) + (525*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(1 - (S1 - S3)^2/S1^2)*Cos[i]^4)/(32*a^5) + (11025*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(1 - (S1 - S3)^2/S1^2)*Cos[i]^4)/(256*a^5) + (15*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(4*a^5*S1) + (75*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(4*a^5*S1) + (1575*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(32*a^5*S1) + (15*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i])/(8*a^5*S1*Sqrt[1 - (S1 - S3)^2/S1^2]) + (75*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i])/(8*a^5*S1*Sqrt[1 - (S1 - S3)^2/S1^2]) + (1575*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i])/(64*a^5*S1*Sqrt[1 - (S1 - S3)^2/S1^2]) - (15*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i])/(4*a^5*S1^2) - (75*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i])/(4*a^5*S1^2) - (1575*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i])/(32*a^5*S1^2) - (105*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(16*a^5*S1) - (525*e^2*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(16*a^5*S1) - (11025*e^4*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(128*a^5*S1) - (315*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^3*Sin[i])/(32*a^5*S1) - (1575*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^3*Sin[i])/(32*a^5*S1) - (33075*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^3*Sin[i])/(256*a^5*S1) + (105*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*(S1 - S3)*Cos[i]^3*Sin[i])/(16*a^5*S1^2) + (525*e^2*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*(S1 - S3)*Cos[i]^3*Sin[i])/(16*a^5*S1^2) + (11025*e^4*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*(S1 - S3)*Cos[i]^3*Sin[i])/(128*a^5*S1^2) - (15*R^4*(S1 - S3)*Sin[i]^2)/(4*a^5*S1^2) - (75*e^2*R^4*(S1 - S3)*Sin[i]^2)/(4*a^5*S1^2) - (1575*e^4*R^4*(S1 - S3)*Sin[i]^2)/(32*a^5*S1^2) + (15*R^4*(S1 - S3)^2*Sin[i]^2)/(4*a^5*S1^3) + (75*e^2*R^4*(S1 - S3)^2*Sin[i]^2)/(4*a^5*S1^3) + (1575*e^4*R^4*(S1 - S3)^2*Sin[i]^2)/(32*a^5*S1^3) + (315*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^2*Sin[i]^2)/(16*a^5*S1^2) + (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^2*Sin[i]^2)/(16*a^5*S1^2) + (33075*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^2*Sin[i]^2)/(128*a^5*S1^2) + (315*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^2*Sin[i]^2)/(32*a^5*S1^2) + (1575*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^2*Sin[i]^2)/(32*a^5*S1^2) + (33075*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^2*Sin[i]^2)/(256*a^5*S1^2) - (315*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Cos[i]^2*Sin[i]^2)/(16*a^5*S1^3) - (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Cos[i]^2*Sin[i]^2)/(16*a^5*S1^3) - (33075*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Cos[i]^2*Sin[i]^2)/(128*a^5*S1^3) - (315*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]*Sin[i]^3)/(16*a^5*S1^3) - (1575*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]*Sin[i]^3)/(16*a^5*S1^3) - (33075*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]*Sin[i]^3)/(128*a^5*S1^3) - (105*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^3*Cos[i]*Sin[i]^3)/(32*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) - (525*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^3*Cos[i]*Sin[i]^3)/(32*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) - (11025*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^3*Cos[i]*Sin[i]^3)/(256*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) + (315*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^3*Cos[i]*Sin[i]^3)/(16*a^5*S1^4) + (1575*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^3*Cos[i]*Sin[i]^3)/(16*a^5*S1^4) + (33075*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^3*Cos[i]*Sin[i]^3)/(128*a^5*S1^4) + (105*R^4*(S1 - S3)^3*Sin[i]^4)/(16*a^5*S1^4) + (525*e^2*R^4*(S1 - S3)^3*Sin[i]^4)/(16*a^5*S1^4) + (11025*e^4*R^4*(S1 - S3)^3*Sin[i]^4)/(128*a^5*S1^4) - (105*R^4*(S1 - S3)^4*Sin[i]^4)/(16*a^5*S1^5) - (525*e^2*R^4*(S1 - S3)^4*Sin[i]^4)/(16*a^5*S1^5) - (11025*e^4*R^4*(S1 - S3)^4*Sin[i]^4)/(128*a^5*S1^5)) + C42*((-7705*e^3*R^4*Cos[i])/(64*a^5*S1) + (7705*e^3*R^4*(S1 - S3)*Cos[i])/(64*a^5*S1^2) + (375*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Cos[i]^2)/(16*a^5) + (1875*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Cos[i]^2)/(16*a^5) + (39375*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Cos[i]^2)/(128*a^5) + (15*R^4*(S1 - S3)*Cos[i]^2)/(2*a^5*S1^2) + (75*e^2*R^4*(S1 - S3)*Cos[i]^2)/(2*a^5*S1^2) - (7705*e^3*R^4*(S1 - S3)*Cos[i]^2)/(64*a^5*S1^2) + (1575*e^4*R^4*(S1 - S3)*Cos[i]^2)/(16*a^5*S1^2) - (15*R^4*(S1 - S3)^2*Cos[i]^2)/(2*a^5*S1^3) - (75*e^2*R^4*(S1 - S3)^2*Cos[i]^2)/(2*a^5*S1^3) + (7705*e^3*R^4*(S1 - S3)^2*Cos[i]^2)/(64*a^5*S1^3) - (1575*e^4*R^4*(S1 - S3)^2*Cos[i]^2)/(16*a^5*S1^3) + (53935*e^3*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]^3)/(128*a^5*S1) + (53935*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]^3)/(128*a^5*S1) - (53935*e^3*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^3)/(128*a^5*S1^2) - (315*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(1 - (S1 - S3)^2/S1^2)*Cos[i]^4)/(8*a^5) - (1575*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(1 - (S1 - S3)^2/S1^2)*Cos[i]^4)/(8*a^5) - (33075*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(1 - (S1 - S3)^2/S1^2)*Cos[i]^4)/(64*a^5) - (315*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^4)/(8*a^5*S1^2) - (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^4)/(8*a^5*S1^2) + (53935*e^3*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^4)/(64*a^5*S1^2) - (33075*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^4)/(64*a^5*S1^2) - (315*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^4)/(16*a^5*S1^2) - (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) + (53935*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^4)/(128*a^5*S1^2) - (33075*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^4)/(128*a^5*S1^2) + (315*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Cos[i]^4)/(8*a^5*S1^3) + (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Cos[i]^4)/(8*a^5*S1^3) - (53935*e^3*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Cos[i]^4)/(64*a^5*S1^3) + (33075*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Cos[i]^4)/(64*a^5*S1^3) - (7705*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i])/(128*a^5*Sqrt[1 - (S1 - S3)^2/S1^2]) - (315*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(8*a^5*S1) - (1575*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(8*a^5*S1) - (7705*e^3*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(64*a^5*S1) - (33075*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(64*a^5*S1) - (315*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i])/(16*a^5*S1*Sqrt[1 - (S1 - S3)^2/S1^2]) - (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]) - (7705*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i])/(128*a^5*S1*Sqrt[1 - (S1 - S3)^2/S1^2]) - (33075*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i])/(128*a^5*S1*Sqrt[1 - (S1 - S3)^2/S1^2]) + (315*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i])/(8*a^5*S1^2) + (1575*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i])/(8*a^5*S1^2) + (7705*e^3*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i])/(64*a^5*S1^2) + (33075*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i])/(64*a^5*S1^2) + (161805*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]^2*Sin[i])/(256*a^5) - (53935*e^3*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^2*Sin[i])/(32*a^5*S1^2) - (53935*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^2*Sin[i])/(128*a^5*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) + (53935*e^3*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^2*Sin[i])/(32*a^5*S1^3) + (315*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(8*a^5*S1) + (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(8*a^5*S1) + (53935*e^3*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(64*a^5*S1) + (33075*e^4*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(64*a^5*S1) + (945*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^3*Sin[i])/(16*a^5*S1) + (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) + (161805*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^3*Sin[i])/(128*a^5*S1) + (99225*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^3*Sin[i])/(128*a^5*S1) - (315*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*(S1 - S3)*Cos[i]^3*Sin[i])/(8*a^5*S1^2) - (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) - (53935*e^3*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*(S1 - S3)*Cos[i]^3*Sin[i])/(64*a^5*S1^2) - (33075*e^4*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*(S1 - S3)*Cos[i]^3*Sin[i])/(64*a^5*S1^2) + (945*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^3*Sin[i])/(8*a^5*S1^3) + (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) - (161805*e^3*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^3*Sin[i])/(64*a^5*S1^3) + (99225*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^3*Sin[i])/(64*a^5*S1^3) + (315*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^3*Cos[i]^3*Sin[i])/(16*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) + (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]) - (53935*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^3*Cos[i]^3*Sin[i])/(128*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) + (33075*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^3*Cos[i]^3*Sin[i])/(128*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) - (945*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^3*Cos[i]^3*Sin[i])/(8*a^5*S1^4) - (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) + (161805*e^3*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^3*Cos[i]^3*Sin[i])/(64*a^5*S1^4) - (99225*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^3*Cos[i]^3*Sin[i])/(64*a^5*S1^4) + (15*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i]^2)/(4*a^5) + (75*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i]^2)/(4*a^5) - (7705*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i]^2)/(128*a^5) + (1575*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i]^2)/(32*a^5) + (375*R^4*(S1 - S3)*Sin[i]^2)/(8*a^5*S1^2) + (1875*e^2*R^4*(S1 - S3)*Sin[i]^2)/(8*a^5*S1^2) + (39375*e^4*R^4*(S1 - S3)*Sin[i]^2)/(64*a^5*S1^2) - (375*R^4*(S1 - S3)^2*Sin[i]^2)/(8*a^5*S1^3) - (1875*e^2*R^4*(S1 - S3)^2*Sin[i]^2)/(8*a^5*S1^3) - (39375*e^4*R^4*(S1 - S3)^2*Sin[i]^2)/(64*a^5*S1^3) - (53935*e^3*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]*Sin[i]^2)/(64*a^5*S1) - (53935*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i]^2)/(64*a^5*S1) + (53935*e^3*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]*Sin[i]^2)/(64*a^5*S1^2) + (161805*e^3*R^4*(S1 - S3)^2*Cos[i]*Sin[i]^2)/(128*a^5*S1^3) - (161805*e^3*R^4*(S1 - S3)^3*Cos[i]*Sin[i]^2)/(128*a^5*S1^4) - (315*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(1 - (S1 - S3)^2/S1^2)*Cos[i]^2*Sin[i]^2)/(8*a^5) - (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) + (53935*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(1 - (S1 - S3)^2/S1^2)*Cos[i]^2*Sin[i]^2)/(64*a^5) - (33075*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(1 - (S1 - S3)^2/S1^2)*Cos[i]^2*Sin[i]^2)/(64*a^5) - (315*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^2*Sin[i]^2)/(4*a^5*S1^2) - (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^2*Sin[i]^2)/(4*a^5*S1^2) - (53935*e^3*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^2*Sin[i]^2)/(16*a^5*S1^2) - (33075*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^2*Sin[i]^2)/(32*a^5*S1^2) - (315*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^2*Sin[i]^2)/(8*a^5*S1^2) - (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) - (53935*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^2*Sin[i]^2)/(32*a^5*S1^2) - (33075*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^2*Sin[i]^2)/(64*a^5*S1^2) + (315*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Cos[i]^2*Sin[i]^2)/(4*a^5*S1^3) + (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) + (53935*e^3*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Cos[i]^2*Sin[i]^2)/(16*a^5*S1^3) + (33075*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Cos[i]^2*Sin[i]^2)/(32*a^5*S1^3) - (315*R^4*(S1 - S3)^3*Cos[i]^2*Sin[i]^2)/(4*a^5*S1^4) - (1575*e^2*R^4*(S1 - S3)^3*Cos[i]^2*Sin[i]^2)/(4*a^5*S1^4) + (53935*e^3*R^4*(S1 - S3)^3*Cos[i]^2*Sin[i]^2)/(32*a^5*S1^4) - (33075*e^4*R^4*(S1 - S3)^3*Cos[i]^2*Sin[i]^2)/(32*a^5*S1^4) + (315*R^4*(S1 - S3)^4*Cos[i]^2*Sin[i]^2)/(4*a^5*S1^5) + (1575*e^2*R^4*(S1 - S3)^4*Cos[i]^2*Sin[i]^2)/(4*a^5*S1^5) - (53935*e^3*R^4*(S1 - S3)^4*Cos[i]^2*Sin[i]^2)/(32*a^5*S1^5) + (33075*e^4*R^4*(S1 - S3)^4*Cos[i]^2*Sin[i]^2)/(32*a^5*S1^5) + (53935*e^3*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Sin[i]^3)/(64*a^5*S1^2) + (53935*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Sin[i]^3)/(256*a^5*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) - (53935*e^3*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Sin[i]^3)/(64*a^5*S1^3) + (315*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]*Sin[i]^3)/(8*a^5*S1) + (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]*Sin[i]^3)/(8*a^5*S1) - (53935*e^3*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]*Sin[i]^3)/(64*a^5*S1) + (33075*e^4*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]*Sin[i]^3)/(64*a^5*S1) + (945*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i]^3)/(16*a^5*S1) + (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) - (161805*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i]^3)/(128*a^5*S1) + (99225*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i]^3)/(128*a^5*S1) - (315*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*(S1 - S3)*Cos[i]*Sin[i]^3)/(8*a^5*S1^2) - (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) + (53935*e^3*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*(S1 - S3)*Cos[i]*Sin[i]^3)/(64*a^5*S1^2) - (33075*e^4*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*(S1 - S3)*Cos[i]*Sin[i]^3)/(64*a^5*S1^2) + (945*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]*Sin[i]^3)/(8*a^5*S1^3) + (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) + (161805*e^3*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]*Sin[i]^3)/(64*a^5*S1^3) + (99225*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]*Sin[i]^3)/(64*a^5*S1^3) + (315*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^3*Cos[i]*Sin[i]^3)/(16*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) + (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]) + (53935*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^3*Cos[i]*Sin[i]^3)/(128*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) + (33075*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^3*Cos[i]*Sin[i]^3)/(128*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) - (945*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^3*Cos[i]*Sin[i]^3)/(8*a^5*S1^4) - (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) - (161805*e^3*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^3*Cos[i]*Sin[i]^3)/(64*a^5*S1^4) - (99225*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^3*Cos[i]*Sin[i]^3)/(64*a^5*S1^4) - (315*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Sin[i]^4)/(8*a^5*S1^2) - (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Sin[i]^4)/(8*a^5*S1^2) + (53935*e^3*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Sin[i]^4)/(64*a^5*S1^2) - (33075*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Sin[i]^4)/(64*a^5*S1^2) - (315*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Sin[i]^4)/(16*a^5*S1^2) - (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) + (53935*e^3*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Sin[i]^4)/(128*a^5*S1^2) - (33075*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Sin[i]^4)/(128*a^5*S1^2) + (315*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Sin[i]^4)/(8*a^5*S1^3) + (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Sin[i]^4)/(8*a^5*S1^3) - (53935*e^3*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Sin[i]^4)/(64*a^5*S1^3) + (33075*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Sin[i]^4)/(64*a^5*S1^3) - (315*R^4*(S1 - S3)^3*Sin[i]^4)/(4*a^5*S1^4) - (1575*e^2*R^4*(S1 - S3)^3*Sin[i]^4)/(4*a^5*S1^4) - (33075*e^4*R^4*(S1 - S3)^3*Sin[i]^4)/(32*a^5*S1^4) + (315*R^4*(S1 - S3)^4*Sin[i]^4)/(4*a^5*S1^5) + (1575*e^2*R^4*(S1 - S3)^4*Sin[i]^4)/(4*a^5*S1^5) + (33075*e^4*R^4*(S1 - S3)^4*Sin[i]^4)/(32*a^5*S1^5))) == 0, Omd - G*m*M*(C22*((-845*e^3*R^2*Cos[i])/(32*a^3*S1) - (845*e^3*R^2*(S1 - S3)*Cos[i]^2)/(32*a^3*S1^2) + (845*e^3*R^2*(S1 - S3)*Sin[i])/(32*a^3*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) - (845*e^3*R^2*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(32*a^3*S1) + (845*e^3*R^2*(S1 - S3)^2*Cos[i]*Sin[i])/(32*a^3*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) + (845*e^3*R^2*(S1 - S3)*Sin[i]^2)/(32*a^3*S1^2)) + C20*((3*R^2*(S1 - S3)*Cos[i]^2)/(2*a^3*S1^2) + (9*e^2*R^2*(S1 - S3)*Cos[i]^2)/(4*a^3*S1^2) + (45*e^4*R^2*(S1 - S3)*Cos[i]^2)/(16*a^3*S1^2) + (3*R^2*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(2*a^3*S1) + (9*e^2*R^2*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(4*a^3*S1) + (45*e^4*R^2*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(16*a^3*S1) - (3*R^2*(S1 - S3)^2*Cos[i]*Sin[i])/(2*a^3*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) - (9*e^2*R^2*(S1 - S3)^2*Cos[i]*Sin[i])/(4*a^3*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) - (45*e^4*R^2*(S1 - S3)^2*Cos[i]*Sin[i])/(16*a^3*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) - (3*R^2*(S1 - S3)*Sin[i]^2)/(2*a^3*S1^2) - (9*e^2*R^2*(S1 - S3)*Sin[i]^2)/(4*a^3*S1^2) - (45*e^4*R^2*(S1 - S3)*Sin[i]^2)/(16*a^3*S1^2)) + C40*((-15*R^4*(S1 - S3)*Cos[i]^2)/(4*a^5*S1^2) - (75*e^2*R^4*(S1 - S3)*Cos[i]^2)/(4*a^5*S1^2) - (1575*e^4*R^4*(S1 - S3)*Cos[i]^2)/(32*a^5*S1^2) + (105*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^4)/(16*a^5*S1^2) + (525*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^4)/(16*a^5*S1^2) + (11025*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^4)/(128*a^5*S1^2) - (15*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(4*a^5*S1) - (75*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(4*a^5*S1) - (1575*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(32*a^5*S1) + (15*R^4*(S1 - S3)^2*Cos[i]*Sin[i])/(4*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) + (75*e^2*R^4*(S1 - S3)^2*Cos[i]*Sin[i])/(4*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) + (1575*e^4*R^4*(S1 - S3)^2*Cos[i]*Sin[i])/(32*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) + (105*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(16*a^5*S1) + (525*e^2*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(16*a^5*S1) + (11025*e^4*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(128*a^5*S1) - (315*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^3*Sin[i])/(16*a^5*S1^3) - (1575*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^3*Sin[i])/(16*a^5*S1^3) - (33075*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^3*Sin[i])/(128*a^5*S1^3) + (15*R^4*(S1 - S3)*Sin[i]^2)/(4*a^5*S1^2) + (75*e^2*R^4*(S1 - S3)*Sin[i]^2)/(4*a^5*S1^2) + (1575*e^4*R^4*(S1 - S3)*Sin[i]^2)/(32*a^5*S1^2) - (315*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^2*Sin[i]^2)/(16*a^5*S1^2) - (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^2*Sin[i]^2)/(16*a^5*S1^2) - (33075*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^2*Sin[i]^2)/(128*a^5*S1^2) + (315*R^4*(S1 - S3)^3*Cos[i]^2*Sin[i]^2)/(16*a^5*S1^4) + (1575*e^2*R^4*(S1 - S3)^3*Cos[i]^2*Sin[i]^2)/(16*a^5*S1^4) + (33075*e^4*R^4*(S1 - S3)^3*Cos[i]^2*Sin[i]^2)/(128*a^5*S1^4) + (315*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]*Sin[i]^3)/(16*a^5*S1^3) + (1575*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]*Sin[i]^3)/(16*a^5*S1^3) + (33075*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]*Sin[i]^3)/(128*a^5*S1^3) - (105*R^4*(S1 - S3)^4*Cos[i]*Sin[i]^3)/(16*a^5*S1^5*Sqrt[1 - (S1 - S3)^2/S1^2]) - (525*e^2*R^4*(S1 - S3)^4*Cos[i]*Sin[i]^3)/(16*a^5*S1^5*Sqrt[1 - (S1 - S3)^2/S1^2]) - (11025*e^4*R^4*(S1 - S3)^4*Cos[i]*Sin[i]^3)/(128*a^5*S1^5*Sqrt[1 - (S1 - S3)^2/S1^2]) - (105*R^4*(S1 - S3)^3*Sin[i]^4)/(16*a^5*S1^4) - (525*e^2*R^4*(S1 - S3)^3*Sin[i]^4)/(16*a^5*S1^4) - (11025*e^4*R^4*(S1 - S3)^3*Sin[i]^4)/(128*a^5*S1^4)) + C42*((7705*e^3*R^4*Cos[i])/(64*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) + (7705*e^3*R^4*(S1 - S3)*Cos[i]^2)/(64*a^5*S1^2) + (33075*e^4*R^4*(S1 - S3)*Cos[i]^2)/(64*a^5*S1^2) - (53935*e^3*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]^3)/(128*a^5*S1) + (53935*e^3*R^4*(S1 - S3)^2*Cos[i]^3)/(64*a^5*S1^3) - (315*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^4)/(8*a^5*S1^2) - (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^4)/(8*a^5*S1^2) - (53935*e^3*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^4)/(64*a^5*S1^2) - (33075*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^4)/(64*a^5*S1^2) - (315*R^4*(S1 - S3)^3*Cos[i]^4)/(8*a^5*S1^4) - (1575*e^2*R^4*(S1 - S3)^3*Cos[i]^4)/(8*a^5*S1^4) + (53935*e^3*R^4*(S1 - S3)^3*Cos[i]^4)/(64*a^5*S1^4) - (33075*e^4*R^4*(S1 - S3)^3*Cos[i]^4)/(64*a^5*S1^4) - (7705*e^3*R^4*(S1 - S3)*Sin[i])/(64*a^5*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) + (315*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(8*a^5*S1) + (1575*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(8*a^5*S1) + (7705*e^3*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(64*a^5*S1) + (33075*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(64*a^5*S1) - (315*R^4*(S1 - S3)^2*Cos[i]*Sin[i])/(8*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) - (1575*e^2*R^4*(S1 - S3)^2*Cos[i]*Sin[i])/(8*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) - (7705*e^3*R^4*(S1 - S3)^2*Cos[i]*Sin[i])/(64*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) - (33075*e^4*R^4*(S1 - S3)^2*Cos[i]*Sin[i])/(64*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) + (377545*e^3*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^2*Sin[i])/(128*a^5*S1^2) - (53935*e^3*R^4*(S1 - S3)^3*Cos[i]^2*Sin[i])/(64*a^5*S1^4*Sqrt[1 - (S1 - S3)^2/S1^2]) - (315*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(8*a^5*S1) - (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(8*a^5*S1) - (53935*e^3*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(64*a^5*S1) - (33075*e^4*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(64*a^5*S1) + (161805*e^3*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^3*Sin[i])/(32*a^5*S1^3) + (315*R^4*(S1 - S3)^4*Cos[i]^3*Sin[i])/(8*a^5*S1^5*Sqrt[1 - (S1 - S3)^2/S1^2]) + (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]) - (53935*e^3*R^4*(S1 - S3)^4*Cos[i]^3*Sin[i])/(64*a^5*S1^5*Sqrt[1 - (S1 - S3)^2/S1^2]) + (33075*e^4*R^4*(S1 - S3)^4*Cos[i]^3*Sin[i])/(64*a^5*S1^5*Sqrt[1 - (S1 - S3)^2/S1^2]) - (315*R^4*(S1 - S3)*Sin[i]^2)/(8*a^5*S1^2) - (1575*e^2*R^4*(S1 - S3)*Sin[i]^2)/(8*a^5*S1^2) - (7705*e^3*R^4*(S1 - S3)*Sin[i]^2)/(64*a^5*S1^2) - (33075*e^4*R^4*(S1 - S3)*Sin[i]^2)/(64*a^5*S1^2) + (53935*e^3*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]*Sin[i]^2)/(64*a^5*S1) - (377545*e^3*R^4*(S1 - S3)^2*Cos[i]*Sin[i]^2)/(128*a^5*S1^3) + (161805*e^3*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^2*Sin[i]^2)/(32*a^5*S1^2) - (161805*e^3*R^4*(S1 - S3)^3*Cos[i]^2*Sin[i]^2)/(32*a^5*S1^4) - (53935*e^3*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Sin[i]^3)/(64*a^5*S1^2) + (53935*e^3*R^4*(S1 - S3)^3*Sin[i]^3)/(128*a^5*S1^4*Sqrt[1 - (S1 - S3)^2/S1^2]) - (315*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]*Sin[i]^3)/(8*a^5*S1) - (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]*Sin[i]^3)/(8*a^5*S1) + (53935*e^3*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]*Sin[i]^3)/(64*a^5*S1) - (33075*e^4*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]*Sin[i]^3)/(64*a^5*S1) - (161805*e^3*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]*Sin[i]^3)/(32*a^5*S1^3) + (315*R^4*(S1 - S3)^4*Cos[i]*Sin[i]^3)/(8*a^5*S1^5*Sqrt[1 - (S1 - S3)^2/S1^2]) + (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]) + (53935*e^3*R^4*(S1 - S3)^4*Cos[i]*Sin[i]^3)/(64*a^5*S1^5*Sqrt[1 - (S1 - S3)^2/S1^2]) + (33075*e^4*R^4*(S1 - S3)^4*Cos[i]*Sin[i]^3)/(64*a^5*S1^5*Sqrt[1 - (S1 - S3)^2/S1^2]) + (315*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Sin[i]^4)/(8*a^5*S1^2) + (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Sin[i]^4)/(8*a^5*S1^2) - (53935*e^3*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Sin[i]^4)/(64*a^5*S1^2) + (33075*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Sin[i]^4)/(64*a^5*S1^2) + (315*R^4*(S1 - S3)^3*Sin[i]^4)/(8*a^5*S1^4) + (1575*e^2*R^4*(S1 - S3)^3*Sin[i]^4)/(8*a^5*S1^4) + (53935*e^3*R^4*(S1 - S3)^3*Sin[i]^4)/(64*a^5*S1^4) + (33075*e^4*R^4*(S1 - S3)^3*Sin[i]^4)/(64*a^5*S1^4)) + C44*((-53935*e^3*R^4*Cos[i])/(16*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) - (53935*e^3*R^4*(S1 - S3)*Cos[i]^2)/(8*a^5*S1^2) + (33075*e^4*R^4*(S1 - S3)*Cos[i]^2)/(16*a^5*S1^2) + (53935*e^3*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]^3)/(16*a^5*S1) + (53935*e^3*R^4*(S1 - S3)^2*Cos[i]^3)/(16*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) + (53935*e^3*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^4)/(8*a^5*S1^2) - (33075*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^4)/(16*a^5*S1^2) - (315*R^4*(S1 - S3)^3*Cos[i]^4)/(2*a^5*S1^4) - (1575*e^2*R^4*(S1 - S3)^3*Cos[i]^4)/(2*a^5*S1^4) + (53935*e^3*R^4*(S1 - S3)^3*Cos[i]^4)/(8*a^5*S1^4) - (33075*e^4*R^4*(S1 - S3)^3*Cos[i]^4)/(16*a^5*S1^4) + (53935*e^3*R^4*(S1 - S3)*Sin[i])/(16*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) - (53935*e^3*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(8*a^5*S1) + (33075*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]) + (53935*e^3*R^4*(S1 - S3)^2*Cos[i]*Sin[i])/(8*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) - (33075*e^4*R^4*(S1 - S3)^2*Cos[i]*Sin[i])/(16*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) - (53935*e^3*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^2*Sin[i])/(16*a^5*S1^2) - (53935*e^3*R^4*(S1 - S3)^3*Cos[i]^2*Sin[i])/(16*a^5*S1^4*Sqrt[1 - (S1 - S3)^2/S1^2]) - (315*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^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) + (53935*e^3*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(8*a^5*S1) - (33075*e^4*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(16*a^5*S1) + (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]) - (53935*e^3*R^4*(S1 - S3)^4*Cos[i]^3*Sin[i])/(8*a^5*S1^5*Sqrt[1 - (S1 - S3)^2/S1^2]) + (33075*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) + (53935*e^3*R^4*(S1 - S3)*Sin[i]^2)/(8*a^5*S1^2) - (33075*e^4*R^4*(S1 - S3)*Sin[i]^2)/(16*a^5*S1^2) + (53935*e^3*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]*Sin[i]^2)/(16*a^5*S1) + (53935*e^3*R^4*(S1 - S3)^2*Cos[i]*Sin[i]^2)/(16*a^5*S1^3) - (53935*e^3*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Sin[i]^3)/(16*a^5*S1^2) - (53935*e^3*R^4*(S1 - S3)^3*Sin[i]^3)/(16*a^5*S1^4*Sqrt[1 - (S1 - S3)^2/S1^2]) - (315*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]*Sin[i]^3)/(2*a^5*S1) - (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]*Sin[i]^3)/(2*a^5*S1) + (53935*e^3*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]*Sin[i]^3)/(8*a^5*S1) - (33075*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]) - (53935*e^3*R^4*(S1 - S3)^4*Cos[i]*Sin[i]^3)/(8*a^5*S1^5*Sqrt[1 - (S1 - S3)^2/S1^2]) + (33075*e^4*R^4*(S1 - S3)^4*Cos[i]*Sin[i]^3)/(16*a^5*S1^5*Sqrt[1 - (S1 - S3)^2/S1^2]) + (315*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Sin[i]^4)/(2*a^5*S1^2) + (1575*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Sin[i]^4)/(2*a^5*S1^2) - (53935*e^3*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Sin[i]^4)/(8*a^5*S1^2) + (33075*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Sin[i]^4)/(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) - (53935*e^3*R^4*(S1 - S3)^3*Sin[i]^4)/(8*a^5*S1^4) + (33075*e^4*R^4*(S1 - S3)^3*Sin[i]^4)/(16*a^5*S1^4))) == 0}