function ab=r_assoc_legendre(N)
%R_ASSOC_LEGENDRE Recurrence coefficients for the associated
%Legendre weight function.
%   R_ASSOC_LEGENDRE(N) generates the first N recurrence coefficients
%   for the associated Legendre weight function

global mc mp iq idelta irout abj pi2
mc=2; mp=0; iq=1; idelta=2; irout=1; Mmax=500;
abj=r_jacobi01(500); pi2=pi^2;
eps0=.5e-6;
[ab,Mcap,kount]=mcdis(N,eps0,@quad_assoc_legendre,Mmax);
[Mcap,kount]
for k=1:9
  beta=(k+1)^2/((2*k+1)*(2*k+3))
end
