forB11 = -((n^2*Csc[eps - i]*(C11*(a/(2*R) - (a*e^2)/(4*R) - (a*e^4)/(128*R))*Sin[eps] + C20*(-3/4 - (9*e^2)/8 - (45*e^4)/32)*Sin[2*eps] + C22*((3/2 - (15*e^2)/4 + (39*e^4)/32)*Sin[eps] + (3/4 - (15*e^2)/8 + (39*e^4)/64)*Sin[2*eps]) + (C31*R*((-3/64 - (3*e^2)/32 - (717*e^4)/4096)*Sin[eps] + (-15/16 - (15*e^2)/8 - (3585*e^4)/1024)*Sin[2*eps] + (-135/64 - (135*e^2)/32 - (32265*e^4)/4096)*Sin[3*eps]))/a + (C33*R*((225/32 - (675*e^2)/16 + (95175*e^4)/2048)*Sin[eps] + (45/8 - (135*e^2)/4 + (19035*e^4)/512)*Sin[2*eps] + (45/32 - (135*e^2)/16 + (19035*e^4)/2048)*Sin[3*eps]))/a + (C42*R^2*((-15/32 - (15*e^2)/32 - (975*e^4)/512)*Sin[eps] + (-15/8 - (15*e^2)/8 - (975*e^4)/128)*Sin[2*eps] + (-315/32 - (315*e^2)/32 - (20475*e^4)/512)*Sin[3*eps] + (-105/16 - (105*e^2)/16 - (6825*e^4)/256)*Sin[4*eps]))/a^2 + (C40*R^2*((15/64 + (75*e^2)/64 + (1575*e^4)/512)*Sin[2*eps] + (105/128 + (525*e^2)/128 + (11025*e^4)/1024)*Sin[4*eps]))/a^2 + (C44*R^2*((735/16 - (8085*e^2)/16 + (146265*e^4)/128)*Sin[eps] + (735/16 - (8085*e^2)/16 + (146265*e^4)/128)*Sin[2*eps] + (315/16 - (3465*e^2)/16 + (62685*e^4)/128)*Sin[3*eps] + (105/32 - (1155*e^2)/32 + (20895*e^4)/256)*Sin[4*eps]))/a^2))/Omd)