numerov1.m

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

function f=numerov1(x0,x1,N,P,mult,E,Y0,DY0);
%> The file <numerov1.m> integrates the 1D Schrodinger equation 
%> for an 'arbitrary' potential using the Numerov method.
%> It produces the function values and derivatives at the second end 
%> point as a function of energy (to satisfy boundary condition),
%> it does not keep the solution function.
%> Call: numerov1(x0,x1,N,'pot(X)',M,E,Y0,DY0);
%> Input: (x0,x1) = end points of integration interval.
%> N = number of steps (N subintervals of length h).
%> P = 'pot(X)' if the potential function is stored in <pot.m>.
%> M = multiplying factor for potential, incl sign.
%> E = energy lattice as a row vector.
%> Y0,DY0 = starting values in x0 , row vectors length(E). 
%> Output: f = [Y;DY] = final values of function and derivative in x1,
%> size = [2,length(E)];
%>
%>