计算方法--辛普森算法实验报告.docVIP

PAGE  PAGE 5 昆 明 理 工 大 学 计 算 方 法 学院：信息工程与自动化 专业：计算机科学与技术 班级：2006级计科（三）班 姓名：凌意平 学号：200830405906 昆明理工大学信息工程与自动化学院学生实验报告 课程名称：计算方法 开课实验室：2301 2008年12月4日 年级、专业、班计算机科学与技术 063班学号200830405906姓名凌意平成绩实验项目名称用复合辛普生公式计算指导教师郭 洪教师评语 教师签名： 2008年 12 月 日 实验目的 用复合辛普生公式计算。 实验原理 由插值型求积公式，，当n为不同的数值时，所表示的公式就不一样，所求的值也不一样。当n＝1时，为梯形求积公式，当n＝2时，有三点式，称为辛普生公式，公式如下： ，当n＝4时，有5点式，称为牛顿柯特斯公式。 程序流程图 四、所用仪器、材料（设备名称、型号、规格等） 所用仪器：兼容计算机一台 所用软件：Visual c++ 五、程序代码 #include stdio.h #include math.h//头文件 float f(float x)//定义函数 {return sqrt(x);}//利用主函数传来的实参，之间返回其平方根 void main()//主函数 {float a,b,t;//定义三个浮点型变量 l1:printf(请输入积分的区间下限和上限: a,b\n); //设置标号l1，并让用户输入积分的下限和上限 scanf(%f,%f,a,b);//把用户输入的上下限分别赋给b和a if (ab) //如果a大于b {printf(输入积分的积分区间不对，请重新输入！\n\n); //提示用户重新输入，因为积分的下限大于上限了。 goto l1;}//并让程序跳转到l1处 else t=(b-a)/6*(f(a)+4*f((a+b)/2)+f(b));//利用辛普生公式计算结果 printf(最后结果为:%f\n\n\n,t);//输出最后结果 } 六、程序调试 在运行环境中调试成功后能够得到正确的运行结果，输入积分的上下限，若下限大于上限时，系统将提示错误，并要求用户重新输入，按照题目的要求，我们应输入“2，3”，系统将得到正确的结果。如下图所示。 七、实验总结及教训 这个程序整体来说，非常的简单，我觉得没费很大的力气就求了出来，比前面两个程序都要简单。本来想把程序更加的精简的，本想不用调用函数，直接在程序首部定义函数，如下：#define f(x) sqrt(x),这样程序更加的简洁了，但是没有办法通过编译,老是提示有错误，也改不出来，后来还是决定麻烦一点，定义一个函数的调用。因为害怕积分的上下限被用户输反，所以增加了一个条件判断。 通过次实验的编写，加深了我对复合辛普生公式计算积分的理解，并加深了对公式的记忆。