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

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