doublet1

Contents

lens specification

rd = [0 30.658123 -26.733948 -25.366250  -57.832251 0];
th = [0  2 0.1 1 48.484530 0];
rn = [1 1.5168 1 1.728277 1 1];
sa = [ 0 5.1 5.1 5.1 5.1 0];
cv = curvature(rd);
lp = [2 4];
py = 5;
lambda = 0.58756e-3;
puc =  tan(2*pi/180);
fprintf('puc %g\n',puc);
puc 0.0349208

draw system

drawsys;

paraxial raytrace

pa = [py 0];
pc = [0 puc];
axial = parax(pa,cv,th,rn);
chief = parax(pc,cv,th,rn);
fprintf('\n%3s %12s %12s %12s %12s\n','srf','py','pu','pyc','puc');
for m = 1:length(th)
    fprintf('%3d %12.6f %12.6f %12.6f %12.6f\n',m-1,axial(m,:),chief(m,:));
end
srf           py           pu          pyc          puc
  0     5.000000     0.000000     0.000000     0.034921
  1     5.000000    -0.055567     0.000000     0.023023
  2     4.888866    -0.178792     0.046045     0.034031
  3     4.870986    -0.022533     0.049448     0.020512
  4     4.848453    -0.100000     0.069960     0.034569
  5     0.000000    -0.100000     1.746038     0.034569

Wavefront aberrations

W = ford(axial,chief,cv,th,rn,lambda);
L =

    0.1746


W =

         0         0         0         0         0
    1.0365    0.8877    0.1901    0.2392    0.1024
   19.3395   -6.9106    0.6173    0.4740   -0.0847
  -23.6190    8.1736   -0.7071   -0.5691    0.0985
    3.0315   -2.2004    0.3993    0.2942   -0.1068
   -0.2116   -0.0497    0.4996    0.4383    0.0095