le_symmetryregulareuationsquareroot.cpp

计算机常用数值计算算法与程序 C++版 各章代码

//LE_SymmetryRegularEuationSquareRoot.cpp   
//求解对称正定方程组的平方根法

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

void main()
{ 
	double a[4][4] = 
	{
		{5.0,  7.0,   6.0,   5.0},
		{7.0, 10.0,   8.0,   7.0},
		{6.0,  8.0,  10.0,   9.0},
		{5.0,  7.0,   9.0,  10.0}
	};
    
	double d[4][2] = 
	{
		{23.0,  92.0},
		{32.0, 128.0},
        {33.0, 132.0},
		{31.0, 124.0}
	};
    
	matrix<double> aa(&a[0][0], 4, 4);
	matrix<double> dd(&d[0][0], 4, 2);
		
  	cout.setf(ios::fixed);	//输出数据为定点法
	cout.precision(6);		//精度6位
	
	if(LE_SymmetryRegularEuationSquareRoot(aa, dd) > 0)
	{
		cout << endl;
		cout.width(22);
		cout << "Group 1";
		cout.width(16);
		cout << "Group 2" << endl << endl;

		for(int i = 0; i < 4; i ++)
		{
			cout << " x(" << i << ") = ";
			MatrixLinePrint(dd, i);
		}
	}
	cout << endl;
}