function ab=sr_lindeloef_dahlquist(dig,N)
%SR_LINDELOEF Recurrence coefficients for Lindeloef-Dahlquist 
%   polynomials.
%   AB=SR_LINDELOEF_DAHLQUIST(DIG,N) generates in DIG-digit
%   arithmetic the Nx2 array AB of the first N recurrence 
%   coefficients for Lindeloef-Dahlquist  polynomials orthogonal 
%   on [-Inf,Inf] with respect to the weight function 
%   w(x)=1/(2*cosh(pi*x)).

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