ftrapz.m

来自「数值分析各种公式的解法」· M 代码 · 共 9 行

function t=ftrapz(a,b,n)
h=(b-a)/n;
t=h*(f(a)+f(b))/2;
for i=1:(n-1)
    t=t+h*f(a+h*i);
end
I=t
function f=f(x)
f=exp(sin(x));