A Numerical Approach for Solving a Class of a Singular Boundary Value Problems Arising in Physiology

Abstract

In this paper, two numerical schemes for finding approximate solutions of singular two-point boundary value problems arising in physiology are presented. While the main ingredient of both approaches is the employment of cubic B-splines, the obstacle of singularity has to be removed first. In the first approach, L'Hopital's rule is used to remove the singularity due to the boundary condition (BC) $y'(0) = 0$. In the second approach, the economized Chebyshev polynomial is implemented in the vicinity of the singular point due to the BC $y(0) = A$, where $A$ is a constant. Numerical examples are presented to demonstrate the applicability and efficiency of the methods on one hand and to confirm the second order convergence on the other hand.