A matrix technique-based numerical treatment of a nonlocal singular boundary value problems

The mathematical modeling of the decisive event of astrophysics, physiology, and many other areas of science and technology witness the involvement of singular boundary value problems. The nonlocal boundary conditions are more informative than local boundary conditions (initial conditions and two-point boundary conditions) to evaluate some mathematical models. This article presents a collocation approach-based matrix technique to approximate the solution of the fusion of a class of singular differential equations subject to nonlocal three-point boundary conditions. The proposed strategy utilizes the truncation of the series expansion of a function belonging to (Formula presented.) in terms of Bernoulli polynomials. It transforms the singular boundary value problems into a set of nonlinear algebraic equations, which can be dealt with by any mathematical software. The Lipschitz condition on an equivalent completely continuous nonlinear operator has been used to prove the convergence analysis of the scheme. Some extremely nonlinear test examples are solved and provided in contrast with the exact solution. These numerical results are also examined against some existing numerical techniques to verify the applicability and significance of the proposed methodology. There are a few numerical examples that are application based but do not have exact solutions. In such cases, residual error norm is employed to measure the accuracy of the numerical strategies. The computational data demonstrate the superiority and validity of the proposed technique over existing numerical approaches.