The purpose of this study is to propose a computationally effective algorithm for the numerical evaluation of a fractional-order system of singular Lane-Emden type equations arising in physical problems. The fractional operator considered is in the sense of the Liouville-Caputo derivative. The presented matrix collocation method is based upon a combination of the quasilinearization method (QLM) and the shifted Legendre functions (SLFs) and is called QLM-SLFs method. By applying first the QLM to the nonlinear underlying system, we get a family of linear equations. Hence, a spectral matrix collocation scheme relied on the SLFs is designed to solve the resulting sequence of linear system of equations at very few iterations. The uniform convergence of the shifted Legendre expansion series solution is established. To illustrate the effectiveness of the proposed QLM-SLFs technique in the present paper, three test examples are carried out. The applicability and validity of the proposed method are testified through comparisons with the outcomes of other existing procedures in the literature. The proposed QLM-SLFs method is efficient and easy to implement. The approximation obtained by the method also converges quickly to the solutions of the underlying model problem. In comparison with available existing computational procedures, the QLM-SLFs approach shows that the use of Legendre functions together with QLM provides solutions with high accuracy and exponential convergence rate.
Read full abstract