MATLAB Answers


Plotting E - e*sin[E] = M as E vs M, when given e as 5 different values [0 0.25 0.5 0.75 1]?

Asked by Daniel
on 12 Sep 2013
Latest activity Answered by Youssef Khmou
on 30 Sep 2013

For a question in my homework , we are suppose to first find the fzeros of this function of E when given any values of e, and M. I already figured this out by the following:

function [ E,z ] = kepler( e,M )
%UNTITLED Summary of this function goes here
%   Detailed explanation goes here
f=@(E) E-e*sin(E)-M;

But I cannot for the life or me figure out how I would plot this E vs M when only given a value of e. If anyone can help me out that would be greatly appreciated :)


Log in to comment.

2 Answers

Answer by Walter Roberson
on 12 Sep 2013

The usual way: loop over all the potential values of M for a fixed e, recording the outcome each time. Plot the potential values of M on one axis, and the recorded E values on another axis.

hint: ndgrid(), arrayfun(), and surf()


Show 1 older comment

Another hint: to loop over values of M, use the 'for' statement. See the examples contained in:

doc for

Also they have to be shown as five different curves on the same graph, but I know if you plotted each separately for e value given then E vs M and use "Hold on" command that would work.

I have the same problems too, how can you graph E vs. M for e of values [0 0.25 0.5 0.75 1] with 5 curves in 1 figure?

Log in to comment.

Answer by Youssef Khmou
on 30 Sep 2013

Daniel, i find this issue not clear at all, besides you have to to give some Physical explanations about the constants, M as mean anomaly ,e as eccentricity :

function [ E] = kepler( e,M )
f=@(E) E-e*sin(E)-M;
e=[0 0.25 0.5 0.75 1];
for x=1:N
     for y=1:n
figure, plot(M,E), grid on,


Log in to comment.

Discover what MATLAB® can do for your career.

Opportunities for recent engineering grads.

Apply Today