%
%  w10besselplot.m -- (djm - 14 mar 06)
%

th = 0:pi/40:2*pi;

bz = 1.84118;
r = -bz:bz/100:bz;

[tg,rg] = meshgrid(th,r);
[xg,yg] = pol2cart(tg,rg);

wave = besselj(1,rg).*cos(tg+pi/6);

figure(2);  clf
subplot(2,1,2)
surf(xg,yg,wave);  shading flat;

subplot(2,1,1)
plot(r,real(besselj(1,r)),'b')
axis([-bz bz -0.6 0.6])
title('\bf first-mode')

bz = 5.33144;
r = -bz:bz/100:bz;

[tg,rg] = meshgrid(th,r);
[xg,yg] = pol2cart(tg,rg);

wave = besselj(1,rg).*cos(tg+pi/6);

figure(1);  clf
subplot(2,1,2)
surf(xg,yg,wave);  shading flat;

subplot(2,1,1)
plot(r,real(besselj(1,r)),'b')
axis([-5.34 5.34 -0.6 0.6])
title('\bf second-mode')