%PLOT_ZEROS_HRSQGBINET

N=15
for k=1
%for k=N
%for k=1:N
  i=0; zp=zeros(10,1);
  for a=[5/100:5/100:95/100]
    i=i+1;
    z=zeros_hrsqgbinet(N,a);
    zp(i)=z(k);
  end  
  dmin=1; dmax=0;
  for r=2:19
    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
      r
      [zp(r) zp(r-1)]
    end
  end
  %[dmin dmax]
  j=(1:19)';
  plot(5*j/100,zp,'r');set(gca,'FontSize',14)
  axis([0 1 .02 .05])
%  axis([0 1 23.88 24.04])
%  axis([0 1 0 30])
  xlabel('alpha')
  ylabel('zeros')
  hold on
end
