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