No label defined (Q4372965)

A fourth-order finite difference method is used to solve the singular two-point boundary value problem \[ \sum^2_{k=0} (-1)^k {d^k \over dx^k} \left[p_k(x) {d^ku \over dx^k} \right]= F(x,u),\;0<x\leq 1 \] \[ \sum^3_{k=0} \alpha_{\ell_k} \left({d^ku \over dx^k} \right)_{x \to 0} =f_\ell, \quad \sum^3_{k=0} \beta_{\ell_k} \left({d^ku \over dx^k} \right)_{x=1} =g_\ell,\;\ell= 0,1. \] Newton's method is used to linearize the problem. For each such linear equation the method of series expansion in a small neighbourhood of the point \(x=0\) is employed and the difference method is used in the rest of the interval. The consistency and stability of the method is studied and some numerical examples are given to illustrate the method.