{
  // Creazione di una canvas con griglia:
  gROOT->Reset();
  c1 = new TCanvas("c1", "Fit dei dati con errori su x e y", 200, 10, 700, 500);
  c1->SetGrid();
  Int_t n = 10;
  // Inizializzazione di 4 vettori:
  // x[] e ex[] sono le misure in ascissa x_i e i loro errori
  Float_t x[n]  = {0.0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.65, 0.75, 0.95, 1.2};
  Float_t ex[n] = {0.05, 0.05, 0.04, 0.05, 0.04, 0.05, 0.06, 0.06, 0.05, 0.05};
  // y[] e ey[] sono le misure in ordinata y_i e i loro errori
  Float_t y[n]  = {0.1, 0.3, 0.55, 0.67, 0.75, 0.79, 0.71, 0.63, 0.45, 0.35};
  Float_t ey[n] = {0.04, 0.03, 0.03, 0.02, 0.02, 0.02, 0.03, 0.03, 0.04, 0.04};
  // Creazione di un grafico x-y con errori
  gr = new TGraphErrors(n,x,y,ex,ey);
  gr->SetTitle("Fit della funzione");
  gr->SetMarkerColor(4);
  gr->SetMarkerStyle(21);
  gr->Draw("AP");
  gStyle->SetOptFit(1111);
  // Definizione di una nuova funzione
  TF1 *funz = new TF1("funz","[0]*x*exp([1]*x+[2])",0.5,100);
  funz->SetParameters(5.,-1.,0.);
  // Fit in un intervallo di interesse, tra -3 e 7
  gr->Fit("funz","","",-3,7);
  c1->Update();

  c1->Update();
}