forB31 = (-3*n^2*Csc[eps - i]*(C11*((9*a*e^2)/(16*R) - (9*a*e^4)/(16*R))*Sin[eps] + C20*(-1/4 - (3*e^2)/8 - (15*e^4)/32)*Sin[2*eps] + C22*((533*e^4*Sin[eps])/32 + (533*e^4*Sin[2*eps])/64) + (C31*R*(((-53*e^2)/512 - (39*e^4)/1024)*Sin[eps] + ((-265*e^2)/128 - (195*e^4)/256)*Sin[2*eps] + ((-2385*e^2)/512 - (1755*e^4)/1024)*Sin[3*eps]))/a + (C42*R^2*((-2775*e^4*Sin[eps])/256 - (2775*e^4*Sin[2*eps])/64 - (58275*e^4*Sin[3*eps])/256 - (19425*e^4*Sin[4*eps])/128))/a^2 + (C40*R^2*((5/64 + (25*e^2)/64 + (525*e^4)/512)*Sin[2*eps] + (35/128 + (175*e^2)/128 + (3675*e^4)/1024)*Sin[4*eps]))/a^2))/Omd