2D Wavepacket Time Evolution

Simulates a Gaussian wavepacket in a square domain.

This script simulates the time evolution of a quantum mechanical wavepacket on the unit square with Dirichlet boundary conditions. Spatial discretization is with a Chebyshev-tau spectral method and time discretization is with a second order Crank-Nicholson method.