function ab=sr_associated_legendre(dig,N)
%SR_ASSOCIATED_LEGENDRE Recurrence coefficients for associated
%Legendre polynomials.
%   AB=SR_ASSOCIATED_LEGENDRE(DIG,N) generates in DIG-digit
%   arithmetic the Nx2 array AB of the first N recurrence 
%   coefficients for associated Legendre polynomials 

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