Question

使用matlab编写好梯形公式的函数M文件后如何使用它去计算积分？

已知
复化梯形公式M文件：
function y=trap(‘f’,a,b,n)
h=(b-a)/n;
k=n-1;
sum=0;
for i=1:k
sum=sum+feval(f,a+i*h);
end
y=h/2*( feval (f,a)+2*sum+ feval (f,b));
请问这个应该怎么做？

Answer 1

被积函数，积分上下限当参数输进去，具体怎么做你要把m文件发上来
f是被积函数，a，b是积分上下限，n设的越大计算精度越高，f写成内联函数（inline），句柄函数都行