复化梯形公式与辛普森公式的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('正在退出，请稍候。。。');