Numerical solution of fractional differential equations using the generalized block pulse operational matrix

The Riemann–Liouville fractional integral for repeated fractional integration is expanded in block pulse functions to yield the block pulse operational matrices for the fractional order integration. Also, the generalized block pulse operational matrices of differentiation are derived. Based on the above results we propose a way to solve the fractional differential equations. The method is computationally attractive and applications are demonstrated through illustrative examples.