def Efn(Ms,eccs):
Ms = np.atleast_1d(Ms)<EOL>eccs = np.atleast_1d(eccs)<EOL>unit = np.floor(Ms / (<NUM_LIT:2>*np.pi))<EOL>Es = EFN((Ms % (<NUM_LIT:2>*np.pi)),eccs)<EOL>Es += unit*(<NUM_LIT:2>*np.pi)<EOL>return Es<EOL>
works for -2pi < Ms < 2pi, e <= 0.97