Von Neumann stability analysis with SymPy: the leapfrog Courant limit

Afterwards you can derive a finite difference scheme's amplification factor with SymPy, find its stability limit, and expand its numerical dispersion.