stepp.m

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

%> The file <stepp.m> calculates the transmission coefficient for a 
%> potential step using the Numerov algorithm.
%> The scattering potential is constant V = V0 > 0 in x < 0, 
%> constant V = 0 in x > 1. 
%> The potential in [0,1] is input as a string like 'cos(pi*X/2)'. 
%> The particle comes in from x = + infinity, we start the integration
%> from the transmitted wave in x = 0, wave number wave2=sqrt(2*(E-V0)), 
%> and integrate to x = 1, wave number wave1=sqrt(2*E), where we adapt the 
%> incoming and reflected waves and renormalize to find the transmission
%> coefficient. This is displayed as a function of E-V. 
%> For comparision the result for a 'pure step' is displayed first.
%> Reference: Messiah Ch III.
%>
%>