function [U] = huygens_source(A,Y,X,jsrc,isrc,dy,dx,lambda)

r = sqrt( (isrc*dx-X).^2 + (jsrc*dy-Y).^2 );
k = 2*pi/lambda;
cos_kszi=abs(jsrc*dy-Y)./r;
U = dx/sqrt(lambda)*A./sqrt(r).*exp(-1i*pi/4+1i*k*r).*(1+cos_kszi)/2;
U(jsrc,isrc) = 0; 







