transm.m

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

%> The file <transm.m> calculates the transmission coefficient in 1D 
%> potential scattering using the Numerov algorithm for integration. 
%> The scattering potential is nonzero only on [0,1]. 
%> A lattice of values for the energy E is chosen in the CW.
%> The particle comes in from x = + infinity, we start the integration
%> from the transmitted wave in x = 0, wave number sqrt(2*E). 
%> At x = 1 we adapt the incoming and reflected waves and renormalize 
%> to find the transmission coefficient, which is displayed as a 
%> function of wave number. Forbidden wave number bands are displayed
%> for a choice of unit cell spacings.
%>
%>