MATLAB 复合梯形公式和复合辛普森公式代码逐行解释

clear;：清除命令窗口中的所有变量。
clc;：清空命令窗口中的所有内容。
h=0.001;：将步长h设为0.001。
n=1/h;：计算数据点的数量，即n=1000。
t=0;：将t的初始值设为0。
for i=1:n-1：循环i从1到n-1。
t=t+fun41(i*h);：将t加上fun41(i*h)的值。
end：结束循环。
T=h/2*(0+2*t+fun41(1));：使用复合梯形公式计算积分近似值T。
T=vpa(T,10)：将T精确到小数点后10位。
s1=0;：将s1的初始值设为0。
s2=0;：将s2的初始值设为0。
for i=0:n-1：循环i从0到n-1。
s1=s1+fun41(h/2+i*h);：将s1加上fun41(h/2+i*h)的值。
end：结束循环。
for i=1:n-1：循环i从1到n-1。
s2=s2+fun41(i*h);：将s2加上fun41(i*h)的值。
end：结束循环。
S=h/6*(0+4*s1+2*s2+fun41(1));：使用复合辛普森公式计算积分近似值S。
S=vpa(S,10)：将S精确到小数点后10位。