  In this paper the singular Emden-Fowler equation of fractional order is
introduced and a computational method is proposed for its numerical solution.
For the approximation of the solutions we have used Boubaker polynomials and
defined the formulation for its fractional derivative operational matrix. This
tool was not used yet, however, this area has not found many practical
applications yet, and here introduced for the first time. The operational
matrix of the Caputo fractional derivative tool converts these problems to a
system of algebraic equations whose solutions are simple and easy to compute.
Numerical examples are examined to prove the validity and the effectiveness of
the proposed method to find approximate and precise solutions.
