Return x,y,z values of a shape
[x,y,z] = ieShape(type,nSamp,varargin)
We calculate shapes like circles, ellipses, maybe other stuff later.
Example:
Positions of the Airy disk zero crossing for an f# 5.6 lens at 550 nm.
The units of x,y are meters
fNumber = 2.0; wavelength = 550;
radius = (2.44*fNumber*wavelength*10^-9);
[x,y] = ieShape('circle',200, radius); plot(x,y,'.');
axis equal; grid on