function [a,b] = coeflege(n)

% Coefficients of a=\alpha and b=\beta of the Legendre polynomials
 
if( n <=1), disp( ' n must be > 1'); return;  end

for k=1:n, a(k) = 0;  b(k) = 0; end;  b(1)=2;

for k=2:n, b(k) = 1/(4-1/(k-1)^2); end