function z=zeros_gbinet(dig,N,a)
%Zeros of the generalized Binet polynomial. 
%   Z=ZEROS_GBINET(N,A) computes the N zeros in
%   ascending order of the Binet polynomial of 
%   degree N with parameter A, 0 < A <= 1.

z=zeros(N,1);
ab=sr_gbinet(dig,6,N,a);
ab=double(ab);
xw=gauss(N,ab);
z=xw(:,1);
