le_toeplitzeuationlevinson.cpp

vc++实现线性方程组求解 1全选主元高斯消元法 2全选主元高斯-约当消元法 3三对角方程组的追赶法 4一般带型方程组求解 5对称方程组的分解法 6对称正定方程组的平方根法

//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;    
}