请教1阶贝塞尔函数绘图的问题
我自己完成的部分代码求出1阶贝塞尔函数的根，关于后面的部分不知如何进行下去为好，请指点一下clear all; maxv=2; maxs=10;
theta=linspace(0,2*pi,100); t=0.4;
=root_bes(maxv,maxs); 求bessel1阶函数的根
function =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 编辑 ]
页:
[1]