legen_pw.m

利用高斯-勒让德多项式计算 sin(x+y)在矩形区域的积分

% 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;
    pbb
    
   % break; 
end
pb=[1 0]; 
if n==1, %1次勒让德多项式系数
   pn=pb; 
   % break; 
end
for i=2:n;
   pn=  ( (2*i-1)*[pb,0] - (i-1)*[0, 0, pbb] )/i;;%利用递推三项式求出n次勒让德多项式系数
   pbb=pb;
   pb=pn;
end