{

/*
TH1F *h1 = new TH1F("h1","test",3,0,3);
h1->Fill(0.5,1);
h1->Fill(1.5,5);
h1->Fill(2.5,3);
*/

float x[11] = {2.00,2.10,2.20,2.30,2.40,2.50,2.60,2.70,2.80,2.90,3.00};
float y[11] = {2.78,3.29,3.29,3.33,3.23,3.69,3.46,3.87,3.62,3.40,3.99};
float ex[11];
float ey[11];
for (int i=0;i<11;++i) {
  ey[i] = 0.3;
  ex[i] = 0.0; // This is needed!
}



TH2F *hframe = new TH2F("hframe","",300,0,4,300,0,5);
hframe->Draw();

TGraphErrors *g = new TGraphErrors(11,x,y,ex,ey);
g->SetMarkerColor(4);
g->SetMarkerStyle(21);
g->SetMarkerSize(.7);

TF1 *f1 = new TF1("f1","[0]+[1]*x",0,5);

TMinuit *gMinuit = new TMinuit(2);
g->Fit("f1","R");
//g->Fit("pol1");

g->Draw("P");
Double_t matrix[2][2];
Int_t ndim=2;
gMinuit->mnemat(&matrix[0][0],ndim);

cout << "matrix[0][0]: " << matrix[0][0] << endl;

}