poly_itg.m

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

% poly_itg(p) integrates a polynomial p which is 
% in a power series. The result is also a power series.
% Copyright S. Nakamura, 1995
function py = poly_itg(p)%求多项式积分后的系数
n=length(p);
py = [p.*[n:-1:1].^(-1),0];