% SSOBZEROS Symbolic analogue of sobzeros.m
%
function z=ssobzeros(dig,n,N,B)
if(n<1|n>N), error('n out of range'), end
syms H
digits(dig)
B=vpa(B);
H=vpa(zeros(n));
for i=1:n
  for j=1:n
    if i==1
      H(i,j)=B(j,j);
    elseif j==i-1
      H(i,j)=vpa(1);
    elseif j>=i
      H(i,j)=B(j-i+1,j);
    end
  end
end
z=eig(H);
[zs,I]=sort(subs(z));
z=z(I);
