Numerical integration of an orbit and its concomitant variations by recurrent power series