Numerical Solution of some Mathematical Models in Physiology

For a class of singular boundary value arising in Physiology, a finite difference scheme is proposed. Based on the asymptotic expansion of the solution in the neighborhood of the singular points, smoothing variable substitutions are introduced in order to improve the convergence order of the finite diference methods. Numerical results are presented and discussed.