Reverted to the equinoctial elements by Broucke and Cefola (1972) using tan(i/2) rather than sin(i/2), since these seem both more prevalent, and since this was the reference given in the comments in the code. Conversion between keplerian and equinoctial elements implemented, with a couple of tests.
Also added function to check for equivalence for Keplerian elements when degenerated (i and/or e equal to zero), but this functions needs unit tests.
Added a separate module for constants, which will not allow redefinition or deletion of existing constants (but imports to other packages can be overwritten, obviously).
Reverted to the equinoctial elements by Broucke and Cefola (1972) using tan(i/2) rather than sin(i/2), since these seem both more prevalent, and since this was the reference given in the comments in the code. Conversion between keplerian and equinoctial elements implemented, with a couple of tests.
Also added function to check for equivalence for Keplerian elements when degenerated (i and/or e equal to zero), but this functions needs unit tests.
Added a separate module for constants, which will not allow redefinition or deletion of existing constants (but imports to other packages can be overwritten, obviously).