%
% Optical Imaging and Spectroscopy
%
% David J. Brady
% Duke University
% www.opticalimaging.org
%
% 
% pupil function
%
%
function pup=P(x)
A= 1; %cm
pup=(1+sign(A/2-sqrt(x.^2)))/2;
return;