//ODE_LinearBoundaryValude.cpp	
//定步長全區(qū)間積分的歐拉法

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

void main(void)
{	
	double a(2.0), b(3.0), ya(0.0), yb(0.0);
		//y[3] = {-1.0, 0.0, 1.0}, t(0.0), h(0.01);
	valarray<double> y(11);
	//matrix<double> vz(3,11);
	
  	cout.setf(ios::fixed);	//輸出數(shù)據(jù)為定點法
	cout.precision(6);		//精度6位

	ODE_LinearBoundaryValude(a,b,ya,yb,11,y);

	for(int i=0; i<11; i++)
		cout << "    y(" << i << ") = " << y[i] << endl;
	cout << endl;
}

//計算二階微分方程組函數(shù)
template <class _Ty>
void FunctionValueLBV(_Ty x, valarray<_Ty>& z)
{
	z[0]=-1.0; 
	z[1]=0.0;
    z[2]=2.0/(x*x);
	z[3]=1.0/x;
}