Clamp mean anomaly of closed orbits to -pi..pi
Add orbitStateAtTime, improve tests

Add helper function orbitStateAtTime.

Add a few test cases.

Split tests to three files.