integral.cpp

一个比较小的代码

/************************************************************
  Copyright (C), 2007, 庞佩璐. Integral, SWUST.
  FileName: Integral.cpp
  Author:庞佩璐       Version :1.0          Date: 2007.1.10
  Description:积分计算器
  History:    历史修改记录
    <author>	<time>		<version >   
     庞佩璐		2007/1/10      1.0		 
***********************************************************/
#include<iostream>
using namespace std;
                      //对函数进行声明 
double Area(double );   
double X(double);
double F(double);

/*************************************************
  函数名称:      main
  函数功能:      求所给函数的积分
  Called:        none
  Called By:     主函数
  输入参数:          
		a:   积分函数的下限双精度数值
		b:   积分函数的上限双精度数值
		n:   要分的梯形的个数双精度数值
  输出:       
        sum:    要求的函数的积分值
  返回值:        无
  其它说明:      无  
*************************************************/
double main()
{
	extern double a,b,n,x,h,i;          //声明全局变量
	double Sum,k;
	cout << "----------欢迎您使用积分计算器----------" << endl;
	cout << "-                                      -" << endl;
	cout << "-  要求积分的函数为:" ;
	cout << "F(x) = 2 * X * X" << "   -" <<endl;
	cout << "-                                      -" << endl;
	cout << "----------------------------------------" << endl;
	for (k=0;k<=10;k++)
	{
	cout << "-----------------------------------------" << endl;
	cout << "请输入要求积分函数的下限:" ;
	                      //输入积分上下限的值和所分梯形的个数
	cin >> a ;
	cout << "请输入要求积分函数的上限:" ;
	cin >> b ;
	cout << "请输入要分梯形的个数:" ;
	cin >> n ;
	cout << "-----------------------------------------" << endl;

	                //判断a,b的大小
	if (a > b)
	{
		cout << "--------请正确输入积分上下限的值!--------" << endl;
	}
	else
	{
		for (i = 1; i <= n ;i ++)
		{
		   Sum = Area(i);   //积分求和式
		}
		cout << "---      该函数的积分值为:" ;
	    cout << Sum ;
		cout << "      ---" << endl;
	}
	}
	cout << "---------------谢谢您的使用--------------" << endl;
	return 0;
}
double a,b,n,x,h,i;   //定义全局变量

/*************************************************
  函数名称:      X
  函数功能:      提供求积分函数的变量
  Called:        none
  Called By:     调用函数
  输入参数:      无    
  输出:          无
  返回值:        x
  其它说明:      无  
*************************************************/
double X(double i)
{
	x = a + ( i - 1 ) * h;
	return x;
}

/*************************************************
  函数名称:      F
  函数功能:      提供所求积分函数
  Called:        none
  Called By:     调用函数
  输入参数:      无    
  输出:          无
  返回值:        d
  其它说明:      无  
*************************************************/
double F(double x)
{
	double d;
	d = 2 * x * x ;
	return d;
}

/*************************************************
  函数名称:      Area
  函数功能:      求体形面积的函数
  Called:        none
  Called By:     调用函数
  输入参数:      无    
  输出:          无
  返回值:        g
  其它说明:      无  
*************************************************/
double Area(double i)
{
	double g;
	h = (b - a) / n ;
	g = ((double)1/2) * (F(X(i)) + F(X(i+1)) * h);
	return g;
}