要用matlab對壹組離散數據（5000個）進行定積分求值，積分範圍是1000到2000，數據組成的圖像如下。

%復化cotes公式求積分值

%[a,b]為積分區間

%n是等分區間份數

load('C:\Users\baiyanglong\Desktop\n.mat')

a=1000;b=1500;

C=0;

h=4;

n=(b-a)/h;

for?i=1:(n-1)

x0=a+i*h;

C=C+14*nlp(x0);

end

for?k=0:(n-1)

x0=a+h*k;

s=32*nlp(x0+h*1/4)+12*nlp(x0+h*1/2)+32*nlp(x0+h*3/4);

C=C+s;

end

C=C+7*(nlp(a)+nlp(b));

C=C*h/90;

C=double(C);

采用復化cotes公式求積分值

結果積分結果為C=1.7581參考資料查數值分析教程。

其實可以簡略試算壹下根據梯形公式

sum(nlp(1000:1500))=1.7753

梯形公式誤差比較大