We use B-spline functions to develop a numerical method for computing approximations to the solution of non-linear singular boundary value problems associated with physiology science. The original differential equation is ...