In the early 17th century Johannes Kepler developed the equation, now often referred to as kepler's equation for the solution of the mean anomaly (M) of a planet in an elliptical orbit as a function of it's eccentricity (e) and it's eccentric anomaly (E).
E - e*sin[E] = M
If you are interested in exactly what this equation represents please feel free to look deeper into orbital mechanics. Curtis's Orbital Mechanics for Engineers is an excellent place to start. For this assignement however the important thing to notice is that given E and e, solving for M is quite trivial. Given M and e, solving for E can prove quite difficult. In fact this problem is considered transcendental, meaning E can not be represented as a finite set of algebraic terms of e and M. As such this problem is almost always solved using some numerical method.
For this assignment I would like you create a function called kepler which solves the kepler equation for E given e and M (your function should accept M and e as inputs, and return E). The actual solution should be performed using the MATLAB built in function fzero() to solve the following form of Kepler's equation:
0 = E - e*sin[E] – M
We use this form because fzero solves for roots of the input function. Recall that fzero() is a function function, and that it that it expects that the function argument which you pass it will only accept one argument. This means that fzero(f(x)) is acceptable, but fzero(f(x,y)) is not. As such you might have to use the shared workspace property of nested functions to your advantage to complete this assignment.
Please plot E vs. M for e of values [0 0.25 0.5 0.75 1]. Please include all 5 curves in 1 figure axis. Please label axis and include a legend.