%hm bessel
[x,y]=meshdom(-12:.6:12,-12:.6:12);
r=sqrt(x.^2+y.^2);
z=bessel(0,r);
m=[-45 60];
mesh(z,m)