//EigenvalueVectorRealSymmetryJacobi.cpp  
//實對稱陣特征值及特征向量雅可比法

#include <iostream>				//輸入輸出流頭文件
#include "Matrix.h"				//矩陣類及相關函數頭文件
#include "EigenvalueVector.h"	//特征值及特征向量頭文件
using namespace std;			//名字空間

void main()						// 定義控制臺應用程序的入口點
{
	double a[3][3] = 
	{
		{ 2.0, -1.0,  0.0},
		{-1.0,  2.0, -1.0},
		{ 0.0, -1.0,  2.0}
	};
	
	matrix<double> v(3,3), da(&a[0][0],3,3);
	double eps = FLOATERROR;
	
  	cout.setf(ios::fixed);	//輸出數據為定點法
	cout.precision(6);		//精度6位

	int k = EigenvalueVectorRealSymmetryJacobi(da, v, eps, 100);
	if(k > 0)
	{
		cout << "MATRIX A IS: " << endl;
		MatrixLinePrint(da);
		cout << endl;

		cout << "MATRIX V IS: " << endl;
		MatrixLinePrint(v);
		cout << endl;
	}

	double dmb[3][3] = 
	{
		{  3.0, -3.0, -2.0 },
		{ -3.0,  8.0,  4.0 },
		{  -2.0, 4.0,  3.0 } 
	};

}