辛普森积分法求圆周率matlab代码

下面是使用辛普森积分法计算圆周率的 Matlab 代码：

% 定义函数
f = @(x) 4 ./ (1 + x.^2);

% 辛普森积分法计算圆周率
a = 0;
b = 1;
n = 10000; % 积分区间等分数
h = (b - a) / n; % 步长
x = a:h:b; % 求出等分节点
y = f(x);
pi_approx = h / 3 * (y(1) + y(end) + 4*sum(y(2:2:end-1)) + 2*sum(y(3:2:end-2)));

% 输出结果
disp(pi_approx);

其中，`n` 表示积分区间等分数，可以根据需要进行调整。在上面的代码中，我们取 `n=10000` 进行计算，得到的圆周率近似值约为 3.141592653589793。