function ab=sr_sqabel_type(dig,N)
%SR_SQABEL_TYPE Recurrence coefficients for squared 
%SR_SQABEL_TYPE Recurrence coefficients for squared 
%Abel-type polynomials.
%   AB=SR_SQABEL_TYPE(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)=2*cosh(pi*x/2)*[x/(2*sinh(pi*x/2))]^2 on 
%   [-Inf,Inf].

digits(dig);
ab=vpa(zeros(N,2));
ab(1,2)=vpa(4/pi);
for k=2:N
  if 2*floor((k-1)/2)~=k-1
    ab(k,2)=k^2;
  else
    ab(k,2)=k^2-1;
  end
end
