fhtx.m

来自「用于求积分」· M 代码 · 共 26 行

function [T1 T2 n h]=fhtx(f,a,b,n,eps)
h=(b-a)/n;
F=feval(f,a)+feval(f,b);
T1=((b-a)/2)*F;
F1=0;
for i=1:(n-1)
    x=a+i*h;
   F1=F1+feval(f,x);
   T2=h*F1+T1/n;
   c=T2-T1;
   if abs(c)<=eps
       i
       break;
       
   end
end 
T1=T2;return;
 
T2