Matlab 梯形规则

我正在完成一个班级的作业。我们将按照流程图查找陷阱规则代码的值。我相信问题出在我的主代码上

我不确定我的功能代码或主代码是否有问题，如有任何帮助，将不胜感激

当我运行该部分时，它将函数显示为答案

以下是我的主脚本代码：

剩下的是我的梯形规则代码

返回f作为trap_haskell的输出参数，它是trap_haskell本身的输入函数。代码中的变量I实际上存储了积分，因此只需更改函数定义的输出变量即可返回积分：

%//       ------ Change here
%//       |
%//       V
function [I] = trap_haskell(f, a, b, n)

---

有什么问题？请阅读，非常感谢rayryeng，这解决了问题！

function [f] = trap_haskell(f, a, b, n)


x = a ;
h = (b - a) / n ;
s = f (a) ;

for k=1:1:n-1 
    x = x + h ;
    s = s + 2 * f(x) ;
end

s = s + f(b) ;
I = (b - a) * s / (2 * n) ;

end

%//       ------ Change here
%//       |
%//       V
function [I] = trap_haskell(f, a, b, n)