
No new physics here, but this may serve as a benchmark for wavepacket dynamics in three dimensions and also as a reference implementation for iso-surface plots in three dimensions.
In higher dimensionality, the physics of (uncoupled!) harmonic oscillators remains essentially unchanged. This case represents an example of separable potentials where the multi-dimensional wavefunctions can be written as products of one-dimensional wavefunctions.
In the animation, we see an example of a squeezed state in three dimensions. While the widths of the Gaussian (bell-) shape of the wavepacket oscillate in time, its center (expectation value of position vector) follows a Kepler ellipsis.