r=0:0.01:1;
[XX,YY]=meshgrid(-1,0.01,1);
[Q,R]=cart2pol(XX,YY);
R(find(R>1))=NaN;
y=inline('besselj(0,x)','x');
x(1)=fzero(y,3)
meshc(XX,YY,y(x(1)*R))
for k=1:3
x(k+1)=fzero(y,x(k)+2)
figure
meshc(XX,YY,y(x(k+1)*R))
end
大佬能不能看看，这段代码能运行吗，得到的应该是几张圆形区域的本征函数图
