clear;clc;m=16;h=1;T1=0+fun411h2for i=2m h=h2; n=1h; t=0; for j=12n-1 t=t+fun41jh; end Ti=Ti-12+ht;梯形公式endfor i=1m-1 for j=mi+1 Tj=4^i4^i-1Tj-14^i-1Tj-1; 通过不断

clear; %清空命令窗口中的所有变量 clc; %清空命令窗口中的所有输出 m=16; %设置迭代次数为16次 h=1; %初始步长为1 T(1)=(0+fun41(1))h/2 %使用梯形公式计算第一次迭代的结果，其中fun41为自定义的函数 for i=2:m %进行16次迭代 h=h/2; %每次迭代步长缩小一半 n=1/h; %计算每次迭代需要计算的点的个数 t=0; %初始化t for j=1:2:n-1 %计算每次迭代中的奇数点 t=t+fun41(jh); %累加每个奇数点的函数值 end T(i)=T(i-1)/2+h*t;%使用梯形公式计算当前迭代的结果 end for i=1:m-1 %进行Richardson外推法迭代 for j=m:i+1 T(j)=4^i/(4^i-1)*T(j)-1/(4^i-1)*T(j-1); %通过外推公式计算T表的对角线上的元素 end end vpa(T(m),10) %将计算结果保留10位小数并输出