A new finite difference method for a class of singular two-point boundary value problems. (Q1399819)

A new finite difference method based on a uniform mesh is given for a class of singular two-point boundary value problems of the form \((x^{\alpha}y')'=f(x,y)\), \(y(0)=A\), \(y(1)=B,\) \(\alpha \in (0,1)\). Assumptions on \(f\) are continuity on \([0,1]\times \mathbb{R}\) and \(\partial f/\partial y\) exists and is non-negative. The method is shown to be fourth order convergent. For \(\alpha=0\), the method reduces to the fourth order Numerov method. Two numerical examples are given.