//ODE_GillVariationalStep.cpp	
//變步長(zhǎng)積分一步基爾法

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

void main(void)
{	
	double q[3] = {0.0, 0.0, 0.0};
	double y[3] = {0.0, 1.0, 1.0};
	double t(0.0), h(0.1), eps(FLOATERROR);
	valarray<double> vy(y, 3), vq(q, 3);
			
  	cout.setf(ios::fixed);	//輸出數(shù)據(jù)為定點(diǎn)法
	cout.precision(6);		//精度6位

	cout << "t = " << t << endl;
	for(int i=0; i<3; i++)
		cout << "    y(" << i << ") = " << y[i] << "\t";
	cout << endl;

	for(i=1; i<11; i++)
	{
		ODE_GillVariationalStep(t,h,vy,eps,vq);
		t = h * i;
		cout << "t = " << t << endl;
		for(int j=0; j<3; j++)
			cout << "    y(" << j << ") = " << vy[j] << "\t";
		cout << endl;
    }
}

//計(jì)算微分方程組中各方程右端函數(shù)值
template <class _Ty>
void FunctionValueGVS(_Ty t, valarray<_Ty>& y, valarray<_Ty>& d)
{
	int n = y.size();	//微分方程組中方程的個(gè)數(shù)，也是未知函數(shù)的個(gè)數(shù)

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