{

TH1F *h1;   

int ncx=100; 
h1 = new TH1F("h1","h1",ncx,0,ncx);


float x[ncx],y[ncx],z[ncx];
for (int i=0; i<ncx;i++) {
  x[i]= i+0.5;
  y[i]= ( x[i]- ncx *0.5 )*( x[i]- ncx *0.5 )/20;
  z[i]=100*sin( x[i]/10)*sin( x[i]/10);
}

TGraph *gr1 = new TGraph(ncx,x,y);
TGraph *gr2 = new TGraph(ncx,x,z);
h1->GetListOfFunctions()->Add(gr1); 
h1->GetListOfFunctions()->Add(gr2); 
TF1 *f1 = new TF1("f1","[0]*cos(x/[1])*sin(x/5)",0,100);
f1->SetParameter(0,100);
f1->SetParameter(1,10);

h1->GetListOfFunctions()->Add(f1); 

float xx;
for (int i=0; i<ncx;i++) {
  
  xx = h1->GetBinCenter(i);
  h1->Fill(xx,xx);
}

h1->Draw();
return;
}