drawsys

%%function drawsys(ap,th,cv,lp)
% draws a lens system

ts = cumsum(th);
for k=lp
    drawlens(sa(k),th(k),cv(k:k+1),ts(k));
    hold on;
end
axis equal;
hold off;
clear k ts;