In[]:=
n1=1.7231;(*1:SF10glass(semi-infiniteinputmedium)*)​​n2=2.7039+3.7651;(*2:Tilayer*)​​d2=2.00;(*nm:Tithickness*)​​n3=0.18377+3.4313;(*3:Aulayer*)​​d3=50.00;(*nm:Authickness*)​​n4=1.00028;(*4:air(semi-infiniteoutputmedium)*)
In[]:=
λ0=632.8;(*nm:laserwavelengthinair*)
In[]:=
k0:=
2π
λ0
;​​kx:=k0n1Sin[θDegree];
In[]:=
β1:=
2
(k0n1)
-
2
(kx)
;​​β2:=
2
(k0n2)
-
2
(kx)
;​​β3:=
2
(k0n3)
-
2
(kx)
;​​β4:=
2
(k0n4)
-
2
(kx)
;
In[]:=
r12:=
2
(n2)
β1-
2
(n1)
β2
2
(n2)
β1+
2
(n1)
β2
;​​r23:=
2
(n3)
β2-
2
(n2)
β3
2
(n3)
β2+
2
(n2)
β3
;​​r34:=
2
(n4)
β3-
2
(n3)
β4
2
(n4)
β3+
2
(n3)
β4
;
In[]:=
Rair=
2
Abs
r12+r23Exp[2β2d2]+r34Exp[2(β2d2+β3d3)]+r12r23r34Exp[2β3d3]
1+r12r23Exp[2β2d2]+r12r34Exp[2(β2d2+β3d3)]+r23r34Exp[2β3d3]

;
In[]:=
n4=1.3317;(*4:water(semi-infiniteoutputmedium)*)
In[]:=
Rwater=
2
Abs
r12+r23Exp[2β2d2]+r34Exp[2(β2d2+β3d3)]+r12r23r34Exp[2β3d3]
1+r12r23Exp[2β2d2]+r12r34Exp[2(β2d2+β3d3)]+r23r34Exp[2β3d3]

;
In[]:=
Show[{Plot[Rair,{θ,30,50},PlotRangeAll,PlotStyleRed,AxesLabel{"θ (°)","R"},PlotLegends{"Air"},AxesOrigin{29,0}],Plot[Rwater,{θ,45,70},PlotRangeAll,PlotStyleBlue,PlotLegends{"Water"}]}]
Out[]=
Air
Water