  Nonadiabatic behavior of metastable systems modeled by anharmonic
Hamiltonians is reproduced by the Fokker-Planck and imaginary time Schrodinger
equation scheme with subsequent symplectic integration. Example solutions
capture ergodicity breaking, reassure the H-theorem of global stability [M.
Shiino, Phys. Rev. A, Vol 36, pp. 2393-2411 (1987)], and reproduce spatially
extended response under alternate source fields.
