void MedieXY()
{
  // carica la libreria
  gSystem->Load("LibEsp3_C.so");

  // legge i dati e crea il grafico con errori
  TGraphErrors *g=LoadGraphErrors("CalibPotR6.dat");

  // recupera il numero di punti
  Int_t N=g->GetN();
  
  // calcola la media
  Double_t mx=0, my=0;
  for (Int_t i=0; i<N; i++) {
    mx += g->GetX()[i];
    my += g->GetY()[i];
  }
  mx /= N;
  my /= N;
  
  // stampa i risultati
  printf("media delle X = %f\nmedia delle Y = %f\n",mx,my);
}