CL_kp_M2Egene — Solution of generalized kepler's equation (ex, ey, M to E)
E = CL_kp_M2Egene(ex,ey,M)
x component of eccentricity vector (1xN)
y component of eccentricity vector (1xN)
mean anomaly [rad] (1xN)
eccentric anomaly [rad] (1xN)
kep = [7078.e3 ; 0.001 ; 0.001 ; %pi/8 ; 0.1 ; 3*%pi/8]; [pos_car,vel_car] = CL_oe_kep2car(kep); cir = CL_oe_car2cir(pos_car,vel_car); ex = cir(2); ey = cir(3); E = CL_kp_M2Egene(ex,ey,3*%pi/8) // Now using CL_kp_M2E : E = CL_kp_M2E(0.001,3*%pi/8) // results slightly differ because convergence criteria is not the same.