复化梯形公式与辛普森公式的matlab程序
【程序代码】
cclc;
disp('1.复化梯形公式求解');
disp('2.simpson公式求解');
disp('请进行选择:');
c=input('');
ifc==1
clc;
disp(吃春房搞生随类味巴减乱'复化梯形公式');
disp('请输入积分下限');
a=input('a=');
disp('请输入积分上限');
b=input('b=');
disp('请输入等分的数目');
n=input('n=');
h=(b-a)/n;
s1=0;
fori=1:n-1
s1=s1+fun1(i*h);
end
disp('复化梯形公式的结果:');
T=h表乱第话汉即/2*(fun1(a)+2*s1+fun1(b))
elseifc==2
clc;
disp('simpson公式');
disp('医请输入积分下限');
a=input('a=');
disp('请输入积分上限');
b=input('b=');
disp('请输入等分的数目');
n=input('n=');
h=(b-a)/n;
s2=0;
fori=0:n-1
s2=s2+fun1((i+0.5)*h);
end
唱省行工三传曲算光disp('辛普森公式的结果:');
S=h/6很善开个*(fun1(a)+4*s2础举续作牛原抗衣切内+2*s1+fun1(b))
end
end
disp('菜单选项');
disp('1.继续运算');
disp('2.退出程序!');
p=input('');
i察fp==1
(fuhua);
早备福elseifp==2
disp('正在退出,请稍候。。。');
(exit);
end
end