The aim of this paper is studying the fractional order nonlinear differential equations arising in electrohydrodynamics flow of a fluid in a circular cylindrical conduit. For this purpose Lucas polynomials are employed via collocation method for reducing the nonlinear singular differential equation of fractional order to an algebraic system. The obtained system is solved by Newton iteration method. The method is computationally simple and attractive.
