Off[General::spell]

複素Fresnel反射

関数

Rad[θ_]:=((θ π)/180.)                           (* ラジアン換算 *)
θt[θ_]:=ArcSin[n Sin[θ]]   (* Snellの法則でθtを計算 *)
rs[θ_]:=-((Sin[θ-θt[θ]])/(Sin[θ+θt[θ]]))      (* Fresnelの式 s偏光 *)
rp[θ_]:=((Tan[θ-θt[θ]])/(Tan[θ+θt[θ]]))        (* Fresnelの式 p偏光 *)

入力

n=1.33;                                                        (*ここに,屈折率を代入する*)

計算

p偏光反射率

Plot[Abs[rp[Rad[θ]]],{θ,0,90},PlotRange->All]