evolut.m

来自「求解量子力学的薛定谔方程」· M 代码 · 共 12 行

%> The file <evolut.m> calculates the evolution of a wave packet on 
%> the unit interval using an iterative method for the propagation. 
%> Algorithm = Visscher scheme. System = particle in a box with
%> boundary conditions: u=0  at x=0,1 and an additional potential.
%> Initial wave function = Gaussian multiplied by plane wave.
%> Reference: Visscher: Computers in Physics, nov/dec 1991. 
%> There are a number of choices for the  potential and you can easily define
%> more of them. Note that the results do tend to give little more 
%> than random waves sloshing about unless the spectrum is almost harmonic. 
%>
%>