#include <TMatrixDSym.h>
#include <TMatrixDSymEigen.h>
#include <TVectorD.h>

void matrix2(){

  const int N = 2;
  
  double e[N*N] = {
    1., 2.,
    3., 4.
  };

  TMatrixDSym m(N, e);
  TMatrixDSymEigen me(m);

  TVectorD eigenval = me.GetEigenValues();
  TMatrixD eigenvec = me.GetEigenVectors();

  m.Print();
  eigenval.Print();
  eigenvec.Print();

}