function  [z,y] = drawlens(ap,th,cv,zstart)
%
% DRAWLENS(AP,TH,CV)
%
%

if (nargin<4)
   zstart = 0;
end

u = linspace(-ap,ap,60);

z = [ sag(u,cv(1)) th+sag(-u,cv(2)) sag(-ap,cv(1)) ] + zstart;
y = [ u  -u -ap ];
if (nargout==0)
  plot(z,y)
end

Published with MATLAB® 7.10