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

void diago_1e(){
  TMatrixD m(2,2);
  m[0][0] = 1.0334001; m[0][1] = 0.1874000;
  m[1][0] = 0.1874000; m[1][1] = 0.6066000;

  m[0][0] += 1.6465800-3.9477999;
  m[1][1] += 1.6465800-3.1635399;

  std::cout << m[0][0] << std::endl;
  std::cout << m[1][1] << std::endl;

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