Scattering over Potential Step

V[x_,V0_:1]:=V0*UnitStep[x];

ψ[x_,En_,V0_]:=Module[
{k=Sqrt[2En],q=Sqrt[2(En-V0)],R,T},
R=(k-q)/(k+q);
T=1+R;
((Exp[I k x]+R Exp[-I k x])UnitStep[-x]+
T Exp[I q x]UnitStep[x])
]

Manipulate[
Plot[
With[{val=ψ[x,En,V0]},
{Abs[val]^2,Re[val],Im[val],En,V[x,V0]}]//Evaluate,
{x,-10,10},
PlotStyle->{RGBColor[1,.47,0],RGBColor[.38,.61,.51],ColorData["HTML","SlateBlue"],Gray,Lighter[Orange,.75]},
Filling->{5->0},
FillingStyle->Lighter[Orange,.75],ImageSize->{500,350},PlotRange->{-2,5.5}],
{{En,0.5,"particle energy"},0,5,Appearance->"Labeled"},
{{V0,1,"potential energy step"},0.5,5,Appearance->"Labeled"},
SaveDefinitions->True
]

THIS NOTEBOOK IS THE SOURCE CODE FROM

"Scattering over Potential Step" from The Wolfram Demonstrations Project
 http://demonstrations.wolfram.com/ScatteringOverPotentialStep/

Contributed by: Mahn-Soo Choi

A full-function Wolfram Mathematica 6 system is required to edit this notebook.
GET WOLFRAM MATHEMATICA 6 »

© The Wolfram Demonstrations Project & Contributors   |   Terms of Use   |   Make a new version of this Demonstration »