function F = rectf(kx,ky,Lx,Ly,x0,y0)
    F = Lx*Ly*mysinc(Lx*kx).*mysinc(Ly*ky).*exp(-1i*2*pi*(x0*kx-y0*ky));

end