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