The time-dependent solution of the Schrödinger equation, for a time-dependent potential, can be obtained with the Crank-Nicholson method. Some examples are shown below.
The Gaussian wave packet in this example maintains its shape, even if the harmonic oscillator potential has a nontrivial time dependence; is there a proof for that?