% This function calculates the transfer function for imaging
% through an aperture. 
function  phi = pup(x,y)
dg=sqrt(x.^2+y.^2);
phi=0.5*(1-sign(dg-1))-0.5*(1-sign(dg-0.99));