{

  float *l_TabularX = new float[6];
  float *l_TabularDeltaX = new float[6];
  float *l_TabularY = new float[6];
  float *l_TabularDeltaY = new float[6];

  l_TabularX[0] = 0.325667;
  l_TabularDeltaX[0] = 0;
  l_TabularY[0] = 23.0751;
  l_TabularDeltaY[0] = 0;
  
  
  l_TabularX[1] = 1.04715;
  l_TabularDeltaX[1] = 0;
  l_TabularY[1] = 75.2812;
  l_TabularDeltaY[1] = 0;

  
  l_TabularX[2] = 0.46265;
  l_TabularDeltaX[2] = 0;
  l_TabularY[2] = 32.9052;
  l_TabularDeltaY[2] = 0;


  l_TabularX[3] = 0.684093;
  l_TabularDeltaX[3] = 0;
  l_TabularY[3] = 49.0763;
  l_TabularDeltaY[3] = 0;

  
  l_TabularX[4] = 1.59671;
  l_TabularDeltaX[4] = 0;
  l_TabularY[4] = 114.905;
  l_TabularDeltaY[4] = 0;

  int n = 5;
  
  TGraphErrors *gr = new TGraphErrors(n, l_TabularX , l_TabularY , l_TabularDeltaX, l_TabularDeltaY);

  float min = 0.325667;
  float max = 1.59671;
  
  TF1 *f1 = new TF1("f1", "[0]*x+[1]", min, max);
  //fit the graph with this function
  gr->Fit("f1");

  TVirtualFitter *fitter = TVirtualFitter::GetFitter(); 
  float Error0 = sqrt(fitter->GetCovarianceMatrixElement(0,0));
  float Error1 = sqrt(fitter->GetCovarianceMatrixElement(1,1));

  cout<<"compare "<<f1->GetParError(0)<<" et "<<Error0<<endl;
  cout<<"compare "<<f1->GetParError(1)<<" et "<<Error1<<endl;

}