function y = kepler(E, Me, e)

y = E - e*sin(E) - Me;

end

