%PLOT_ZEROS_LSRBINET

N=15
for k=1
%for k=N
%for k=1:N
  i=0; zp=zeros(30,1);
  for c=[1/10:1/10:3]
    i=i+1;
    z=zeros_lsrbinet(N,c);
    zp(i)=z(k);
  end  
  dmin=1; dmax=0;
  for r=2:30
    d=(zp(r)-zp(r-1))/zp(r);
    if abs(d)<=dmin, dmin=abs(d); end
    if abs(d)>dmax, dmax=abs(d); end
    if d<=0
      disp('wrong');  
      [k r d] 
    end
  end
  [dmin dmax]
  j=(1:30)';
  plot(j/10,zp,'r');set(gca,'FontSize',14)
  axis([0 1 .06 .1])
%  axis([0 1 47.94 48.04])
%  axis([0 1 0 60])
  xlabel('c')
  ylabel('zeros')
  hold on
end
