In this study, a new approach based on B-spline scaling functions and Galerkin method has been developed for solving nonlinear singular boundary value problems associated with astrophysics science. The method is computationally attractive and application is demonstrated with illustrative example.
