Added kepler_pos to numerically calculate position(t) assuming a keplerian orbit