// Integral.h		數(shù)值積分頭文件
// Ver 1.0.0.0
// 版權(quán)所有(C) 2002
// 版權(quán)所有(C) 何渝, 2002
// 最后修改: 2002.5.31.

#ifndef _INTEGRAL_H		//避免多次編譯
#define _INTEGRAL_H

#include <comm.h>		//公共頭文件
#include <random.h>		//隨機(jī)算法頭文件
#include <valarray>		//模板類valarray的標(biāo)準(zhǔn)頭文件

using namespace std;	//名字空間


//變步長梯形法求積
template <class _Ty>
_Ty IntegralTrapezia(_Ty a, _Ty b, _Ty eps);

//變步長辛卜生法求積
template <class _Ty>
_Ty IntegralSimpson1D(_Ty a, _Ty b, _Ty eps);

//自適應(yīng)梯形法求積
template <class _Ty>
_Ty IntegralTrapeziaSelfAdapt(_Ty a, _Ty b, _Ty eps, _Ty d);

//龍貝格法求積
template <class _Ty>
_Ty IntegralRomberg(_Ty a, _Ty b, _Ty eps);

//一維連分式法求積
template <class _Ty>
_Ty IntegralFraction1D(_Ty a, _Ty b, _Ty eps);

//高振蕩函數(shù)法求積
template <class _Ty>
void IntegralSurge(_Ty a, _Ty b, int m, valarray<_Ty>& fa, 
								valarray<_Ty>& fb, valarray<_Ty>& s);
//勒讓德-高斯法求積
template <class _Ty>
_Ty IntegralLegendreGauss(_Ty a, _Ty b, _Ty eps);

//拉蓋爾-高斯法求積
template <class _Ty>
void IntegralLaguerreGauss(_Ty& dValue);

//埃爾米特-高斯法求積
template <class _Ty>
void IntegralHermiteGauss(_Ty& dValue);

//切比雪夫法求積
template <class _Ty>
_Ty IntegralChebyshev(_Ty a, _Ty b, _Ty eps);

//蒙特卡洛法求積
template <class _Ty >
_Ty IntegralMonteCarlo1D(_Ty a, _Ty b);

//二重變步長辛卜生法求積
template <class _Ty>
_Ty IntegralSimpson2D(_Ty a, _Ty b, _Ty eps);

//多重高斯法求積
template <class _Ty >
_Ty IntegralGaussMD(valarray<_Ty>& js);

//二重連分式法求積
template <class _Ty >
_Ty IntegralFraction2D(_Ty a, _Ty b, _Ty eps);

//多重蒙特卡洛法求積
template <class _Ty >
_Ty IntegralMonteCarlo2D(valarray<_Ty>& a, valarray<_Ty>& b);

#include "Integral.inl"		//類及相關(guān)函數(shù)的定義頭文件

#endif		// _INTEGRAL_H