In this study, a numerical solution of singular nonlinear differential equations, stemming from biology and physiology problems, is proposed. The methodology is based on the shifted Chebyshev polynomials operational matrix of derivative and collocation. To assess the accuracy of the method, five numerical problems, such as the human head, Oxygen diffusion and Bessel differential equation, were solved. The numerical results were compared with other existed methods in tables for verification. Manuscript profile