le_toeplitzeuationlevinson.cpp

来自「计算机常用数值算法与程序(C++版) Chap07 线性方程组求解」· C++ 代码 · 共 26 行

//LE_ToeplitzEuationLevinson().cpp   
//求解托伯利兹方程组的列文逊方法 

#include <iostream>				//输入输出流头文件
#include "LinearEquation.h"		//线性方程(组)求解头文件

void main()
{
	int i;
    double t[6] = {6.0, 5.0, 4.0, 3.0, 2.0, 1.0};
    double b[6] = {11.0, 9.0, 9.0, 9.0, 13.0, 17.0};
    
    valarray<double> xx(6);
	valarray<double> tt(t, 6);
	valarray<double> bb(b, 6);
		
	if (LE_ToeplitzEuationLevinson(tt, bb, xx) > 0)
		for (i = 0; i < 6; i ++)
		{
			cout << "x(" << i << ") = " << xx[i] << endl;
		}

	cout << endl;    
}