The anomaly system should be described in `Orbit` docstring, and it should be possible to give `mean_anomaly` as input to the init function