Couple of the variational iteration method and fractional-order Legendre functions method for fractional differential equations.

We present a new numerical method to get the approximate solutions of fractional differential equations. A new operational matrix of integration for fractional-order Legendre functions (FLFs) is first derived. Then a modified variational iteration formula which can avoid "noise terms" is constructed.