计算方法实验报告(5)
学生姓名杨贤邦学号指导教师吴明芬实验时间2014.4.16地点综合实验大楼203
实验题目数值积分方法
实验目的●利用复化梯形、辛普森公式和龙贝格数值积分公式计算定积分的
近似植。
实验内容●梯形、辛普森、柯特斯法及其Matlab实现;
●变步长的梯形、辛普森、柯特斯法及其Matlab实现。
●题目由同学从学习材料中任意选两题
算法分析梯形:function y=jifeng_tixing(a,b,n,fun)
fa=feval(fun,a);
fb=feval(fun,b);
s=0;
h=(b-a)/n;
for k=1:n-1
xk=a+k*h;
s=feval(fun,xk)+s;
end
y=(h/2)*(fa+fb+2*s);
辛普生:function y=jifeng_xingpu(a,b,n,fun) fa=feval(fun,a);
fb=feval(fun,b);
h=(b-a)/n;
s=0;
s2=feval(fun,a+0.5*h);
for k=1:n-1
xk=a+k*h;
s=feval(fun,xk)+s;
s2=feval(fun,xk+(h/2))+s2;
end