河北工业大学计算机软件技术基础(VC)课程设计报告
学院 信息工程学院 班级 姓名 _ 学号 _ 成绩 __ ____
一、题目:
设计采用梯形法和辛普生法求定积分(36)
二、设计思路
1、总体设计
1)分析程序的功能
本题目的功能是对梯形法和辛普生法,在不同区间数下计算所得的定积分的值,进行精度比较。
2)系统总体结构:设计程序的组成模块,简述各模块功能。
该程序共分为以下几个模块
模块一:各函数原型的声明。
模块二:主函数。
模块三:各函数的定义。
包括两个数学函数y1=1+x*x、y2=1+x+x*x+x*x*x的定义和两个函数指针
double integralt(double ,double ,int ,double(*f)(double))
double integrals(double ,double ,int ,double(*f)(double))
的定义。
2、各功能模块的设计:说明各功能模块的实现方法
模块一:对各种函数进行声明。
模块二:求梯形法和辛普森法,在不同区间数下计算所得的定积分的值。
模块三:将各函数写出来。
3、设计中的主要困难及解决方案
在这部分论述设计中遇到的主要困难及解决方案。
困难:函数指针的应用。解决方案:仔细阅读课本,以及与同学之间的讨论,和向老师求助。
4、你所设计的程序最终完成的功能
1)说明你编制的程序能完成的功能
在数学上求一个函数与x轴在一定范围内所围的面积即求定积分,对梯形法和辛普森法求定积分的比较。
2)准备的测试数据及运行结果
三、程序清单
#include
double integralt(double ,double ,int ,double(*f)(double));//用梯形法求定积分的声明.
double integrals(double ,double ,int ,double(*f)(double));//用辛普森法求定积分的声.明.
double f1(double x);//函数y1=1+x*x的声明.
double f2(double x);//函数y2=1+x+x*x+x*x*x的声明.
void main()
{
double a,b,intesum1,intesum2,intesum3,intesum4;//对求定积分的定义.
cout<
cin>>a>>b;
int n[7]={2,10,100,1000,5000,20000,50000};//n的不同取值.
//下面是对第一行不同n的值的输出.
cout<
for(int i=0;i<7;i++)
cout<
cout<
//下面是对n取不同值时,用梯形法对f1求定积分.
cout<
for(i=0;i<7;i++)
{
intesum1=integralt(a,b,n[i],f1);
cout<
}
cout<
//下面是对n取不同值时,用梯形法对f2求定积分.
cout<
for(i=0;i<7;i++)
{
intesum2=integralt(a,b,n[i],f2);
cout<
}
cout<
//下面是对n取不同值时,用用辛普森法对f1求定积分.
cout<
for(i=0;i<7;i++)
{
intesum3=integrals(a,b,n[i],f1);
cout<
}
cout<
//下面是对n取不同值时,用用辛普森法对f2求定积分.
cout<
for(i=0;i<7;i++)
{
intesum4=integrals(a,b,n[i],f2);
cout<
}
cout<
}
double f1(double x) //定义函数y1=1+x*x.
{
double y1;
y1=1+x*x;
return y1;
}
double f2(double x) //定义函数y2=1+x+x*x+x*x*x.
{
double y2;
y2=1+x+x*x+x*x*x;
return y2;
}
double integrals(double a,doubl