tixing.m

为配套王能超版本数值分析简明教程编写,方便学生教师演示!

function s=tixing(f1,a,b,n)
f=inline(f1);
h=(b-a)/n;
s=0;
s=h/2*(f(a)+f(b));
for i=1:n-1
    s=s+h*f(a+i*h);
end
fprintf('%13.8f',s);