Numerical analysis of time fractional three dimensional diffusion equation

The three dimensional diffusion equations were extended to the scope of fractional order derivative. The fractional operator used here is in Caputo sense. The resulting equation was solved using two numerical approaches: The forward in time and central in space method and the Crank-Nicholson method. The stability analysis of both methods was studied, and the study showed that the Crank-Nicholson method is unconditionally stable while the forward method is stable if some conditions are satisfied