The Lane–Emden–Fowler-type equations along with boundary conditions form nonlinear singular boundary value problems. Analytical solutions to this type of problem are not easy, which draws the attention of the researchers toward their numerical solutions. Many methods have been used to provide numerical solutions to this type of problem. In this study, we have proposed the quasilinear Bessel polynomial collocation method (QBPCM) for the solution of these problems. The accuracy and efficiency of the QBPCM have been demonstrated by solving several nonlinear singular models that arise in real-life problems. It has been shown that the numerical results obtained by the proposed method have excellent agreement with the exact solutions and better accuracy than other methods.