Test: Free Particle
wp_sod = QuantumDynamics1D(m)
wp_sod.setup_grid(xmin, xmax, ngrid)
wp_sod.set_potential(np.zeros(wp_sod.xgrid.shape))
psi0 = wp.get_wave_function(wp_sod.xgrid, 0.0)
wp_sod.propagate(dt, nstep, psi0, method="SOD")
anim_numerical = wp_sod.animate(nstep=nstep, delay=50)