Solving the time-dependent Schrödinger equation using finite difference methods

We solve the time-dependent Schrequation in one and two dimensions using the finite difference approximation. The evolution is carried out using the method of lines. The illustrative cases include: the particle in a box and the harmonic oscillator in one and two dimensions. As non-standard examples we evolve two solitons and show the time-dependent solitonic behavior in one dimension and the stabilization of an atomic gas model in two dimensions. The codes used to generate the results in this manuscript are freely available under request, and we expect this material could help students to have a better grasp of the solution of partial differential equations related to dynamical systems.