《实验二、数值积分》由会员分享,可在线阅读,更多相关《实验二、数值积分(2页珍藏版)》请在金锄头文库上搜索。
1、实验二、实验二、数值积分数值积分一、 【实验目的与要求】通过实验,观察复化梯形公式和复化辛普森公式随区间数增加各自误差n 减少的规律二、 【实验内容和步骤】计算定积分dxeIx201取,16, 8 , 4 , 2n步骤 1:用复化梯形公式步骤 2:用复化辛普森公式 步骤 3:与精确值进行比较。006994. 4I三、 【注意事项】1、求函数值需先运用建立函数句柄,再运用命令;feval例如:求函数在的函数值可用下面的命令:2)(xxf4xf=(x) x.2;y=feval(f,4)2、进行结果分析;参考程序:function s=agui_traps(fname,a,b,n)%fname 为被
2、积函数,a,b 分别为下上界,n 为等分数h=(b-a)/n;fa=feval(fname,a);fb=feval(fname,b);f=feavl(fname,a+h:h:b-h+0.001*h);t=h*(0.5(fa+fb)+sum(f);function s=agui_simpson(fname,a,b,n)%fname 为被积函数,a,b 分别为下上界,n 为等分数h=(b-a)/n;fa=feval(fname,a);fb=feval(fname,b);s=fa-fb;x=a;for i=1:nx=x+h/2;s=s+4*feval(fname,x);x=x+h/2;s=s+2*feval(fname,x);ends=s*h/6