Numerical Solution of Some Fractional Diffusion Equations

This Demonstration shows numerical solutions of the fractional diffusion equation by means of weighted average methods (or methods). The boundary conditions specify that the solution equals zero at and . Four different initial conditions can be chosen. When the weight parameter equals 1/2, the numerical method is the fractional Crank–Nicolson method. When the "normal solution" checkbox is checked, the normal diffusion solution is also plotted.


  • [Snapshot]
  • [Snapshot]
  • [Snapshot]


The fractional diffusion equation is , where the is the Riemann–Liouville derivative of order , , is solved numerically by means of the fractional difference weighted-average methods (or methods) discussed in [1]. When , the equation is just the normal diffusion equation and these methods become normal weighted-average methods [2]. Depending on the value of the weight factor , the method can be fully implicit when , or fully explicit when . The fractional Crank–Nicolson method for , or for other values of , is a generic -weighted-average method. It is also determined whether the stability criterion of equation (35) of [1] is also satisfied. The size of the timestep and the number of timesteps can be chosen. The diffusion constant and the space discretization are fixed to 1 and with , respectively. When the "normal solution" checkbox is checked, the normal diffusion solution () is also plotted (orange points). The broken lines are to guide the eye. The values of the numerical solutions are shown when the mouse pointer is over the points.
[1] S. B. Yuste, "Weighted Average Finite Difference Methods for Fractional Diffusion Equations," Journal of Computational Physics, 216, 2006 pp. 264–274.
[2] K. W. Morton and D. F. Mayers, Numerical Solution of Partial Differential Equations, Cambridge: Cambridge University Press, 1994.
    • Share:

Embed Interactive Demonstration New!

Just copy and paste this snippet of JavaScript code into your website or blog to put the live Demonstration on your site. More details »

Files require Wolfram CDF Player or Mathematica.

Mathematica »
The #1 tool for creating Demonstrations
and anything technical.
Wolfram|Alpha »
Explore anything with the first
computational knowledge engine.
MathWorld »
The web's most extensive
mathematics resource.
Course Assistant Apps »
An app for every course—
right in the palm of your hand.
Wolfram Blog »
Read our views on math,
science, and technology.
Computable Document Format »
The format that makes Demonstrations
(and any information) easy to share and
interact with.
STEM Initiative »
Programs & resources for
educators, schools & students.
Computerbasedmath.org »
Join the initiative for modernizing
math education.
Step-by-Step Solutions »
Walk through homework problems one step at a time, with hints to help along the way.
Wolfram Problem Generator »
Unlimited random practice problems and answers with built-in step-by-step solutions. Practice online or make a printable study sheet.
Wolfram Language »
Knowledge-based programming for everyone.
Powered by Wolfram Mathematica © 2018 Wolfram Demonstrations Project & Contributors  |  Terms of Use  |  Privacy Policy  |  RSS Give us your feedback
Note: To run this Demonstration you need Mathematica 7+ or the free Mathematica Player 7EX
Download or upgrade to Mathematica Player 7EX
I already have Mathematica Player or Mathematica 7+