%
% Optical Imaging and Spectroscopy
%
% David J. Brady
% Duke University
% www.opticalimaging.org
%
%
% OTF function for a unit aperture circular pupil
%
%
function H=otf(u);
H=(2/pi)*(acos(u)-u.*sqrt(1-u.^2));
return