We aim to present numerical methods based on Legendre wavelets and quasilinearization technique for fractional Lane-Emden type equations and time-fractional Fisher’s equation. The Lane-Emden equation is a second order singular non-linear ordinary differential equation, which is useful for modelling many astrophysical phenomena such as the distribution of stars in star clusters and star formation in molecular clouds. The Fisher’s equation is a non-linear reaction-diffusion equation that models the spread of mutant genes in a population. We start with a brief discussion of the purpose of studying fractional differential equations. Then some practical aspects of wavelets are explained. We also give introductory definitions and properties of fractional calculus and Legendre wavelets. Using Legendre wavelets and quasilinearization technique, we derive numerical methods for fractional Lane-Emden type equations and time-fractional Fisher’s equation. Moreover, the convergence analysis of both methods is studied. Some problems are solved to evaluate the efficiency of the proposed methods. Test problems show that the proposed methods are very effective.