%PLOT_ZEROS_LSSRBINET

global dig
N=30
%for k=1
%for k=N
%for k=1:N
for k=ceil((N+2)/2):N
  k
  i=0; zp=zeros(30,1);
  for c=[1/10:1/10:3]
    if c<1/2
      dig=160;
    elseif c<1
      dig-124;
    elseif c<3/2
      dig=108;
    elseif c<2
      dig=96;
    elseif c<5/2
      dig=92;
    elseif c<3
      dig=84;
    end
    i=i+1;
    z=zeros_lssrbinet(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',24)
%  axis([0 1 .06 .1])
%  axis([0 1 47.94 48.04])
%  axis([0 1 0 60])
  xlabel('c')
  ylabel('zeros')
  hold on
end
