//IntegralSurge.cpp		 高振蕩函數法求積

#include <iostream>		//輸入輸出流頭文件
#include "Integral.h"	//積分頭文件
using namespace std;	//名字空間

void main(void)
{	
	double eps = FLOATERROR, a(0), b(6.2831852);
    double fa[4] = {0, 1, 0, -3};
	double fb[4] = {6.2831852, 1, -6.2831852, -3};
	int m(30);
	valarray<double> ffa(fa, 4);
	valarray<double> ffb(fb, 4);
	valarray<double> s(2);
	
	cout << "IntegralSurge()" << endl << endl;

	IntegralSurge(a, b, m, ffa, ffb, s);

	cout << "s(0) = " << s[0] << "\t s(1) = " << s[1] << endl << endl;
}