This study deals with elaborating a novel framework for acquiring the approximate solution of the nonlinear singular boundary value problems of fractional order arising in biology. The quesilinearization technique is applied to reduce the given nonlinear problem to a sequence of linear problems. We modify the resulting set of fractional order differential equations at the singular point then treat this set of boundary value problems by using Bernoulli wavelets via spectral collocation approach. The method is computationally attractive, and applications are demonstrated through illustrative examples.
