[julia-users] Re: fzero

[digxx]

Almost forgot:
r is defined by:
g_r(a)=besselj1(a)-a/2*besselj0(a);

N=500;
rr=zeros(1,1,N);
r=zeros(1,1,N);
index=zeros(5,N);

for i=1:N
r[1,1,i]=fzero(g_r,[pi+(2*i-1),2*pi+(2*i-1)]);
end

eps=0.1;
for i=1:N
if ~(i in index)
ind=find( abs(r-r[1,1,i]) .< eps );
rr[1,1,i]=mean(r[ind]);
index[1:length(ind),i]=ind;
end
end

ind=find(rr .!= 0)
r=zeros(1,1,length(ind))
r[1,1,:]=rr[ind]