sys11 = {-n - omd - Omd + S1/C - G*m*M*(C11*((R*Cos[i])/(2*a^2*S1) - (e^2*R*Cos[i])/(4*a^2*S1) - (e^4*R*Cos[i])/(128*a^2*S1) - (R*(S1 - S3)*Cos[i])/(2*a^2*S1^2) + (e^2*R*(S1 - S3)*Cos[i])/(4*a^2*S1^2) + (e^4*R*(S1 - S3)*Cos[i])/(128*a^2*S1^2) + (R*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i])/(4*a^2*Sqrt[1 - (S1 - S3)^2/S1^2]) - (e^2*R*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i])/(8*a^2*Sqrt[1 - (S1 - S3)^2/S1^2]) - (e^4*R*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i])/(256*a^2*Sqrt[1 - (S1 - S3)^2/S1^2])) + C22*((3*R^2*Cos[i])/(2*a^3*S1) - (15*e^2*R^2*Cos[i])/(4*a^3*S1) + (39*e^4*R^2*Cos[i])/(32*a^3*S1) - (3*R^2*(S1 - S3)*Cos[i])/(2*a^3*S1^2) + (15*e^2*R^2*(S1 - S3)*Cos[i])/(4*a^3*S1^2) - (39*e^4*R^2*(S1 - S3)*Cos[i])/(32*a^3*S1^2) + (3*R^2*(S1 - S3)*Cos[i]^2)/(2*a^3*S1^2) - (15*e^2*R^2*(S1 - S3)*Cos[i]^2)/(4*a^3*S1^2) + (39*e^4*R^2*(S1 - S3)*Cos[i]^2)/(32*a^3*S1^2) - (3*R^2*(S1 - S3)^2*Cos[i]^2)/(2*a^3*S1^3) + (15*e^2*R^2*(S1 - S3)^2*Cos[i]^2)/(4*a^3*S1^3) - (39*e^4*R^2*(S1 - S3)^2*Cos[i]^2)/(32*a^3*S1^3) + (3*R^2*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i])/(4*a^3*Sqrt[1 - (S1 - S3)^2/S1^2]) - (15*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]) + (39*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]) + (3*R^2*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(2*a^3*S1) - (15*e^2*R^2*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(4*a^3*S1) + (39*e^4*R^2*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(32*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]) - (15*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]) + (39*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]) - (3*R^2*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i])/(2*a^3*S1^2) + (15*e^2*R^2*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i])/(4*a^3*S1^2) - (39*e^4*R^2*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i])/(32*a^3*S1^2) + (3*R^2*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i]^2)/(4*a^3) - (15*e^2*R^2*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i]^2)/(8*a^3) + (39*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)) + C33*((45*R^3*Cos[i])/(8*a^4*S1) - (135*e^2*R^3*Cos[i])/(4*a^4*S1) + (19035*e^4*R^3*Cos[i])/(512*a^4*S1) - (45*R^3*(S1 - S3)*Cos[i])/(8*a^4*S1^2) + (135*e^2*R^3*(S1 - S3)*Cos[i])/(4*a^4*S1^2) - (19035*e^4*R^3*(S1 - S3)*Cos[i])/(512*a^4*S1^2) + (45*R^3*(S1 - S3)*Cos[i]^2)/(4*a^4*S1^2) - (135*e^2*R^3*(S1 - S3)*Cos[i]^2)/(2*a^4*S1^2) + (19035*e^4*R^3*(S1 - S3)*Cos[i]^2)/(256*a^4*S1^2) - (45*R^3*(S1 - S3)^2*Cos[i]^2)/(4*a^4*S1^3) + (135*e^2*R^3*(S1 - S3)^2*Cos[i]^2)/(2*a^4*S1^3) - (19035*e^4*R^3*(S1 - S3)^2*Cos[i]^2)/(256*a^4*S1^3) + (45*R^3*(S1 - S3)^2*Cos[i]^3)/(8*a^4*S1^3) - (135*e^2*R^3*(S1 - S3)^2*Cos[i]^3)/(4*a^4*S1^3) + (19035*e^4*R^3*(S1 - S3)^2*Cos[i]^3)/(512*a^4*S1^3) - (45*R^3*(S1 - S3)^3*Cos[i]^3)/(8*a^4*S1^4) + (135*e^2*R^3*(S1 - S3)^3*Cos[i]^3)/(4*a^4*S1^4) - (19035*e^4*R^3*(S1 - S3)^3*Cos[i]^3)/(512*a^4*S1^4) + (45*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i])/(16*a^4*Sqrt[1 - (S1 - S3)^2/S1^2]) - (135*e^2*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]) + (19035*e^4*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i])/(1024*a^4*Sqrt[1 - (S1 - S3)^2/S1^2]) + (45*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(4*a^4*S1) - (135*e^2*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(2*a^4*S1) + (19035*e^4*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(256*a^4*S1) + (45*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i])/(8*a^4*S1*Sqrt[1 - (S1 - S3)^2/S1^2]) - (135*e^2*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i])/(4*a^4*S1*Sqrt[1 - (S1 - S3)^2/S1^2]) + (19035*e^4*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i])/(512*a^4*S1*Sqrt[1 - (S1 - S3)^2/S1^2]) - (45*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i])/(4*a^4*S1^2) + (135*e^2*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i])/(2*a^4*S1^2) - (19035*e^4*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i])/(256*a^4*S1^2) + (45*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^2*Sin[i])/(4*a^4*S1^2) - (135*e^2*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^2*Sin[i])/(2*a^4*S1^2) + (19035*e^4*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^2*Sin[i])/(256*a^4*S1^2) + (45*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]) - (135*e^2*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^2*Sin[i])/(8*a^4*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) + (19035*e^4*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^2*Sin[i])/(1024*a^4*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) - (45*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^2*Sin[i])/(4*a^4*S1^3) + (135*e^2*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^2*Sin[i])/(2*a^4*S1^3) - (19035*e^4*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^2*Sin[i])/(256*a^4*S1^3) + (45*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i]^2)/(8*a^4) - (135*e^2*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i]^2)/(4*a^4) + (19035*e^4*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i]^2)/(512*a^4) + (45*R^3*(1 - (S1 - S3)^2/S1^2)*Cos[i]*Sin[i]^2)/(8*a^4*S1) - (135*e^2*R^3*(1 - (S1 - S3)^2/S1^2)*Cos[i]*Sin[i]^2)/(4*a^4*S1) + (19035*e^4*R^3*(1 - (S1 - S3)^2/S1^2)*Cos[i]*Sin[i]^2)/(512*a^4*S1) + (45*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i]^2)/(8*a^4*S1) - (135*e^2*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i]^2)/(4*a^4*S1) + (19035*e^4*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i]^2)/(512*a^4*S1) - (45*R^3*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]*Sin[i]^2)/(8*a^4*S1^2) + (135*e^2*R^3*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]*Sin[i]^2)/(4*a^4*S1^2) - (19035*e^4*R^3*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]*Sin[i]^2)/(512*a^4*S1^2) + (45*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sqrt[1 - (S1 - S3)^2/S1^2]*Sin[i]^3)/(16*a^4) - (135*e^2*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sqrt[1 - (S1 - S3)^2/S1^2]*Sin[i]^3)/(8*a^4) + (19035*e^4*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sqrt[1 - (S1 - S3)^2/S1^2]*Sin[i]^3)/(1024*a^4)) + C31*((-3*R^3*Cos[i])/(4*a^4*S1) - (3*e^2*R^3*Cos[i])/(2*a^4*S1) - (717*e^4*R^3*Cos[i])/(256*a^4*S1) + (3*R^3*(S1 - S3)*Cos[i])/(4*a^4*S1^2) + (3*e^2*R^3*(S1 - S3)*Cos[i])/(2*a^4*S1^2) + (717*e^4*R^3*(S1 - S3)*Cos[i])/(256*a^4*S1^2) + (15*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Cos[i]^2)/(16*a^4) + (15*e^2*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Cos[i]^2)/(8*a^4) + (3585*e^4*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Cos[i]^2)/(1024*a^4) + (45*R^3*(1 - (S1 - S3)^2/S1^2)*Cos[i]^3)/(16*a^4*S1) + (45*e^2*R^3*(1 - (S1 - S3)^2/S1^2)*Cos[i]^3)/(8*a^4*S1) + (10755*e^4*R^3*(1 - (S1 - S3)^2/S1^2)*Cos[i]^3)/(1024*a^4*S1) + (45*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]^3)/(16*a^4*S1) + (45*e^2*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]^3)/(8*a^4*S1) + (10755*e^4*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]^3)/(1024*a^4*S1) - (45*R^3*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^3)/(16*a^4*S1^2) - (45*e^2*R^3*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^3)/(8*a^4*S1^2) - (10755*e^4*R^3*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^3)/(1024*a^4*S1^2) - (3*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]) - (3*e^2*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i])/(4*a^4*Sqrt[1 - (S1 - S3)^2/S1^2]) - (717*e^4*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i])/(512*a^4*Sqrt[1 - (S1 - S3)^2/S1^2]) - (15*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(8*a^4*S1) - (15*e^2*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(4*a^4*S1) - (3585*e^4*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(512*a^4*S1) - (15*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]) - (15*e^2*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i])/(8*a^4*S1*Sqrt[1 - (S1 - S3)^2/S1^2]) - (3585*e^4*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i])/(1024*a^4*S1*Sqrt[1 - (S1 - S3)^2/S1^2]) + (15*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i])/(8*a^4*S1^2) + (15*e^2*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i])/(4*a^4*S1^2) + (3585*e^4*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i])/(512*a^4*S1^2) + (135*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) + (135*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])/(16*a^4) + (32265*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])/(2048*a^4) - (45*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^2*Sin[i])/(4*a^4*S1^2) - (45*e^2*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^2*Sin[i])/(2*a^4*S1^2) - (10755*e^4*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^2*Sin[i])/(256*a^4*S1^2) - (45*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]) - (45*e^2*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^2*Sin[i])/(8*a^4*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) - (10755*e^4*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^2*Sin[i])/(1024*a^4*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) + (45*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^2*Sin[i])/(4*a^4*S1^3) + (45*e^2*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^2*Sin[i])/(2*a^4*S1^3) + (10755*e^4*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^2*Sin[i])/(256*a^4*S1^3) + (15*R^3*(S1 - S3)*Sin[i]^2)/(8*a^4*S1^2) + (15*e^2*R^3*(S1 - S3)*Sin[i]^2)/(4*a^4*S1^2) + (3585*e^4*R^3*(S1 - S3)*Sin[i]^2)/(512*a^4*S1^2) - (15*R^3*(S1 - S3)^2*Sin[i]^2)/(8*a^4*S1^3) - (15*e^2*R^3*(S1 - S3)^2*Sin[i]^2)/(4*a^4*S1^3) - (3585*e^4*R^3*(S1 - S3)^2*Sin[i]^2)/(512*a^4*S1^3) - (45*R^3*(1 - (S1 - S3)^2/S1^2)*Cos[i]*Sin[i]^2)/(8*a^4*S1) - (45*e^2*R^3*(1 - (S1 - S3)^2/S1^2)*Cos[i]*Sin[i]^2)/(4*a^4*S1) - (10755*e^4*R^3*(1 - (S1 - S3)^2/S1^2)*Cos[i]*Sin[i]^2)/(512*a^4*S1) - (45*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i]^2)/(8*a^4*S1) - (45*e^2*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i]^2)/(4*a^4*S1) - (10755*e^4*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i]^2)/(512*a^4*S1) + (45*R^3*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]*Sin[i]^2)/(8*a^4*S1^2) + (45*e^2*R^3*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]*Sin[i]^2)/(4*a^4*S1^2) + (10755*e^4*R^3*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]*Sin[i]^2)/(512*a^4*S1^2) + (135*R^3*(S1 - S3)^2*Cos[i]*Sin[i]^2)/(16*a^4*S1^3) + (135*e^2*R^3*(S1 - S3)^2*Cos[i]*Sin[i]^2)/(8*a^4*S1^3) + (32265*e^4*R^3*(S1 - S3)^2*Cos[i]*Sin[i]^2)/(1024*a^4*S1^3) - (135*R^3*(S1 - S3)^3*Cos[i]*Sin[i]^2)/(16*a^4*S1^4) - (135*e^2*R^3*(S1 - S3)^3*Cos[i]*Sin[i]^2)/(8*a^4*S1^4) - (32265*e^4*R^3*(S1 - S3)^3*Cos[i]*Sin[i]^2)/(1024*a^4*S1^4) + (45*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Sin[i]^3)/(8*a^4*S1^2) + (45*e^2*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Sin[i]^3)/(4*a^4*S1^2) + (10755*e^4*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Sin[i]^3)/(512*a^4*S1^2) + (45*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]) + (45*e^2*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Sin[i]^3)/(16*a^4*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) + (10755*e^4*R^3*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Sin[i]^3)/(2048*a^4*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) - (45*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Sin[i]^3)/(8*a^4*S1^3) - (45*e^2*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Sin[i]^3)/(4*a^4*S1^3) - (10755*e^4*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Sin[i]^3)/(512*a^4*S1^3)) + C44*((525*R^4*Cos[i])/(4*a^5*S1) - (735*e^2*R^4*Cos[i])/(4*a^5*S1) + (34545*e^4*R^4*Cos[i])/(32*a^5*S1) - (525*R^4*(S1 - S3)*Cos[i])/(4*a^5*S1^2) + (735*e^2*R^4*(S1 - S3)*Cos[i])/(4*a^5*S1^2) - (34545*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) + (105*R^4*(S1 - S3)*Cos[i]^2)/(4*a^5*S1^2) - (7875*e^2*R^4*(S1 - S3)*Cos[i]^2)/(4*a^5*S1^2) - (20265*e^4*R^4*(S1 - S3)*Cos[i]^2)/(32*a^5*S1^2) - (105*R^4*(S1 - S3)^2*Cos[i]^2)/(4*a^5*S1^3) + (7875*e^2*R^4*(S1 - S3)^2*Cos[i]^2)/(4*a^5*S1^3) + (20265*e^4*R^4*(S1 - S3)^2*Cos[i]^2)/(32*a^5*S1^3) - (105*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]^3)/(a^5*S1) - (105*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]^3)/(a^5*S1) - (6825*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]^3)/(16*a^5*S1) - (105*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]^3)/(a^5*S1) - (105*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]^3)/(a^5*S1) - (6825*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]^3)/(16*a^5*S1) + (105*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^3)/(a^5*S1^2) + (105*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^3)/(a^5*S1^2) + (6825*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^3)/(16*a^5*S1^2) - (945*R^4*(S1 - S3)^2*Cos[i]^3)/(4*a^5*S1^3) - (4725*e^2*R^4*(S1 - S3)^2*Cos[i]^3)/(4*a^5*S1^3) + (21735*e^4*R^4*(S1 - S3)^2*Cos[i]^3)/(32*a^5*S1^3) + (945*R^4*(S1 - S3)^3*Cos[i]^3)/(4*a^5*S1^4) + (4725*e^2*R^4*(S1 - S3)^3*Cos[i]^3)/(4*a^5*S1^4) - (21735*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) + (105*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^4)/(a^5*S1^2) + (1365*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^4)/(a^5*S1^2) + (13125*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^4)/(4*a^5*S1^2) + (105*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^4)/(2*a^5*S1^2) + (1365*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) + (13125*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^4)/(8*a^5*S1^2) - (105*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Cos[i]^4)/(a^5*S1^3) - (1365*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Cos[i]^4)/(a^5*S1^3) - (13125*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Cos[i]^4)/(4*a^5*S1^3) + (315*R^4*(S1 - S3)^3*Cos[i]^4)/(4*a^5*S1^4) + (6615*e^2*R^4*(S1 - S3)^3*Cos[i]^4)/(4*a^5*S1^4) + (164745*e^4*R^4*(S1 - S3)^3*Cos[i]^4)/(32*a^5*S1^4) - (315*R^4*(S1 - S3)^4*Cos[i]^4)/(4*a^5*S1^5) - (6615*e^2*R^4*(S1 - S3)^4*Cos[i]^4)/(4*a^5*S1^5) - (164745*e^4*R^4*(S1 - S3)^4*Cos[i]^4)/(32*a^5*S1^5) + (525*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]) - (735*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]) + (34545*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]) + (525*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(4*a^5*S1) - (5775*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(4*a^5*S1) + (23835*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(32*a^5*S1) + (525*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]) - (5775*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]) + (23835*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]) - (525*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i])/(4*a^5*S1^2) + (5775*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i])/(4*a^5*S1^2) - (23835*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i])/(32*a^5*S1^2) - (315*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])/(2*a^5) - (315*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])/(2*a^5) - (20475*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])/(32*a^5) - (105*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^2*Sin[i])/(2*a^5*S1^2) - (3885*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^2*Sin[i])/(2*a^5*S1^2) + (49035*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^2*Sin[i])/(16*a^5*S1^2) - (105*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]) - (3885*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]) + (49035*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]) + (105*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^2*Sin[i])/(2*a^5*S1^3) + (3885*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^2*Sin[i])/(2*a^5*S1^3) - (49035*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^2*Sin[i])/(16*a^5*S1^3) - (105*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(2*a^5*S1) + (1155*e^2*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(2*a^5*S1) + (19425*e^4*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(16*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])/(4*a^5*S1) + (3465*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) + (58275*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) + (105*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*(S1 - S3)*Cos[i]^3*Sin[i])/(2*a^5*S1^2) - (1155*e^2*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*(S1 - S3)*Cos[i]^3*Sin[i])/(2*a^5*S1^2) - (19425*e^4*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*(S1 - S3)*Cos[i]^3*Sin[i])/(16*a^5*S1^2) - (315*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^3*Sin[i])/(4*a^5*S1^3) + (3465*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^3*Sin[i])/(4*a^5*S1^3) + (179235*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^3*Sin[i])/(32*a^5*S1^3) - (105*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^3*Cos[i]^3*Sin[i])/(8*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) + (1155*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^3*Cos[i]^3*Sin[i])/(8*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) + (59745*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]) + (315*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^3*Cos[i]^3*Sin[i])/(4*a^5*S1^4) - (3465*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^3*Cos[i]^3*Sin[i])/(4*a^5*S1^4) - (179235*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^3*Cos[i]^3*Sin[i])/(32*a^5*S1^4) + (105*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i]^2)/(8*a^5) - (7875*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i]^2)/(8*a^5) - (20265*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) - (105*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]*Sin[i]^2)/(4*a^5*S1) - (3885*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]*Sin[i]^2)/(4*a^5*S1) + (49035*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]*Sin[i]^2)/(32*a^5*S1) - (105*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i]^2)/(4*a^5*S1) - (3885*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) + (49035*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) + (105*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]*Sin[i]^2)/(4*a^5*S1^2) + (3885*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]*Sin[i]^2)/(4*a^5*S1^2) - (49035*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]*Sin[i]^2)/(32*a^5*S1^2) - (315*R^4*(S1 - S3)^2*Cos[i]*Sin[i]^2)/(a^5*S1^3) - (315*e^2*R^4*(S1 - S3)^2*Cos[i]*Sin[i]^2)/(a^5*S1^3) - (20475*e^4*R^4*(S1 - S3)^2*Cos[i]*Sin[i]^2)/(16*a^5*S1^3) + (315*R^4*(S1 - S3)^3*Cos[i]*Sin[i]^2)/(a^5*S1^4) + (315*e^2*R^4*(S1 - S3)^3*Cos[i]*Sin[i]^2)/(a^5*S1^4) + (20475*e^4*R^4*(S1 - S3)^3*Cos[i]*Sin[i]^2)/(16*a^5*S1^4) + (105*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 + (1365*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 + (13125*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)/(4*a^5) + (1155*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^2*Sin[i]^2)/(4*a^5*S1^2) + (7455*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^2*Sin[i]^2)/(4*a^5*S1^2) + (272685*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^2*Sin[i]^2)/(32*a^5*S1^2) + (1155*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) + (7455*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) + (272685*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) - (1155*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Cos[i]^2*Sin[i]^2)/(4*a^5*S1^3) - (7455*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Cos[i]^2*Sin[i]^2)/(4*a^5*S1^3) - (272685*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Cos[i]^2*Sin[i]^2)/(32*a^5*S1^3) + (210*R^4*(S1 - S3)^3*Cos[i]^2*Sin[i]^2)/(a^5*S1^4) + (2730*e^2*R^4*(S1 - S3)^3*Cos[i]^2*Sin[i]^2)/(a^5*S1^4) + (13125*e^4*R^4*(S1 - S3)^3*Cos[i]^2*Sin[i]^2)/(2*a^5*S1^4) - (210*R^4*(S1 - S3)^4*Cos[i]^2*Sin[i]^2)/(a^5*S1^5) - (2730*e^2*R^4*(S1 - S3)^4*Cos[i]^2*Sin[i]^2)/(a^5*S1^5) - (13125*e^4*R^4*(S1 - S3)^4*Cos[i]^2*Sin[i]^2)/(2*a^5*S1^5) - (945*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) - (4725*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) + (21735*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) - (210*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Sin[i]^3)/(a^5*S1^2) - (210*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Sin[i]^3)/(a^5*S1^2) - (6825*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Sin[i]^3)/(8*a^5*S1^2) - (105*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Sin[i]^3)/(2*a^5*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) - (105*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Sin[i]^3)/(2*a^5*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) - (6825*e^4*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]) + (210*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Sin[i]^3)/(a^5*S1^3) + (210*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Sin[i]^3)/(a^5*S1^3) + (6825*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Sin[i]^3)/(8*a^5*S1^3) - (105*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]*Sin[i]^3)/(4*a^5*S1) + (1155*e^2*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]*Sin[i]^3)/(4*a^5*S1) + (59745*e^4*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]*Sin[i]^3)/(32*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]*Sin[i]^3)/(8*a^5*S1) + (3465*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)/(8*a^5*S1) + (179235*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) + (105*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*(S1 - S3)*Cos[i]*Sin[i]^3)/(4*a^5*S1^2) - (1155*e^2*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*(S1 - S3)*Cos[i]*Sin[i]^3)/(4*a^5*S1^2) - (59745*e^4*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*(S1 - S3)*Cos[i]*Sin[i]^3)/(32*a^5*S1^2) - (315*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]*Sin[i]^3)/(2*a^5*S1^3) + (3465*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]*Sin[i]^3)/(2*a^5*S1^3) + (58275*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]*Sin[i]^3)/(16*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)/(4*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) + (1155*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]) + (19425*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]) + (315*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^3*Cos[i]*Sin[i]^3)/(2*a^5*S1^4) - (3465*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^3*Cos[i]*Sin[i]^3)/(2*a^5*S1^4) - (58275*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^3*Cos[i]*Sin[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)*Sin[i]^4)/(8*a^5) + (6615*e^2*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) + (164745*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) + (105*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Sin[i]^4)/(a^5*S1^2) + (1365*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Sin[i]^4)/(a^5*S1^2) + (13125*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Sin[i]^4)/(4*a^5*S1^2) + (105*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Sin[i]^4)/(2*a^5*S1^2) + (1365*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) + (13125*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Sin[i]^4)/(8*a^5*S1^2) - (105*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Sin[i]^4)/(a^5*S1^3) - (1365*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Sin[i]^4)/(a^5*S1^3) - (13125*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Sin[i]^4)/(4*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*((-15*R^4*Cos[i])/(4*a^5*S1) - (15*e^2*R^4*Cos[i])/(4*a^5*S1) - (975*e^4*R^4*Cos[i])/(64*a^5*S1) + (15*R^4*(S1 - S3)*Cos[i])/(4*a^5*S1^2) + (15*e^2*R^4*(S1 - S3)*Cos[i])/(4*a^5*S1^2) + (975*e^4*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)/(4*a^5*S1^2) + (135*e^2*R^4*(S1 - S3)*Cos[i]^2)/(4*a^5*S1^2) + (5325*e^4*R^4*(S1 - S3)*Cos[i]^2)/(64*a^5*S1^2) - (15*R^4*(S1 - S3)^2*Cos[i]^2)/(4*a^5*S1^3) - (135*e^2*R^4*(S1 - S3)^2*Cos[i]^2)/(4*a^5*S1^3) - (5325*e^4*R^4*(S1 - S3)^2*Cos[i]^2)/(64*a^5*S1^3) + (105*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]^3)/(8*a^5*S1) + (105*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]^3)/(8*a^5*S1) + (6825*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]^3)/(128*a^5*S1) + (105*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]^3)/(8*a^5*S1) + (105*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]^3)/(8*a^5*S1) + (6825*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]^3)/(128*a^5*S1) - (105*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^3)/(8*a^5*S1^2) - (105*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^3)/(8*a^5*S1^2) - (6825*e^4*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) - (105*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^4)/(8*a^5*S1^2) - (1365*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^4)/(8*a^5*S1^2) - (13125*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^4)/(32*a^5*S1^2) - (105*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^4)/(16*a^5*S1^2) - (1365*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) - (13125*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Cos[i]^4)/(64*a^5*S1^2) + (105*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Cos[i]^4)/(8*a^5*S1^3) + (1365*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Cos[i]^4)/(8*a^5*S1^3) + (13125*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Cos[i]^4)/(32*a^5*S1^3) - (15*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]) - (15*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]) - (975*e^4*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]) - (345*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(8*a^5*S1) - (1605*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(8*a^5*S1) - (17025*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(32*a^5*S1) - (345*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]) - (1605*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]) - (17025*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]) + (345*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i])/(8*a^5*S1^2) + (1605*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i])/(8*a^5*S1^2) + (17025*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]*Sin[i])/(32*a^5*S1^2) + (315*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) + (315*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])/(16*a^5) + (20475*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])/(256*a^5) - (105*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^2*Sin[i])/(2*a^5*S1^2) - (105*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^2*Sin[i])/(2*a^5*S1^2) - (6825*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^2*Sin[i])/(32*a^5*S1^2) - (105*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]) - (105*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]) - (6825*e^4*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]) + (105*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^2*Sin[i])/(2*a^5*S1^3) + (105*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^2*Sin[i])/(2*a^5*S1^3) + (6825*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^2*Sin[i])/(32*a^5*S1^3) + (525*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(8*a^5*S1) + (1785*e^2*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(8*a^5*S1) + (9975*e^4*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(16*a^5*S1) + (1575*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) + (5355*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) + (29925*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) - (525*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*(S1 - S3)*Cos[i]^3*Sin[i])/(8*a^5*S1^2) - (1785*e^2*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*(S1 - S3)*Cos[i]^3*Sin[i])/(8*a^5*S1^2) - (9975*e^4*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*(S1 - S3)*Cos[i]^3*Sin[i])/(16*a^5*S1^2) + (315*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^3*Sin[i])/(8*a^5*S1^3) + (4095*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^3*Sin[i])/(8*a^5*S1^3) + (39375*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^3*Sin[i])/(32*a^5*S1^3) + (105*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]) + (1365*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]) + (13125*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]) - (315*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^3*Cos[i]^3*Sin[i])/(8*a^5*S1^4) - (4095*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^3*Cos[i]^3*Sin[i])/(8*a^5*S1^4) - (39375*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^3*Cos[i]^3*Sin[i])/(32*a^5*S1^4) + (15*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i]^2)/(8*a^5) + (135*e^2*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i]^2)/(8*a^5) + (5325*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*Sin[i]^2)/(128*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) - (105*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]*Sin[i]^2)/(4*a^5*S1) - (105*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]*Sin[i]^2)/(4*a^5*S1) - (6825*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]*Sin[i]^2)/(64*a^5*S1) - (105*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i]^2)/(4*a^5*S1) - (105*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) - (6825*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)*Cos[i]*Sin[i]^2)/(64*a^5*S1) + (105*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]*Sin[i]^2)/(4*a^5*S1^2) + (105*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]*Sin[i]^2)/(4*a^5*S1^2) + (6825*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]*Sin[i]^2)/(64*a^5*S1^2) + (315*R^4*(S1 - S3)^2*Cos[i]*Sin[i]^2)/(8*a^5*S1^3) + (315*e^2*R^4*(S1 - S3)^2*Cos[i]*Sin[i]^2)/(8*a^5*S1^3) + (20475*e^4*R^4*(S1 - S3)^2*Cos[i]*Sin[i]^2)/(128*a^5*S1^3) - (315*R^4*(S1 - S3)^3*Cos[i]*Sin[i]^2)/(8*a^5*S1^4) - (315*e^2*R^4*(S1 - S3)^3*Cos[i]*Sin[i]^2)/(8*a^5*S1^4) - (20475*e^4*R^4*(S1 - S3)^3*Cos[i]*Sin[i]^2)/(128*a^5*S1^4) - (105*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) - (1365*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) - (13125*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)/(32*a^5) - (735*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^2*Sin[i]^2)/(4*a^5*S1^2) - (1995*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^2*Sin[i]^2)/(4*a^5*S1^2) - (46725*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^2*Sin[i]^2)/(32*a^5*S1^2) - (735*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) - (1995*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) - (46725*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) + (735*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Cos[i]^2*Sin[i]^2)/(4*a^5*S1^3) + (1995*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Cos[i]^2*Sin[i]^2)/(4*a^5*S1^3) + (46725*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Cos[i]^2*Sin[i]^2)/(32*a^5*S1^3) - (105*R^4*(S1 - S3)^3*Cos[i]^2*Sin[i]^2)/(4*a^5*S1^4) - (1365*e^2*R^4*(S1 - S3)^3*Cos[i]^2*Sin[i]^2)/(4*a^5*S1^4) - (13125*e^4*R^4*(S1 - S3)^3*Cos[i]^2*Sin[i]^2)/(16*a^5*S1^4) + (105*R^4*(S1 - S3)^4*Cos[i]^2*Sin[i]^2)/(4*a^5*S1^5) + (1365*e^2*R^4*(S1 - S3)^4*Cos[i]^2*Sin[i]^2)/(4*a^5*S1^5) + (13125*e^4*R^4*(S1 - S3)^4*Cos[i]^2*Sin[i]^2)/(16*a^5*S1^5) + (105*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Sin[i]^3)/(4*a^5*S1^2) + (105*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Sin[i]^3)/(4*a^5*S1^2) + (6825*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Sin[i]^3)/(64*a^5*S1^2) + (105*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]) + (105*e^2*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]) + (6825*e^4*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]) - (105*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Sin[i]^3)/(4*a^5*S1^3) - (105*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Sin[i]^3)/(4*a^5*S1^3) - (6825*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Sin[i]^3)/(64*a^5*S1^3) + (105*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]*Sin[i]^3)/(8*a^5*S1) + (1365*e^2*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]*Sin[i]^3)/(8*a^5*S1) + (13125*e^4*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]*Sin[i]^3)/(32*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]*Sin[i]^3)/(16*a^5*S1) + (4095*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) + (39375*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) - (105*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*(S1 - S3)*Cos[i]*Sin[i]^3)/(8*a^5*S1^2) - (1365*e^2*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*(S1 - S3)*Cos[i]*Sin[i]^3)/(8*a^5*S1^2) - (13125*e^4*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*(S1 - S3)*Cos[i]*Sin[i]^3)/(32*a^5*S1^2) + (1575*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]*Sin[i]^3)/(8*a^5*S1^3) + (5355*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]*Sin[i]^3)/(8*a^5*S1^3) + (29925*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]*Sin[i]^3)/(16*a^5*S1^3) + (525*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]) + (1785*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]) + (9975*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]) - (1575*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^3*Cos[i]*Sin[i]^3)/(8*a^5*S1^4) - (5355*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^3*Cos[i]*Sin[i]^3)/(8*a^5*S1^4) - (29925*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^3*Cos[i]*Sin[i]^3)/(16*a^5*S1^4) - (105*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Sin[i]^4)/(8*a^5*S1^2) - (1365*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Sin[i]^4)/(8*a^5*S1^2) - (13125*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Sin[i]^4)/(32*a^5*S1^2) - (105*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Sin[i]^4)/(16*a^5*S1^2) - (1365*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) - (13125*e^4*R^4*((-2*(S1 - S3))/S1^2 + (2*(S1 - S3)^2)/S1^3)*(S1 - S3)^2*Sin[i]^4)/(64*a^5*S1^2) + (105*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Sin[i]^4)/(8*a^5*S1^3) + (1365*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Sin[i]^4)/(8*a^5*S1^3) + (13125*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)^2*Sin[i]^4)/(32*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*(-(R*Cos[i])/(2*a^2*S1) + (e^2*R*Cos[i])/(4*a^2*S1) + (e^4*R*Cos[i])/(128*a^2*S1) + (R*(S1 - S3)*Sin[i])/(2*a^2*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) - (e^2*R*(S1 - S3)*Sin[i])/(4*a^2*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) - (e^4*R*(S1 - S3)*Sin[i])/(128*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*((-3*R^2*Cos[i])/(2*a^3*S1) + (15*e^2*R^2*Cos[i])/(4*a^3*S1) - (39*e^4*R^2*Cos[i])/(32*a^3*S1) - (3*R^2*(S1 - S3)*Cos[i]^2)/(2*a^3*S1^2) + (15*e^2*R^2*(S1 - S3)*Cos[i]^2)/(4*a^3*S1^2) - (39*e^4*R^2*(S1 - S3)*Cos[i]^2)/(32*a^3*S1^2) + (3*R^2*(S1 - S3)*Sin[i])/(2*a^3*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) - (15*e^2*R^2*(S1 - S3)*Sin[i])/(4*a^3*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) + (39*e^4*R^2*(S1 - S3)*Sin[i])/(32*a^3*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) - (3*R^2*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(2*a^3*S1) + (15*e^2*R^2*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(4*a^3*S1) - (39*e^4*R^2*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(32*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]) - (15*e^2*R^2*(S1 - S3)^2*Cos[i]*Sin[i])/(4*a^3*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) + (39*e^4*R^2*(S1 - S3)^2*Cos[i]*Sin[i])/(32*a^3*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) + (3*R^2*(S1 - S3)*Sin[i]^2)/(2*a^3*S1^2) - (15*e^2*R^2*(S1 - S3)*Sin[i]^2)/(4*a^3*S1^2) + (39*e^4*R^2*(S1 - S3)*Sin[i]^2)/(32*a^3*S1^2)) + C33*((-45*R^3*Cos[i])/(8*a^4*S1) + (135*e^2*R^3*Cos[i])/(4*a^4*S1) - (19035*e^4*R^3*Cos[i])/(512*a^4*S1) - (45*R^3*(S1 - S3)*Cos[i]^2)/(4*a^4*S1^2) + (135*e^2*R^3*(S1 - S3)*Cos[i]^2)/(2*a^4*S1^2) - (19035*e^4*R^3*(S1 - S3)*Cos[i]^2)/(256*a^4*S1^2) - (45*R^3*(S1 - S3)^2*Cos[i]^3)/(8*a^4*S1^3) + (135*e^2*R^3*(S1 - S3)^2*Cos[i]^3)/(4*a^4*S1^3) - (19035*e^4*R^3*(S1 - S3)^2*Cos[i]^3)/(512*a^4*S1^3) + (45*R^3*(S1 - S3)*Sin[i])/(8*a^4*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) - (135*e^2*R^3*(S1 - S3)*Sin[i])/(4*a^4*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) + (19035*e^4*R^3*(S1 - S3)*Sin[i])/(512*a^4*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) - (45*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(4*a^4*S1) + (135*e^2*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(2*a^4*S1) - (19035*e^4*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(256*a^4*S1) + (45*R^3*(S1 - S3)^2*Cos[i]*Sin[i])/(4*a^4*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) - (135*e^2*R^3*(S1 - S3)^2*Cos[i]*Sin[i])/(2*a^4*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) + (19035*e^4*R^3*(S1 - S3)^2*Cos[i]*Sin[i])/(256*a^4*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) - (45*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^2*Sin[i])/(4*a^4*S1^2) + (135*e^2*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^2*Sin[i])/(2*a^4*S1^2) - (19035*e^4*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^2*Sin[i])/(256*a^4*S1^2) + (45*R^3*(S1 - S3)^3*Cos[i]^2*Sin[i])/(8*a^4*S1^4*Sqrt[1 - (S1 - S3)^2/S1^2]) - (135*e^2*R^3*(S1 - S3)^3*Cos[i]^2*Sin[i])/(4*a^4*S1^4*Sqrt[1 - (S1 - S3)^2/S1^2]) + (19035*e^4*R^3*(S1 - S3)^3*Cos[i]^2*Sin[i])/(512*a^4*S1^4*Sqrt[1 - (S1 - S3)^2/S1^2]) + (45*R^3*(S1 - S3)*Sin[i]^2)/(4*a^4*S1^2) - (135*e^2*R^3*(S1 - S3)*Sin[i]^2)/(2*a^4*S1^2) + (19035*e^4*R^3*(S1 - S3)*Sin[i]^2)/(256*a^4*S1^2) - (45*R^3*(1 - (S1 - S3)^2/S1^2)*Cos[i]*Sin[i]^2)/(8*a^4*S1) + (135*e^2*R^3*(1 - (S1 - S3)^2/S1^2)*Cos[i]*Sin[i]^2)/(4*a^4*S1) - (19035*e^4*R^3*(1 - (S1 - S3)^2/S1^2)*Cos[i]*Sin[i]^2)/(512*a^4*S1) + (45*R^3*(S1 - S3)^2*Cos[i]*Sin[i]^2)/(4*a^4*S1^3) - (135*e^2*R^3*(S1 - S3)^2*Cos[i]*Sin[i]^2)/(2*a^4*S1^3) + (19035*e^4*R^3*(S1 - S3)^2*Cos[i]*Sin[i]^2)/(256*a^4*S1^3) + (45*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Sin[i]^3)/(8*a^4*S1^2) - (135*e^2*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Sin[i]^3)/(4*a^4*S1^2) + (19035*e^4*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Sin[i]^3)/(512*a^4*S1^2)) + C31*((3*R^3*Cos[i])/(4*a^4*S1) + (3*e^2*R^3*Cos[i])/(2*a^4*S1) + (717*e^4*R^3*Cos[i])/(256*a^4*S1) + (15*R^3*(S1 - S3)*Cos[i]^2)/(8*a^4*S1^2) + (15*e^2*R^3*(S1 - S3)*Cos[i]^2)/(4*a^4*S1^2) + (3585*e^4*R^3*(S1 - S3)*Cos[i]^2)/(512*a^4*S1^2) - (45*R^3*(1 - (S1 - S3)^2/S1^2)*Cos[i]^3)/(16*a^4*S1) - (45*e^2*R^3*(1 - (S1 - S3)^2/S1^2)*Cos[i]^3)/(8*a^4*S1) - (10755*e^4*R^3*(1 - (S1 - S3)^2/S1^2)*Cos[i]^3)/(1024*a^4*S1) + (45*R^3*(S1 - S3)^2*Cos[i]^3)/(8*a^4*S1^3) + (45*e^2*R^3*(S1 - S3)^2*Cos[i]^3)/(4*a^4*S1^3) + (10755*e^4*R^3*(S1 - S3)^2*Cos[i]^3)/(512*a^4*S1^3) - (3*R^3*(S1 - S3)*Sin[i])/(4*a^4*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) - (3*e^2*R^3*(S1 - S3)*Sin[i])/(2*a^4*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) - (717*e^4*R^3*(S1 - S3)*Sin[i])/(256*a^4*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) + (15*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(8*a^4*S1) + (15*e^2*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(4*a^4*S1) + (3585*e^4*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(512*a^4*S1) - (15*R^3*(S1 - S3)^2*Cos[i]*Sin[i])/(8*a^4*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) - (15*e^2*R^3*(S1 - S3)^2*Cos[i]*Sin[i])/(4*a^4*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) - (3585*e^4*R^3*(S1 - S3)^2*Cos[i]*Sin[i])/(512*a^4*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) + (315*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^2*Sin[i])/(16*a^4*S1^2) + (315*e^2*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^2*Sin[i])/(8*a^4*S1^2) + (75285*e^4*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^2*Sin[i])/(1024*a^4*S1^2) - (45*R^3*(S1 - S3)^3*Cos[i]^2*Sin[i])/(8*a^4*S1^4*Sqrt[1 - (S1 - S3)^2/S1^2]) - (45*e^2*R^3*(S1 - S3)^3*Cos[i]^2*Sin[i])/(4*a^4*S1^4*Sqrt[1 - (S1 - S3)^2/S1^2]) - (10755*e^4*R^3*(S1 - S3)^3*Cos[i]^2*Sin[i])/(512*a^4*S1^4*Sqrt[1 - (S1 - S3)^2/S1^2]) - (15*R^3*(S1 - S3)*Sin[i]^2)/(8*a^4*S1^2) - (15*e^2*R^3*(S1 - S3)*Sin[i]^2)/(4*a^4*S1^2) - (3585*e^4*R^3*(S1 - S3)*Sin[i]^2)/(512*a^4*S1^2) + (45*R^3*(1 - (S1 - S3)^2/S1^2)*Cos[i]*Sin[i]^2)/(8*a^4*S1) + (45*e^2*R^3*(1 - (S1 - S3)^2/S1^2)*Cos[i]*Sin[i]^2)/(4*a^4*S1) + (10755*e^4*R^3*(1 - (S1 - S3)^2/S1^2)*Cos[i]*Sin[i]^2)/(512*a^4*S1) - (315*R^3*(S1 - S3)^2*Cos[i]*Sin[i]^2)/(16*a^4*S1^3) - (315*e^2*R^3*(S1 - S3)^2*Cos[i]*Sin[i]^2)/(8*a^4*S1^3) - (75285*e^4*R^3*(S1 - S3)^2*Cos[i]*Sin[i]^2)/(1024*a^4*S1^3) - (45*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Sin[i]^3)/(8*a^4*S1^2) - (45*e^2*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Sin[i]^3)/(4*a^4*S1^2) - (10755*e^4*R^3*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Sin[i]^3)/(512*a^4*S1^2) + (45*R^3*(S1 - S3)^3*Sin[i]^3)/(16*a^4*S1^4*Sqrt[1 - (S1 - S3)^2/S1^2]) + (45*e^2*R^3*(S1 - S3)^3*Sin[i]^3)/(8*a^4*S1^4*Sqrt[1 - (S1 - S3)^2/S1^2]) + (10755*e^4*R^3*(S1 - S3)^3*Sin[i]^3)/(1024*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*((15*R^4*Cos[i])/(4*a^5*S1) + (15*e^2*R^4*Cos[i])/(4*a^5*S1) + (975*e^4*R^4*Cos[i])/(64*a^5*S1) + (345*R^4*(S1 - S3)*Cos[i]^2)/(8*a^5*S1^2) + (1605*e^2*R^4*(S1 - S3)*Cos[i]^2)/(8*a^5*S1^2) + (17025*e^4*R^4*(S1 - S3)*Cos[i]^2)/(32*a^5*S1^2) - (105*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]^3)/(8*a^5*S1) - (105*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]^3)/(8*a^5*S1) - (6825*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]^3)/(128*a^5*S1) + (105*R^4*(S1 - S3)^2*Cos[i]^3)/(4*a^5*S1^3) + (105*e^2*R^4*(S1 - S3)^2*Cos[i]^3)/(4*a^5*S1^3) + (6825*e^4*R^4*(S1 - S3)^2*Cos[i]^3)/(64*a^5*S1^3) - (525*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^4)/(8*a^5*S1^2) - (1785*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^4)/(8*a^5*S1^2) - (9975*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^4)/(16*a^5*S1^2) - (105*R^4*(S1 - S3)^3*Cos[i]^4)/(8*a^5*S1^4) - (1365*e^2*R^4*(S1 - S3)^3*Cos[i]^4)/(8*a^5*S1^4) - (13125*e^4*R^4*(S1 - S3)^3*Cos[i]^4)/(32*a^5*S1^4) - (15*R^4*(S1 - S3)*Sin[i])/(4*a^5*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) - (15*e^2*R^4*(S1 - S3)*Sin[i])/(4*a^5*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) - (975*e^4*R^4*(S1 - S3)*Sin[i])/(64*a^5*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) + (345*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(8*a^5*S1) + (1605*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(8*a^5*S1) + (17025*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(32*a^5*S1) - (345*R^4*(S1 - S3)^2*Cos[i]*Sin[i])/(8*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) - (1605*e^2*R^4*(S1 - S3)^2*Cos[i]*Sin[i])/(8*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) - (17025*e^4*R^4*(S1 - S3)^2*Cos[i]*Sin[i])/(32*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) + (735*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^2*Sin[i])/(8*a^5*S1^2) + (735*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^2*Sin[i])/(8*a^5*S1^2) + (47775*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^2*Sin[i])/(128*a^5*S1^2) - (105*R^4*(S1 - S3)^3*Cos[i]^2*Sin[i])/(4*a^5*S1^4*Sqrt[1 - (S1 - S3)^2/S1^2]) - (105*e^2*R^4*(S1 - S3)^3*Cos[i]^2*Sin[i])/(4*a^5*S1^4*Sqrt[1 - (S1 - S3)^2/S1^2]) - (6825*e^4*R^4*(S1 - S3)^3*Cos[i]^2*Sin[i])/(64*a^5*S1^4*Sqrt[1 - (S1 - S3)^2/S1^2]) - (525*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(8*a^5*S1) - (1785*e^2*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(8*a^5*S1) - (9975*e^4*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(16*a^5*S1) + (315*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^3*Sin[i])/(2*a^5*S1^3) + (315*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^3*Sin[i])/(2*a^5*S1^3) + (20475*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^3*Sin[i])/(32*a^5*S1^3) + (105*R^4*(S1 - S3)^4*Cos[i]^3*Sin[i])/(8*a^5*S1^5*Sqrt[1 - (S1 - S3)^2/S1^2]) + (1365*e^2*R^4*(S1 - S3)^4*Cos[i]^3*Sin[i])/(8*a^5*S1^5*Sqrt[1 - (S1 - S3)^2/S1^2]) + (13125*e^4*R^4*(S1 - S3)^4*Cos[i]^3*Sin[i])/(32*a^5*S1^5*Sqrt[1 - (S1 - S3)^2/S1^2]) - (345*R^4*(S1 - S3)*Sin[i]^2)/(8*a^5*S1^2) - (1605*e^2*R^4*(S1 - S3)*Sin[i]^2)/(8*a^5*S1^2) - (17025*e^4*R^4*(S1 - S3)*Sin[i]^2)/(32*a^5*S1^2) + (105*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]*Sin[i]^2)/(4*a^5*S1) + (105*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]*Sin[i]^2)/(4*a^5*S1) + (6825*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]*Sin[i]^2)/(64*a^5*S1) - (735*R^4*(S1 - S3)^2*Cos[i]*Sin[i]^2)/(8*a^5*S1^3) - (735*e^2*R^4*(S1 - S3)^2*Cos[i]*Sin[i]^2)/(8*a^5*S1^3) - (47775*e^4*R^4*(S1 - S3)^2*Cos[i]*Sin[i]^2)/(128*a^5*S1^3) + (315*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^2*Sin[i]^2)/(2*a^5*S1^2) + (315*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^2*Sin[i]^2)/(2*a^5*S1^2) + (20475*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*(S1 - S3)^3*Cos[i]^2*Sin[i]^2)/(2*a^5*S1^4) - (315*e^2*R^4*(S1 - S3)^3*Cos[i]^2*Sin[i]^2)/(2*a^5*S1^4) - (20475*e^4*R^4*(S1 - S3)^3*Cos[i]^2*Sin[i]^2)/(32*a^5*S1^4) - (105*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Sin[i]^3)/(4*a^5*S1^2) - (105*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Sin[i]^3)/(4*a^5*S1^2) - (6825*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Sin[i]^3)/(64*a^5*S1^2) + (105*R^4*(S1 - S3)^3*Sin[i]^3)/(8*a^5*S1^4*Sqrt[1 - (S1 - S3)^2/S1^2]) + (105*e^2*R^4*(S1 - S3)^3*Sin[i]^3)/(8*a^5*S1^4*Sqrt[1 - (S1 - S3)^2/S1^2]) + (6825*e^4*R^4*(S1 - S3)^3*Sin[i]^3)/(128*a^5*S1^4*Sqrt[1 - (S1 - S3)^2/S1^2]) - (105*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]*Sin[i]^3)/(8*a^5*S1) - (1365*e^2*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]*Sin[i]^3)/(8*a^5*S1) - (13125*e^4*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]*Sin[i]^3)/(32*a^5*S1) - (315*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]*Sin[i]^3)/(2*a^5*S1^3) - (315*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]*Sin[i]^3)/(2*a^5*S1^3) - (20475*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]*Sin[i]^3)/(32*a^5*S1^3) + (525*R^4*(S1 - S3)^4*Cos[i]*Sin[i]^3)/(8*a^5*S1^5*Sqrt[1 - (S1 - S3)^2/S1^2]) + (1785*e^2*R^4*(S1 - S3)^4*Cos[i]*Sin[i]^3)/(8*a^5*S1^5*Sqrt[1 - (S1 - S3)^2/S1^2]) + (9975*e^4*R^4*(S1 - S3)^4*Cos[i]*Sin[i]^3)/(16*a^5*S1^5*Sqrt[1 - (S1 - S3)^2/S1^2]) + (105*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Sin[i]^4)/(8*a^5*S1^2) + (1365*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Sin[i]^4)/(8*a^5*S1^2) + (13125*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Sin[i]^4)/(32*a^5*S1^2) + (525*R^4*(S1 - S3)^3*Sin[i]^4)/(8*a^5*S1^4) + (1785*e^2*R^4*(S1 - S3)^3*Sin[i]^4)/(8*a^5*S1^4) + (9975*e^4*R^4*(S1 - S3)^3*Sin[i]^4)/(16*a^5*S1^4)) + C44*((-525*R^4*Cos[i])/(4*a^5*S1) + (735*e^2*R^4*Cos[i])/(4*a^5*S1) - (34545*e^4*R^4*Cos[i])/(32*a^5*S1) - (525*R^4*(S1 - S3)*Cos[i]^2)/(4*a^5*S1^2) + (5775*e^2*R^4*(S1 - S3)*Cos[i]^2)/(4*a^5*S1^2) - (23835*e^4*R^4*(S1 - S3)*Cos[i]^2)/(32*a^5*S1^2) + (105*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]^3)/(a^5*S1) + (105*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]^3)/(a^5*S1) + (6825*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]^3)/(16*a^5*S1) + (105*R^4*(S1 - S3)^2*Cos[i]^3)/(4*a^5*S1^3) + (3885*e^2*R^4*(S1 - S3)^2*Cos[i]^3)/(4*a^5*S1^3) - (49035*e^4*R^4*(S1 - S3)^2*Cos[i]^3)/(32*a^5*S1^3) + (105*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^4)/(2*a^5*S1^2) - (1155*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^4)/(2*a^5*S1^2) - (19425*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^4)/(16*a^5*S1^2) + (105*R^4*(S1 - S3)^3*Cos[i]^4)/(4*a^5*S1^4) - (1155*e^2*R^4*(S1 - S3)^3*Cos[i]^4)/(4*a^5*S1^4) - (59745*e^4*R^4*(S1 - S3)^3*Cos[i]^4)/(32*a^5*S1^4) + (525*R^4*(S1 - S3)*Sin[i])/(4*a^5*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) - (735*e^2*R^4*(S1 - S3)*Sin[i])/(4*a^5*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) + (34545*e^4*R^4*(S1 - S3)*Sin[i])/(32*a^5*S1^2*Sqrt[1 - (S1 - S3)^2/S1^2]) - (525*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(4*a^5*S1) + (5775*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(4*a^5*S1) - (23835*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*Cos[i]*Sin[i])/(32*a^5*S1) + (525*R^4*(S1 - S3)^2*Cos[i]*Sin[i])/(4*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) - (5775*e^2*R^4*(S1 - S3)^2*Cos[i]*Sin[i])/(4*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) + (23835*e^4*R^4*(S1 - S3)^2*Cos[i]*Sin[i])/(32*a^5*S1^3*Sqrt[1 - (S1 - S3)^2/S1^2]) - (525*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^2*Sin[i])/(2*a^5*S1^2) + (3255*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^2*Sin[i])/(2*a^5*S1^2) - (34755*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Cos[i]^2*Sin[i])/(8*a^5*S1^2) - (105*R^4*(S1 - S3)^3*Cos[i]^2*Sin[i])/(4*a^5*S1^4*Sqrt[1 - (S1 - S3)^2/S1^2]) - (3885*e^2*R^4*(S1 - S3)^3*Cos[i]^2*Sin[i])/(4*a^5*S1^4*Sqrt[1 - (S1 - S3)^2/S1^2]) + (49035*e^4*R^4*(S1 - S3)^3*Cos[i]^2*Sin[i])/(32*a^5*S1^4*Sqrt[1 - (S1 - S3)^2/S1^2]) + (105*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(2*a^5*S1) - (1155*e^2*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(2*a^5*S1) - (19425*e^4*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]^3*Sin[i])/(16*a^5*S1) - (315*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^3*Sin[i])/(4*a^5*S1^3) + (3465*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^3*Sin[i])/(4*a^5*S1^3) - (62685*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]^3*Sin[i])/(32*a^5*S1^3) - (105*R^4*(S1 - S3)^4*Cos[i]^3*Sin[i])/(4*a^5*S1^5*Sqrt[1 - (S1 - S3)^2/S1^2]) + (1155*e^2*R^4*(S1 - S3)^4*Cos[i]^3*Sin[i])/(4*a^5*S1^5*Sqrt[1 - (S1 - S3)^2/S1^2]) + (59745*e^4*R^4*(S1 - S3)^4*Cos[i]^3*Sin[i])/(32*a^5*S1^5*Sqrt[1 - (S1 - S3)^2/S1^2]) + (525*R^4*(S1 - S3)*Sin[i]^2)/(4*a^5*S1^2) - (5775*e^2*R^4*(S1 - S3)*Sin[i]^2)/(4*a^5*S1^2) + (23835*e^4*R^4*(S1 - S3)*Sin[i]^2)/(32*a^5*S1^2) + (105*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]*Sin[i]^2)/(4*a^5*S1) + (3885*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]*Sin[i]^2)/(4*a^5*S1) - (49035*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*Cos[i]*Sin[i]^2)/(32*a^5*S1) + (525*R^4*(S1 - S3)^2*Cos[i]*Sin[i]^2)/(2*a^5*S1^3) - (3255*e^2*R^4*(S1 - S3)^2*Cos[i]*Sin[i]^2)/(2*a^5*S1^3) + (34755*e^4*R^4*(S1 - S3)^2*Cos[i]*Sin[i]^2)/(8*a^5*S1^3) - (315*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^2*Sin[i]^2)/(4*a^5*S1^2) + (3465*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Cos[i]^2*Sin[i]^2)/(4*a^5*S1^2) - (62685*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*(S1 - S3)^3*Cos[i]^2*Sin[i]^2)/(4*a^5*S1^4) - (3465*e^2*R^4*(S1 - S3)^3*Cos[i]^2*Sin[i]^2)/(4*a^5*S1^4) + (62685*e^4*R^4*(S1 - S3)^3*Cos[i]^2*Sin[i]^2)/(32*a^5*S1^4) - (105*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Sin[i]^3)/(4*a^5*S1^2) - (3885*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Sin[i]^3)/(4*a^5*S1^2) + (49035*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)*Sin[i]^3)/(32*a^5*S1^2) - (105*R^4*(S1 - S3)^3*Sin[i]^3)/(a^5*S1^4*Sqrt[1 - (S1 - S3)^2/S1^2]) - (105*e^2*R^4*(S1 - S3)^3*Sin[i]^3)/(a^5*S1^4*Sqrt[1 - (S1 - S3)^2/S1^2]) - (6825*e^4*R^4*(S1 - S3)^3*Sin[i]^3)/(16*a^5*S1^4*Sqrt[1 - (S1 - S3)^2/S1^2]) + (105*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]*Sin[i]^3)/(4*a^5*S1) - (1155*e^2*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]*Sin[i]^3)/(4*a^5*S1) - (59745*e^4*R^4*(1 - (S1 - S3)^2/S1^2)^(3/2)*Cos[i]*Sin[i]^3)/(32*a^5*S1) + (315*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]*Sin[i]^3)/(4*a^5*S1^3) - (3465*e^2*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]*Sin[i]^3)/(4*a^5*S1^3) + (62685*e^4*R^4*Sqrt[1 - (S1 - S3)^2/S1^2]*(S1 - S3)^2*Cos[i]*Sin[i]^3)/(32*a^5*S1^3) - (105*R^4*(S1 - S3)^4*Cos[i]*Sin[i]^3)/(2*a^5*S1^5*Sqrt[1 - (S1 - S3)^2/S1^2]) + (1155*e^2*R^4*(S1 - S3)^4*Cos[i]*Sin[i]^3)/(2*a^5*S1^5*Sqrt[1 - (S1 - S3)^2/S1^2]) + (19425*e^4*R^4*(S1 - S3)^4*Cos[i]*Sin[i]^3)/(16*a^5*S1^5*Sqrt[1 - (S1 - S3)^2/S1^2]) - (105*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Sin[i]^4)/(4*a^5*S1^2) + (1155*e^2*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Sin[i]^4)/(4*a^5*S1^2) + (59745*e^4*R^4*(1 - (S1 - S3)^2/S1^2)*(S1 - S3)*Sin[i]^4)/(32*a^5*S1^2) - (105*R^4*(S1 - S3)^3*Sin[i]^4)/(2*a^5*S1^4) + (1155*e^2*R^4*(S1 - S3)^3*Sin[i]^4)/(2*a^5*S1^4) + (19425*e^4*R^4*(S1 - S3)^3*Sin[i]^4)/(16*a^5*S1^4))) == 0}