ftrapz.asv

数值分析各种公式的解法

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=