added a way to perform numerical propagation in either cartesian parameters or equinoctial parameters