复化科特斯公式matlab_matlab实现复化Newton-Cotes公式求积分的程序应用和代码

执行函数为mymulNewtonCotes.m

1、使用方法：

Step1：在MATLAB 命令窗口输入被积函数212

0t t e dt ?。 输入应为：ft=@(t)t.*exp(t^2/2)。

Step2：执行函数。输入形式为mymulNewtonCotes(ft,a,b,m,n)；

其中ft —被积函数，此体重ft=@(t)t.*exp(t^2/2)，已经在第一步赋值；

a —积分下限，本题中为0；

b —积分上限，本题中为1；

m —将区间[a,b]等分的子区间数量，本题可选为10；

n —采用的Newton-Cotes 公式的阶数，必须满足n<8，否则积分没法保

证稳定性。

当n=1时，即为复化梯形公式；n=2时，即为复化复化辛普森公式。

所以，分别输入mymulNewtonCotes(ft,0,1,10,1)和

mymulNewtonCotes(ft,0,1,10,2)就可以得到两种方法的积分计算结果。

2、计算结果

而根据积分运算，可得：

221112

110222

220000() 1.648710.64872t t x x t t e dt e d e dx e e e ====-=-=??? 说明复化梯形和复化辛普森公式计算出的结果基本一致，与实际结果相符。