Plots

In[]:=
Clear[t,r,K1,K2,K3]
Constants
In[]:=
aa=12×
-9
10
;(*m*)​​me=9.1093837015×
-31
10
;(*kg*)​​h=1.054571817×
-34
10
;(*Js*)​​ev2J=1.602176634×
-19
10
;(*J*)
In[]:=
r[a_,k1_,k2_,k3_]:=
-((
-ak2

-
ak2

)(
2
k2
-k1k3))+(
-ak2

+
ak2

)(k1k2-k2k3)
-
ak2

(-k1+k2)(k2-k3)+
-ak2

(k1+k2)(k2+k3)
​​K1[En_]:=
2π
2me
h
ev2JEn
​​K2[En_,U1_]:=
2π
2me
h
ev2J(En-U1)
​​K3[En_,U2_]:=
2π
2me
h
ev2J(En-U2)

Case 1: U1 = 1 eV and U2 = 2 eV

In[]:=
U11=1;​​U21=2;
In[]:=
R1[En_]:=
2
Abs[Simplify[r[aa,K1[En],K2[En,U11],K3[En,U21]]]]
In[]:=
Plot[Piecewise[{{0,x<0},{U11,0<x<aa},{U21,x>aa}}],{x,-aa,2aa},PlotRange{{-aa,2aa},{0,5}},Ticks->{{0,{aa,"a"},{2aa,"2a"}},{0,1,2,3,4,5}},ImageSize{{300},{300}},AspectRatio1]​​Rotate[Plot[{R1[x],1-R1[x]},{x,0,5},Ticks{{0,1,2,3,4,5},{0,0.2,0.4,0.6,0.8,1}},ImageSize{{300},{300}},AspectRatio1],0.5π]
Out[]=
Out[]=
GraphicsRow[{Plot[Piecewise[{{0,x<0},{U11,0<x<aa},{U21,x>aa}}],{x,-aa,2aa},PlotRange{{-aa,2aa},{0,5}},Ticks->{{0,{aa,"a"},{2aa,"2a"}},{0,1,2,3,4,5}},ImageSize{{300},{300}},AspectRatio1],​​Rotate[Plot[{R1[x],1-R1[x]},{x,0,5},Ticks{{0,1,2,3,4,5},{0,0.2,0.4,0.6,0.8,1}},ImageSize{{300},{300}},AspectRatio1],0.5π]}]
Out[]=

Failed.
