trapezoidal.m

% 复合梯形公式
% 注释的地方需要修改
clear

a=0; b=pi/6;          % a为积分下界，b为积分上界
n=6;                  % 分割的区间数

syms x
f=(4-sin(x)^2)^0.5;   % 函数表达式

h=(b-a)/n;

for k=0:n              
    xk(k+1)=a+k*h;      
end

fa=subs(f,x,xk(1));
fb=subs(f,x,xk(n+1));

Tn=h/2*(fa+fb);

for k=2:n    
    fxk=subs(f,x,xk(k));
    Tn=Tn+h*fxk;
end

vpa(Tn,7)             % 结果有效数字位数