Please answer part a as pen and paper answer and part b,c,d as MATLAB code

5.9. In celestial mechanics, Kepler's equation M = E - esin(E) relates the mean anomaly M to the eccentric anomaly E of an elliptical orbit of eccentricity e, where 0