A family of variable mesh methods for the estimates of (du/dr) and solution of non-linear two point boundary value problems with singularity

In this article, using three grid points, we discuss variable mesh methods of order two and three for the numerical solution of the nonlinear differential equation u″=f(r,u,u′), 0