function y=RSkernel(x,x2,z2)
y=exp(j*2*pi*sqrt((x2-x).^2+z2^2))*z2 ./((x2-x).^2+z2^2);