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 ...