#include "TMatrixD.h"
#include "TVectorD.h"
#include "TMatrixDEigen.h"

void diago5(){
  TMatrixD m(5,5);
  m[0][0] = -0.0665; m[0][1] = -0.6149; m[0][2] = -0.5154; m[0][3] = -1.6221; m[0][4] = -0.4041;
  m[1][0] = -0.6149; m[1][1] = -0.3248; m[1][2] = -0.5247; m[1][3] = -0.2828; m[1][4] = -0.4770;
  m[2][0] = -0.5154; m[2][1] = -0.5247; m[2][2] = -0.4064; m[2][3] = -0.6198; m[2][4] = -1.9410;
  m[3][0] = -1.6221; m[3][1] = -0.2828; m[3][2] = -0.6198; m[3][3] = -1.0020; m[3][4] = -0.8616;
  m[4][0] = -0.4041; m[4][1] = -0.4770; m[4][2] = -1.9410; m[4][3] = -0.8616; m[4][4] = -0.8183;
  
  m[0][0] += 2*1.6465800;
  m[1][1] += 1.6466-3.9477999;
  m[2][2] += 1.6466-3.1635399;
  m[3][3] += 2*(-3.9477999);
  m[4][4] += -3.9477999-3.1635399;

  TMatrixDEigen *x = new TMatrixDEigen(m);
  TVectorD v(5);
  v = x->GetEigenValues();
  
  m.Print();
  v.Print();
}