%PLOT_ZEROS_GBINET

N=30
%for k=ceil((N+2)/2)
%for k=N
for k=ceil((N+2)/2):N
  i=0; zp=zeros(10,1);
  for a=[5/100:5/100:1]
    i=i+1;
    z=zeros_gbinet(16,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 .5 .8])
%  axis([0 1 78.178 78.195])
 axis([0 1 0 100])
  xlabel('alpha')
  ylabel('zeros')
  hold on
end
