format compact;
A = 10;				% radius

NptsR=51;
NptsT=51;

Nplots=12; Nrows=3;			% number of columns
plotsols=0;
Ncols=Nplots/Nrows;

r=linspace(0,A,NptsR);
t=linspace(0,2*pi,NptsT);
[R,T]=meshgrid(r,t);
X=R.*cos(T);
Y=R.*sin(T);

fs=8;
figure(3)
clf

for n=1:Nrows
 for m=0:Ncols-1
  subplot(Nrows,Ncols,(m+1)+(n-1)*Ncols)
  zmn=besselzero(m,n,1);
  lambda=(zmn(n)/A)^2;
  B=besselj(m,sqrt(lambda)*R).*cos(m*T);
  surfl(X,Y,B)
  shading interp;colormap gray;
  axis tight;
  drawnow
 end
end