A note on fourth order method for doubly singular boundary value problems (Q1953195)

Summary: We present a fourth-order finite difference method for the doubly singular boundary value problem \((p(x)y'(x))' = q(x)f(x, y), ~0 < x \leq 1\), with boundary conditions \(y(0) = A\) (or \(y'(0) = 0\) or \(\lim_{x \to 0}p(x)y'(x) = 0\) and \(\alpha y(1) + \beta y'(1) = \gamma\), where \(\alpha(> 0)\), \(\beta(\geq 0)\), \(\gamma\) and \(A\) are finite constants. Here \(p(0) = 0\) and \(q(x)\) is allowed to be discontinuous at the singular point \(x = 0\). The method is based on uniform mesh. The accuracy of the method is established under quite general conditions and also corroborated through one numerical example.