function z=zeros_sqgbinet(N,a)
%Zeros of the squared generalized Binet polynomial. 
%   Z=ZEROS_SQGBINET(N,A) computes the N zeros in
%   ascending order of the polynomial of degree N
%   orthogonal with respect to the squared generalized
%   Binet weight function with parameter A, 0 < A < 1. 

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