请教1阶贝塞尔函数绘图的问题

wusemm

我自己完成的部分代码求出1阶贝塞尔函数的根，关于后面的部分不知如何进行下去为好，请指点一下
clear all; maxv=2; maxs=10;
theta=linspace(0,2*pi,100); t=0.4;
[r]=root_bes(maxv,maxs); 求bessel1阶函数的根


function [r]=root_bes(maxv,maxs)
r= zeros(maxv-1, maxs); incr = 4.0;
for v=maxv-1
    h = v+1.9*v^(1/3)+1;
    if (v==0)
        r(v+1,1) = fzero(@(x)besselj(v,x),2);
    else
        r(v+1,1) = fzero(@(x)besselj(v,x),h);
    end
    for s=2:maxs
        r(v+1,s) = fzero(@(x)besselj(v,x),r(v+1,s-1)+incr);
    end   
end
r

[ 本帖最后由 ChaChing 于 2010-3-30 21:25 编辑 ]