%PLOT_ZEROS_HRGBINET

N=15
for k=1
%for k=N
%for k=1:N
  i=0; zp=zeros(10,1);
  for a=[5/100:5/100:1]
    i=i+1;
    z=zeros_hrgbinet(N,a);
    zp(i)=z(k);
  end  
  dmin=1; dmax=0;
  for r=2:20
    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:20)';
  plot(5*j/100,zp,'r');set(gca,'FontSize',14)
  axis([0 1 .06 .1])
%  axis([0 1 47.94 48.04])
%  axis([0 1 0 60])
  xlabel('alpha')
  ylabel('zeros')
  hold on
end
