matlab复合辛普森求定积分

Matlab中可以使用`quad`函数进行定积分的求解，也可以使用`quadl`函数进行复合辛普森求解。

`quadl`函数的语法为：

Q = quadl(fun,a,b)
Q = quadl(fun,a,b,tol)
Q = quadl(fun,a,b,tol,trace)
Q = quadl(fun,a,b,tol,trace,P1,P2,...)

其中，`fun`为被积函数的句柄，`a`为积分下限，`b`为积分上限，`tol`为误差容限，`trace`表示是否显示迭代过程，`P1,P2,...`为被积函数的参数。

以下是一个使用`quadl`函数求解复合辛普森求积的示例代码：

fun = @(x) exp(-x.^2/2)/sqrt(2*pi); % 定义被积函数
a = -Inf; % 积分下限
b = 1.96; % 积分上限
tol = 1e-6; % 误差容限
Q = quadl(fun,a,b,tol); % 使用quadl函数求解定积分

其中，被积函数为高斯分布的概率密度函数，积分下限为负无穷，积分上限为1.96，误差容限为1e-6，函数返回的结果`Q`即为积分值。