mhankel

function u = mhankel(rho,W,r)

wopd = W(rho);
u=exp(1i*2*pi*wopd).*besselj(0,2*pi*r*rho).*rho;