数值分析matlab编程实现复化梯形公式

n T_n = C_t_r(f,n,a,b)
%f为函数,n为精度,a,b为区间端点
h=(b-a)/n;
%f=sqrt(x);
for i=0:n
    x(i+1)=a+i*h;
end
T_1=h/2*(1+f(b));
T_2=0;
for j=2:n
    F(j)=h*f(x(j));
    T_2=T_2+F(j);
end
%T_2=sum(F);
T_n=T_1+T_2;`