辛普森公式matlab代码

以下是辛普森公式的 MATLAB 代码：

function I = simpson(f, a, b, n)
% f: 被积函数
% a: 积分下限
% b: 积分上限
% n: 区间数，必须为偶数
h = (b - a) / n;
x = a:h:b;
y = f(x);
I = h / 3 * (y(1) + 4 * sum(y(2:2:end-1)) + 2 * sum(y(3:2:end-2)) + y(end));
end

希望对您有帮助！

相关问题

复化辛普森公式MATLAB代码

复化辛普森公式是一种数值积分方法，用于计算函数在给定区间上的定积分。下面是使用MATLAB实现复化辛普森公式的代码：

function result = composite_simpson(f, a, b, n)
    % f: 被积函数
    % a: 积分下限
    % b: 积分上限
    % n: 划分的子区间数
    
    h = (b - a) / n; % 子区间宽度
    x = a:h:b; % 子区间节点
    
    % 计算每个子区间的积分值
    y = f(x);
    
    % 计算复化辛普森公式的近似积分值
    result = (h/3) * (y(1) + 4*sum(y(2:2:end-1)) + 2*sum(y(3:2:end-2)) + y(end));
end

使用方法：
1. 定义被积函数f(x)，例如：`f = @(x) x.^2;`
2. 调用`composite_simpson`函数计算积分值，例如：`result = composite_simpson(f, a, b, n);`

其中，`f`为被积函数，`a`和`b`为积分区间的下限和上限，`n`为划分的子区间数。返回值`result`为复化辛普森公式的近似积分值。

辛普森公式matlab

辛普森公式是一种数值积分方法，可以用于求解定积分。在MATLAB中，可以使用`integral`函数来实现辛普森公式。例如，要求解函数`f(x) = x^2`在区间`[1, 2]`上的定积分，可以使用以下代码：

f = @(x) x.^2;
integral(f, 1, 2, 'Method', 'Simpson')

其中，`f`为要积分的函数句柄，`1`和`2`分别为积分区间的下限和上限，`'Method', 'Simpson'`表示使用辛普森公式进行积分。运行代码后，MATLAB会输出积分的近似值。