A numerical perturbation expansion method is developed, analysed and implemented for the numerical solution of a second-......