We present a numerical method for obtaining an approximate polynomial solution for nonlinear singular differential equations using the Bessel matrix method based on collocation points. This method provides a system of linear equations with unknown Bessel coefficients. In addition, it presents high accuracy of solutions when exact solutions are polynomials. The reasoning and efficiency of the method are demonstrated via numerical examples. In addition, the method is compared with other techniques

Keywords : Bessel Matrix, Nonlinear, Singular, Numerical Solutions