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

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