function F = circf(kx,ky,diam,x0,y0)
    r = sqrt(kx.^2+ky.^2);
    F = diam^2*myjinc(diam*r).*exp(-1i*2*pi*(x0*kx-y0*ky));

end