clear all; maxv = 10; maxs = 10; j= zeros(maxv, maxs); incr = 4.0; for v=0:maxv-1 h = v+1.9*v^(1/3)+1; if (v==0) j(v+1,1) = fzero(@(x)besselj(v,x),2); else j(v+1,1) = fzero(@(x)besselj(v,x),h); end for s=2:maxs j(v+1,s) = fzero(@(x)besselj(v,x),j(v+1,s-1)+incr); end end j