自己编写一个变步长积分公式（MATLAB）

在MATLAB中编写一个简单的变步长积分函数，你需要模拟矩形法则的过程，但考虑到实际的非均匀步长，我们需要手动计算每个子区间内的积分。这里有一个基础版本的例子，假设我们已经有了一个一维的输入向量（x_values），以及对应值的函数（y_values），我们将对它们进行累加：

function result = my_variable_step_integral(x_values, y_values)
    % 初始化总和为0
    total_integral = 0;
    
    % 遍历x值
    for i = 1 : length(x_values) - 1
        % 计算当前子区间的长度（步长）
        dx = x_values(i+1) - x_values(i);
        
        % 计算并累加该区间的面积（y值乘以dx）
        total_integral = total_integral + y_values(i) * dx;
    end
    
    % 别忘了处理最后一个元素，它可能不是完整区间
    if ~isempty(x_values(end))
        total_integral = total_integral + y_values(end-1) * (x_values(end) - x_values(end-1));
    end

    result = total_integral;
end

在这个函数中，`x_values` 是x轴的非均匀值，`y_values` 是对应的函数值。注意这只是一个基本示例，如果函数复杂或者需要更精确的算法，你可能需要使用数值积分库如`trapz` 或者更高级的 `integral` 函数。