RES =辛普森(Y) Y通过计算积分的近似值
辛普森的1/3规则(单元间距)。辛普森的1/3规则使用
二次interpolants数值积分。计算
积分间距不同,RES乘以间距
增量。

为向量,辛普森(Y)的积分Y矩阵,辛普森(Y)
是一个行向量积分每一列。为一天
数组,辛普森(Y)跨越第一个单体尺寸。

RES =辛普森(X, Y)计算Y对X使用的积分
辛普森的1/3。X和Y都必须相同的向量
必须一个列向量长度,或者X和Y的第一个数组
单体尺寸长度(X)。辛普森是沿着这
维度。注意,必须正确执行的等距的X
1/3和3/8规则。如果X是不等距的,梯形法则
建议(MATLAB TRAPZ)。

RES =辛普森(X, Y,昏暗的)或者辛普森(Y,昏暗)集成了跨维度
昏暗的Y . X的长度必须相同,大小(Y,昏暗))。

RES =辛普森(X, Y,昏暗,规则)可以用来辛普森之间切换的1/3
规则和辛普森的3/8规则。辛普森的3/8规则使用立方interpolants
完成的数值积分。如果昏暗的默认值
需要指定一个空矩阵。

——规则选项

(默认)的1/3为二次interpolants辛普森法则

“3/8”为立方interpolants辛普森法则

例子:
%集成Y = SIN (X)
x = 0:0.2:π;
y = sin (x);
一笔= (y) * 0.2;%矩形规则
b = trapz (x, y);%梯形法则
c =辛普森(x, y, [], 1/3);辛普森%的1/3
d =辛普森(x, y, [],“3/8”);%辛普森的3/8规则
e = cos (x (1)) cos (x(结束));%实际积分
流(“矩形规则:%。15 f \ n”)
流(“梯形法则:%。15 f \ n”, b)
流(“辛普森”年代1/3规则:%。15 f \ n”, c)
流(“辛普森的规则:3/8 %。15 f \ n”, d)
流(“实际积分:%。15 f \ n”, e)

%http://math.fullerton.edu/mathews/n2003/simpson38rule/Simpson38RuleMod/Links/Simpson38RuleMod_lnk_2.html
x1 = linspace (0、2、4);
x2 = linspace (0、2、7);
x4 = linspace (0、2、13);
y = @ (x) 2 + cos(2 *倍根号(x));
格式长
日元= y (x1);res1 =辛普森(x1, y1, [],“3/8”);disp (res1)
y2 = y (x2);它=辛普森(x2, y2, [],“3/8”);disp(它)
y4 = y (x4);res4 =辛普森(x4、y4 [],“3/8”);disp (res4)

类支持输入X万博1manbetx, Y:
浮:双、单

参见总和,cumsum trapz cumtrapz。