ode_eulercontentstep.c

来自「给出计算常微分方程（组）的各种算法的使用示例。」· C语言 代码 · 共 39 行

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

#include <iostream>							//输入输出流头文件
#include "OrdinaryDifferentialEguation.h"	//求解常微分方程(组)头文件
using namespace std;						//名字空间

void main(void)
{	
	double y[3] = {-1.0, 0.0, 1.0}, t(0.0), h(0.01);
	valarray<double> vy(y, 3);
	matrix<double> vz(3,11);

	ODE_EulerContentStep(t,vy,h,11,vz);
	
  	cout.setf(ios::fixed);	//输出数据为定点法
	cout.precision(6);		//精度6位

	for(int i=0; i<11; i++)
	{
		double x = h * i;
		cout << "t = " << x << endl;

		for(int j=0; j<3; j++)
			cout << "    y(" << j << ") = " << vz(j, i) << "\t";
		cout << endl;
    }
}

template <class _Ty>
void FunctionValueECS(_Ty t, valarray<_Ty>& y, valarray<_Ty>& d)
{
    d[0]=y[1];
	d[1]=-y[0]; 
	d[2]=-y[2];

}