zerobes.m

global nu 

nu=0;

for n=1:30
zer(n)=fsolve(@bes,[n*pi],optimset('Display','off'));
end

zer