ode_linearboundaryvalude.cpp

何光渝所著《计算机常用数值算法与程序(C++版)》的随书源代码。 是想学习C++科学计算的人士的必看代码。

//ODE_LinearBoundaryValude.cpp	
//定步长全区间积分的欧拉法

#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);	//输出数据为定点法
	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;
}

//计算二阶微分方程组函数
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;
}