vc0105.cpp

来自「由吕军等人著的《Visual C++ 与面向对象程序设计教程》课后习题的源代码」· C++ 代码 · 共 32 行

// Example 1-5：用梯形法计算定记分
#include <iostream.h>
#include <math.h>
// 定义被积函数
double f(double x)
{
	return sin(x);
}
// 主函数: 用梯形法计算定积分
void main()
{
	double a, b;		// 双精度类型变量: 积分的下限和上限
	double h;		// 双精度类型变量: 积分步长
	double sum;		// 双精度类型变量: 工作变量，最后为积分值
	int n;			// 整型变量 : 积分区间等分数
	int i;			// 整型变量 : 循环工作变量
	// 根据题意确定积分的下限、上限和积分区间等分数
	a = 0.0;
	b = 1.0;
	n = 1000;
	h = (b-a)/n;		// 计算小区间长度
	// 为工作变量赋初值 : 先计算无需循环运算的部分
	sum = (f(a)+f(b))/2;
	// 循环计算公式中的Σ和式
	for(i=1;i<n;i=i+1)
		sum = sum+f(a+i*h);
	// 完成计算，变量 sum 中存放积分结果
	sum = sum*h;
	// 输出计算结果
	cout<<"The result is " << sum << endl;
}