function ab=sr_sqabel(dig,N)
%SR_SQABEL Recurrence coefficients for square 
%Abel polynomials.
%   AB=SR_SQABEL(DIG,N) generates in DIG-digit
%   arithmetic the Nx2 array AB of the first N 
%   recurrence coefficients for polynomials orthogonal
%   with respect to the weight function
%   w(x)=[x/(2*sinh(pi*x))]^2 on [-Inf,Inf].

digits(dig);
ab=vpa(zeros(N,2));
k=2:N; c=[vpa(1/(12*pi)) vpa((k-1).*k.^2.*(k+1)./(4*(2*k-1).*(2*k+1)))]; 
ab(:,2)=c';
