%?平凸透鏡光線追跡
clear;clc;
clear?all;
n=1.5062;%1064nm?波長折射率，k9?玻璃
d=3;?%?透鏡中心厚度
R=100;%?透鏡凸面曲率半徑
Dr=sqrt（R^2-（R-d）^2）;%?透鏡尺寸（最大半徑）
mh=30;
if?Dr>10
	hmax=10;
else
	hmax=Dr;
end
h0=linspace（-hmaxhmaxmh）;
mz=1000;
z0=20;%初始光線與透鏡平面的距離
y=zeros（size（z0））;
theta1=asin（hmax/R）;
theta2=asin（n*hmax/R）;
theta=theta2-theta1;
f=hmax/ta