triplet129

Contents

define lens
show specifications
draw lens

define lens

% radius of curvature
rd = [0 24.110 215.090 -94.81 23.76 0 104.5 -63.89];
idx = find(rd);
cv = zeros(size(rd));
cv(idx) = 1./rd(idx);
th = [0 3.7 4.66 1.6 2.37 6.76 3.5 0];
rn = [1 1.613 1 1.596 1 1 1.613 1];
lp = find(rn>1);
sa = [0 11.4 11.4 7 7 6.6 10.7 10.7];
ast = 6; % location of aperture stop

show specifications

fprintf('%4s%12s%10s%8s%8s\n','n','rd','th','index','sa');
n = length(rd);
for k=1:n
    fprintf('%4d%12.4f%10.4f%8.4f%8.2f\n',k-1,rd(k),th(k),rn(k),sa(k));
end

   n          rd        th   index      sa
   0      0.0000    0.0000  1.0000    0.00
   1     24.1100    3.7000  1.6130   11.40
   2    215.0900    4.6600  1.0000   11.40
   3    -94.8100    1.6000  1.5960    7.00
   4     23.7600    2.3700  1.0000    7.00
   5      0.0000    6.7600  1.0000    6.60
   6    104.5000    3.5000  1.6130   10.70
   7    -63.8900    0.0000  1.0000   10.70

draw lens

drawsys(sa,th,cv,lp);
axis([-1 25 -13 13]);