复化梯形公式与辛普森公式的matlab程序
【程序代码】
cclc;
disp('1.复化梯形公式求解');
disp('2.simpson公式求解');
disp('请进行选择:');
c=input(' ');
if c==1
clc;
disp('复化梯形公式');
disp('请输入积分下限 ');
a=input('a=');
disp('请输入积分上限 ');
b=input('b=');
disp('请输入等分的数目 ');
n=input('n=');
h=(b-a)/n;
s1=0;
for i=1:n-1
s1=s1+fun1(i*h);
end
disp('复化梯形公式的结果:');
T=h/2*(fun1(a)+2*s1+fun1(b))
else if c==2
clc;
disp('simpson公式');
disp('请输入积分下限 ');
a=input('a=');
disp('请输入积分上限 ');
b=input('b=');
disp('请输入等分的数目 ');
n=input('n=');
h=(b-a)/n;
s2=0;
for i=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(' ');
if p==1
(fuhua);
else if p==2
disp('正在退出,请稍候。。。');