平凸透镜光线追迹

% 平凸透镜光线追迹
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