%% kepler_plot plots the Kepler function, with M = 5 and e = 2.
%
a = 0.0;
b = 10.0;
E = linspace ( a, b, 101 );
fE = kepler ( E );

clf ( );
hold ( 'on' );
plot ( E, fE, 'linewidth', 3 );
plot ( [a,b], [0,0], 'k-', 'linewidth', 3 );
hold ( 'off' );

grid ( 'on' );
xlabel ( '<-- E -->' );
ylabel ( '<-- kepler(E) -->' );
title ( 'Kepler function of E, with M = 5, e = 2' );

print ( '-dpng', 'kepler.png' );
