Integrablization of time fractional PDEs

Here we present a numerical method for the solution of time fractional partial differential equations (fPDEs). The method is based on constructing a sequence of integrable approximations to the fPDE, such that in the appropriate limit the original fractional PDE is obtained. Convergence results are obtained in general and examples are given for simple cases such as the fractional heat equation and the fractional diffusion-wave equation.