数值分析matlab编程实现复化辛普森公式

%复化辛普森公式
function 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/6*(f(a)+f(b));
T_2=0;
T_3=0;
for j=2:n
    F(j)=h/3*f(x(j));
    T_2=T_2+F(j);
end
for k=n+2:2*n+1
    F(k)=h/3*(f(x(k))+h/2);
    T_3=T_3+F(k);
end
%T_2=sum(F);
T_n=T_1+T_2+T_3;