transfer.m

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

function f=transfer(x0,x1,N,P,mult,E);
%> The file <transfer.m> integrates the 1D Schrodinger equation 
%> for an 'arbitrary' potential using the Numerov method.
%> It produces the transfer matrix from x0 to x1
%> as a function of energy (to satisfy boundary condition),
%> Call: transfer(x0,x1,N,'pot(X)',M,E);
%> Input: (x0,x1) = end points. Interchanging x0,x1 gives inverse matrix.
%> 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.
%> Output: f = [u(x1); Du(x1); v(x1); Dv(x1)] = transfer matrix 
%> Fundamental solutions (u,v) defined by 
%> u(x0)=1, Du(x0)=0, v(x0)=0, Dv(x0)=1;
%> size of f = [4,length(E)]; 
%>
%>