legen_pw(n)

% function legen_pw(n) finds the coefficients of a 
% Legendre polynomial of order n. 
% Copyright S. Nakamura, 1995
function pn = legen_pw(n)
pbb=[1];  if n==0, pn=pbb; break; end
pb=[1 0]; if n==1, pn=pb;  break; end
for i=2:n;
   pn=  ( (2*i-1)*[pb,0] - (i-1)*[0, 0, pbb] )/i;
   pbb=pb; pb=pn;
end