In this paper we solve initial and boundary value problem for non-homogeneous fractional order partial differential equations. Here we use operational matrix approach to construct approximate solutions using Legendre ...