This paper applies the Polynomial Least Squares Method (PLSM) to the case of fractional Lane-Emden differential equations. PLSM offers an analytical approximate polynomial solution in a straightforward way. A comparison with previously obtained results proves how accurate the method is.