forA32 = -((n*Csc[eps - i]*(C20*(-1/2 - (3*e^2)/4 - (15*e^4)/16)*Sin[2*eps] + C22*(((7*e)/2 - (123*e^3)/16)*Sin[eps] + ((7*e)/4 - (123*e^3)/32)*Sin[2*eps]) + (C42*R^2*(((-45*e)/32 + (15*e^3)/256)*Sin[eps] + ((-45*e)/8 + (15*e^3)/64)*Sin[2*eps] + ((-945*e)/32 + (315*e^3)/256)*Sin[3*eps] + ((-315*e)/16 + (105*e^3)/128)*Sin[4*eps]))/a^2 + (C44*R^2*(((12495*e^2)/16 - (78645*e^4)/16)*Sin[eps] + ((12495*e^2)/16 - (78645*e^4)/16)*Sin[2*eps] + ((5355*e^2)/16 - (33705*e^4)/16)*Sin[3*eps] + ((1785*e^2)/32 - (11235*e^4)/32)*Sin[4*eps]))/a^2 + (C40*R^2*((5/32 + (25*e^2)/32 + (525*e^4)/256)*Sin[2*eps] + (35/64 + (175*e^2)/64 + (3675*e^4)/512)*Sin[4*eps]))/a^2))/(Omd*(1 + (2*omd)/(3*n) + (2*Omd*Cos[eps - i])/(3*n))))