matlab矩形法求定积分

矩形法是一种数值积分方法，用于近似计算定积分。在 MATLAB 中，可以使用以下代码实现矩形法求定积分：

function I = rectangle_integration(f, a, b, n)
% f: 被积函数
% a: 积分下限
% b: 积分上限
% n: 分割数

% 计算每个小矩形的宽度
dx = (b - a) / n;

% 初始化积分结果
I = 0;

% 计算每个小矩形的面积并将它们累加起来
for i = 1:n
    x = a + (i - 0.5) * dx;  % 计算小矩形中心的 x 坐标
    I = I + f(x) * dx;       % 计算小矩形面积并累加到积分结果中
end
end

% 示例：
% f(x) = x^2, 积分区间 [0, 1]
f = @(x) x.^2;
a = 0;
b = 1;
n = 100;
I = rectangle_integration(f, a, b, n);
fprintf('定积分的近似值为 %f\n', I);

在上述示例中，我们计算了函数 $f(x) = x^2$ 在区间 $[0, 1]$ 上的定积分，将区间分成了 100 个小矩形进行近似计算。您可以根据需要修改被积函数、积分区间和分割数进行计算。