Kepler orbit

Solves Kepler’s equation and computes the period, perihelion, aphelion, vis-viva speed, total energy and a table of orbit samples, then checks energy conservation with a numerical integration.