trf1.m

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

function f=trf1(X,V,E);
%> This file calculates the transfer matrix for 1D Schroedinger DE
%> for an interval with a constant potential 
%> Call:  trf1(X,V,E)
%> Input: X=[X1, X2] is the interval as a row 2-vector, X1 < X2
%> V = potential, scalar
%> E = energy lattice, row vector, INCREASING ARGUMENT
%> Output: a column 4-vector [u;Du;v;Dv] at X2
%> where the fundamental solutions (u,v) are defined by
%> u(X1) = Dv(X1) = 1, Du(X1) = v(X1) = 0.
%>
%>