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

void diago(){
  TMatrixD m(3,3);
  m[0][0] =  1.10866  ; m[0][1] = -3.1856000; m[0][2] = -1.0835000;
  m[1][0] = -3.1856000; m[1][1] = -10.715299; m[1][2] = -1.3247000;
  m[2][0] = -1.0835000; m[2][1] = -1.3247000; m[2][2] = -8.4516797;

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