数值分析实验复合辛普森公式报告.docVIP

实验题目：用复合辛普森公式求方程的积 学生姓名：何弯弯 专 业：信息与计算科学 学 号：0908060222 完成日期：2011/11/20 西安科技大学计算机科学与技术学院 实验题目：利用复合辛普森公式求方程的积 学生姓名： 何弯弯 学号： 0908060222 完成日期：2011/10/22 1 实验目的 1、了解复合辛普森公式的方法原理； 2、利用复合辛普森公式求方程的积 2.1 算法原理 原理：将区间[,b]等分成N个子区间[x(k),x(k+1)]（k=0,1,…..,N-1）， h=（b-a）/N，在每个子区间[x(k),x(k+1)]上使用公式，可求得结果。公式 其中是[]得中点，即 2.2 算法步骤 步骤一：输入a,b,h 步骤二：代入辛普森公式求的结果 2.3 程序流程图 3 实验结果分析 4实验心得体会 通过本次实验我熟悉用复化公式求数值积分的全过程，并更加熟悉的掌握了复化辛普生公式以及原理。使我对这门课有了进一步的理解，为以后的学习打下了良好的基础。#includeiostream.h #includemath.h double get(double x) { if(x==0) x=0.000000000001; return sqrt(x)*log(x); } double g(double a,double h,double b) { int i,n; double sum1,sum2=0,f; n=(int)((b-a)/h); sum1=get(a+h/2); for(i=1;i=n-1;i++) //求f(x(k))的和 { sum1+=get(a+i*h+h/2); sum2+=get(a+i*h); } f=4*sum1+2*sum2; return f; } void main() { double a=0,b=1,h;//给定的区间，以及步长 int n; while(true) { cout请输入要给的步长h(h0):; cinh; if(h==-1) break; double I1,f1,f2,f3; n=(int)((b-a)/h); f1=get(a); coutf1endl; f2=g(a,h,b); coutf2endl; f3=get(b); coutf3endl; I1=h/6*(f1+f2+f3);//复化辛普森公式 cout给定步长为h=hendl; cout复化辛普森公式计算结果:I=I1endl; cout所分成n个区间，n=nendl; } } 结束 输出I1，h,n i i=1,n,1 输入a,b,h 开始