#include <iostream>

#include "TLinearFitter.h"
#include "TVector.h"

using namespace std;

void test()
{
  Double_t a[] = {1.0, 3.0, 4.0, 2.0, 0.0, 3.0, 0.0, 1.0, 1.0, 1.0, 1.0, 0.0};
  Double_t y[] = {2.0, 1.0, 0.0, 8.0};

  TLinearFitter *lf = new TLinearFitter( 1 );
  lf->StoreData(1);
  lf->SetFormula( "x[0]++x[1]++x[2]" );
  lf->AssignData(4, 3, a, y);
  lf->EvalRobust();

  TVectorD result(3);
  lf->GetParameters(result);
  cout << "results: " << result[0] << " " << result[1] << " " << result[2] << endl;
}