%
% Optical Imaging and Spectroscopy
%
% David J. Brady
% Duke University
% www.opticalimaging.org
%
%
function otF = otf(u)
%otf function
u=abs(u);
otF=(2/pi)*real(acos(u)-u.*sqrt(1-u.^2));
return