void grafico(){

  TGraphErrors *gl=new TGraphErrors("output.dat","%lg %lg");
  TF1 *retta0=new TF1("retta0","[0]*sin([2]*x+[1])+[3]",1.,1024.);
  
  gl->Draw("AP");
  retta0->Fit("retta0");
  retta0->SetParameters(0.24,1.,1.,1.);
  retta0->Draw("same");
  
}