Newton-Cotes公式

科茨系数为

C

?n?kn?1???t?t?1???t?k?1??t?k?1???t?n?dt n?k!?n?k!?0n?k

Matlab代码

function NewtonCotes

for n=1:6

for k=0:n

Cotes(n,k+1)=CotesCoeffsCal(n,k);

end;

end;

return;

function Cotes=CotesCoeffsCal(n,k)

syms t;

a=t*ones(1,n+1);

b=0:1:n;

a=a-b;

a(k+1)=1;

Cotes=int(prod(a),t,0,n);

Cotes=(-1)^(n-k)/n/prod(1:k)/prod(1:n-k)*Cotes;

return;

Cotes =

[ 1/2, 1/2, 0, 0, 0, 0, 0]

[ 1/6, 2/3, 1/6, 0, 0, 0, 0]

[ 1/8, 3/8, 3/8, 1/8, 0, 0, 0]

[ 7/90, 16/45, 2/15, 16/45, 7/90, 0, 0]

[ 19/288, 25/96, 25/144, 25/144, 25/96, 19/288, 0]

[ 41/840, 9/35, 9/280, 34/105, 9/280, 9/35, 41/840]

Newton-Cotes公式.doc下载

久久建筑网m.kkreddy.com提供大量：建筑图纸、施工方案、工程书籍、建筑论文、合同表格、标准规范、CAD图纸等内容。