points for a direct solution of stiff second order ODEs. Power series approximation method
was adopted as basis function using interpolation and collocation techniques at two
selected grid points. We investigate some basic properties of our method and found it to
converge. Stability region plotted shows that the method is A-stable and can handle stiff
problems. The method was then tested on some highly stiff problems and numerical results …