//g++ `root-config --cflags --glibs` -o main main.cpp
#include <iostream>
#include "TROOT.h"
#include "TMatrixT.h"
#include "TVectorT.h"
#include "TDecompSVD.h"
#include "RVersion.h"

using namespace std;

int main(int argc, char** argv)
{	
	cout << ROOT_RELEASE << endl;
	

	double a = 1.689339828220179e-03;
	double b = 1.689339828220178e-03;
	double c = 1.5e-4;
	TMatrixD m(4, 3);
		
	m(0,0) =  a; m(0,1) =  a; m(0,2) = -c;
	m(1,0) =  a; m(1,1) =  a; m(1,2) =  c;
	m(2,0) = -b; m(2,1) = -b; m(2,2) =  c;
	m(3,0) = -b; m(3,1) = -b; m(3,2) = -c;
	
	m.Print();
	
	
	TDecompSVD s(m);
	
	s.Decompose();
	s.Print();
	
	return 0; 
}