//QR方法求矩阵特征值
/******
数据
-1 2 1
2 4 -1
1 1 -6
******/
#include
#include
#include
using namespace std;
#include
#define N 3 //矩阵的维数
#define NUM 15 //QR分解次数
#define SNUM 13 //用来控制输出格式
void Print(long double A[N][N],long double Q[N][N],long double R[N][N]); //矩阵输出
void QR(long double A[N][N],long double Q[N][N],long double R[N][N]); //QR分解
void Multiplicate(long double A[N][N],long double R[N][N],long double Q[N][N]); //迭代,获得下一次矩阵A=QR
int main()
{
int i,j;
long double A[N][N];
long double Q[N][N]={0};
long double R[N][N]={0};
cout<
cout<
cout<
for(i=0;i
for(j=0;j
cin>>A[i][j];
cout<
cout<
cout<
for(i=1