//ODE_BothSidesInterval.cpp	
//全區間積分連分式法

#include <iostream>							//輸入輸出流頭文件
#include "OrdinaryDifferentialEguation.h"	//求解常微分方程(組)頭文件
using namespace std;						//名字空間

void main(void)
{	
	double y[2] = {1.0, 0.0}, t(0.0), h(0.1), eps(FLOATERROR);
	valarray<double> vy(y, 2);
	matrix<double> vz(2,11);
	
  	cout.setf(ios::fixed);	//輸出數據為定點法
	cout.precision(6);		//精度6位

	ODE_BothSidesInterval(t,h,vy,eps,11,vz);

	for(int i=0; i<11; i++)
	{
		double x = h * i;
		cout << "t = " << x << endl;
		for(int j=0; j<2; j++)
			cout << "    y(" << j << ") = " << vz(j, i) << "\t";
		cout << endl;
    }
}

//計算微分方程組中各方程右端函數值
template <class _Ty>
void FunctionValueBSI(_Ty t, valarray<_Ty>& y, valarray<_Ty>& d)
{

	int n = y.size();	//微分方程組中方程的個數，也是未知函數的個數

    d[0] = -y[1];
	d[1] = y[0]; 
}