clear;clc;h=0001; h为步长可分别令h=1010010001n=1h; t=0; for i=1n-1 t=t+fun41ih;endT=h20+2t+fun411;T=vpaT10 以上为复合梯形公式 以下为复合辛普森公式s1=0;s2=0;for i=0n-1 s1=s1+fun41h2+ih;endfor i=1n-1 s2=s2+fun41ih;

clear; %清空工作区变量 clc; %清空命令窗口

h=0.001; %设定步长为0.001 n=1/h; %计算步数 t=0; %初始化t为0

% 循环计算f(x)的值并累加 for i=1:n-1 t=t+fun41(i*h); end

% 计算复合梯形公式的近似值 T=h/2*(0+2*t+fun41(1)); T=vpa(T,10) %将结果保留10位有效数字输出

% 初始化s1和s2为0 s1=0; s2=0;

% 循环计算f(x)的值并累加 for i=0:n-1 s1=s1+fun41(h/2+i*h); end

for i=1:n-1 s2=s2+fun41(i*h); end

% 计算复合辛普森公式的近似值 S=h/6*(0+4s1+2s2+fun41(1)); S=vpa(S,10) %将结果保留10位有效数字输出