function z=zeros_hrgbinet(N,a)
%Zeros of the half-range generalized Binet polynomial. 
%   Z=ZEROS_HRGBINET(N,A) computes the N zeros in
%   ascending order of the half-range generalized Binet 
%   polynomial of degree N with parameter A, 0 < A <= 1.

z=zeros(N,1);
ab=sr_hrgbinet(16,16,N,a);
ab=double(ab);
xw=gauss(N,ab);
z=xw(:,1);
