integralsimpson1d.cpp

计算机常用数值计算算法与程序 C++版 各章代码

//IntegralSimpson1D.cpp		 变步长辛卜生法求积

#include <iostream>			//输入输出流头文件
#include "Integral.h"		//积分头文件
using namespace std;		//名字空间

void main(void)
{	
	double eps = FLOATERROR, a(0), b(1);
                
	cout << "IntegralSimpson1D()" << endl << endl;

	double t = IntegralSimpson1D(a, b, eps);

	cout << "t = " << t << endl << endl;
}

// 计算被积函数值
template <class _Ty>
_Ty FunctionValueIS1D(_Ty x)
{
	_Ty y;

	y = log(1.0 + x) / (1.0 + x * x);
    
	return(y);
}