Optical Imaging and Spectroscopy

David J. Brady
Duke University
www.opticalimaging.org

Figures 2.7 and 2.8

thetaPerp1=Pi/10;
thetaPerp2=-Pi/10;
n1=1.0;
n2Max=1.8;
thetaAMin=thetaPerp1-Pi/2;
thetaAMax=thetaPerp1+ArcSin[(n2Max/n1)Sin[thetaPerp2-thetaPerp1+ArcSin[n1/n2Max]]];
thetab[thetaa_,n2_]:=thetaPerp2+ArcSin[(n2/n1) Sin[ArcSin[(n1/n2)Sin[thetaa-thetaPerp1]]+thetaPerp1-thetaPerp2]];

figure 2.7

Plot[{-((180 thetab[(1/180) (-x) π,1.1`])/π),-((180 thetab[(1/180) (-x) π,1.3`])/π),-((180 thetab[(1/180) (-x) π,1.5`])/π),-((180 thetab[(1/180) (-x) π,1.7`])/π)},{x,-((180 thetaAMax)/π),-((180 thetaAMin)/π)},Frame->True,Axes->False,FrameLabel->{"Incident Angle","Transmitted  angle"},BaseStyle->{FontSize->14}]

f2 flint glass

b1=1.34533359;
b2=2.09073176 10^(-1);
b3=9.37357162 10^(-1);
c1=9.97743871 10^(-3);
c2=4.70450767 10^(-2);
c3=1.11886764 10^2;
n[lambda_]:=Sqrt[1+b1 lambda^2/(lambda^2-c1)+b2 lambda^2/(lambda^2-c2)+b3 lambda^2/(lambda^2-c3)];

Plot[n[(l/1000)],{l,500,700},Frame->True,Axes->False,FrameLabel->{"wavelength (nm)","Index of Refraction"},BaseStyle->{FontSize->14}]

figure 2.8

Plot[-((180 thetab[0,n[(l/1000)]])/π),{l,400,800},Frame->True,Axes->False,FrameLabel->{"wavelength (nm)","Transmitted angle"},BaseStyle->{FontSize->14}]

sf14 glass

b1=1.69182538;
b2=0.285919934;
b3= 1.12595145;
c1= 0.0133151542;
c2= 0.0612647445;
c3= 118.405242;
n[lambda_]:=Sqrt[1+b1 lambda^2/(lambda^2-c1)+b2 lambda^2/(lambda^2-c2)+b3 lambda^2/(lambda^2-c3)];