Function propagate 

pub fn propagate(
    initial: &PropState,
    config: &IntegratorConfig,
    duration: f64,
) -> Vec<PropState>

Propagate an orbit from initial state for a given duration.

Returns a vector of PropState at each time step (including initial state).