《复化Simpson公式.docx》由会员分享,可在线阅读,更多相关《复化Simpson公式.docx(4页珍藏版)》请在第壹文秘上搜索。
1、2023-2023(1)专业课程实践论文复化Simpson公式肖云龙,0818180214,R数学08-2班一、算法理论SimPSon公式是一个数值积分公式,在计算一些多项式函数(三次或三次一下)的定积分时会得出精确值。但容易验证它对于/(X)=/通常是不准确的,因此,SimPSon公式实际上具有三次代数精度.提高阶的途径并不总能取得满意的效果。为了改善求积的精度,通常采用复化求积法。Simpson公式将定积分=/x)dx的积分区间a,b分隔为等分,各节点为,c,buxi-a+jh,J=OJ,,*6=复化梯形公式的形式记子区间/的中点为X1那么复化Simpson公式:二、算法框图include
2、开始三、算法程序includedoublefunction(doublex)doubles;s=x(4+x*x);returns;doubleReiterationOfSimpson(doublea,doubleb,doublen,doublef(doublex)(doubleh,fa,fb,xk,xj;h=(b-a);fa=f(a);fb=f(b);doubles1=0.0;doubles2=0.0;for(intk=l;kn;k+)xk=a+k*h;sl=sl+f(xk);for(intj=O;jn;j+)xj=a+(j+0.5)*h;s2=s2+f(xj);doublesn;sn=h6*
3、(fa+fb+2*s1+4*s2);returnsn;)main()doublea,b,Result,n;coutM请输入积分下限:a;coutM请输入积分上限:b;COUt”请输入分割区间数n:*n;CoUt复化Simpson公式计算结果:;Result=ReiterationOfSimpson(a,b,n,function);coutResultendl;return0;四、算法实现例1.用Simpson公式求积分/以X?解:运行程序(1)请输入积分下限:输入口的值为0,回车。(2)请输入积分上限:输入方的值为1,回车。请输入分割区间数:输入的值为10,回车。(4)显示结果:复化SimPSon公式计算结果:0.111572例2.用Simpson公式求积分:解:运行程序(1)请输入积分下限:输入a的值为0,回车。(2)请输入积分上限:输入人的值为1,回车。(3)请输入分割区间数n:输入的值为10,回车。(4)显示结果:复化SimPSon公式计算结果:0.333333