{
  double f=1. ; // frequenza
  double omega=2.*TMath::Pi()*f; // omega
  double x0=0., v0=1.; // condizioni iniziali
  double E = 0.5*(v0*v0+omega*omega*x0*x0); // energia dell'oscillatore
  double A = sqrt(2.*E)/omega; // ampiezza dell'oscillazione
  double Phi = atan2(v0,x0*omega);   // fase dell'oscillazione al tempo 0
  double t,x,v;
  TGraph grafico;
  int punto=0;
  ifstream in("pippo3.dat");
  while (1) {
    in >> t >> x >> v;
    if ( in.eof() ) break;
    double err = fabs(x-A*cos(omega*t-Phi));
    grafico.SetPoint(punto,t,err);
    punto++;
  }
  in.close();
  grafico.Draw("AL");
}