numerov2.m

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

function f=numerov2(x0,x1,N,P,mult,E,Y0,DY0);
%>  The file <numerov2.m> integrates the 1D Schrodinger equation 
%>  for an 'arbitrary' potential using the Numerov method.
%>  It produces the function values in the integration interval (x0,x1)
%>  for a row vector of energy values.
%>  Call: numerov2(x0,x1,N,'pot(X)',M,E,Y0,DY0);
%>  Input: (x0,x1) = 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.
%>  X0,X1 = starting values in N = x(1),x(2), row vectors size(E).
%>  Output: matrix of dimension (N, length(E));
%>
%>