wkb.m

求解量子力学的薛定谔方程

function W = wkb(X,Y,Emax,N);
%> This file finds the semiclassial WKB energy eigenvalues.
%> Call: f = wkb(X,Y,Emax,N)
%> Input: X = lattice in a space interval, Y=pot(X), both row vectors;
%> 	Emax = maximal energy, can be set = 0 or a little bit larger to
%> get a weakly bound state which ends up with a positive energy
%> due to the errors in method and calculation. Emin is automatically 
%> set to min(Y)! 
%>   N = number of points in energy lattice.
%> Output: W = row vector of eigenvalues in increasing order.
%> This procedure should be used with caution unless there is a
%> single minimum of the potential.
%> Reference: For the theory see Messiah Ch 6.
%>
%>