Numerical solution of second order Painlevé differential equation