int rutherford ()
{
const Int_t dim = 185;
const Int_t nrighe = 21;
Double_t x1[nrighe], ex1[nrighe], y1[nrighe], ey1[nrighe];
  
x1[0]=0.000000;ex1[0]=0;y1[0]=422.800000;ey1[0]=10.851728;
x1[1]=5.000000;ex1[1]=0;y1[1]=392.833333;ey1[1]=15.410134;
x1[2]=10.000000;ex1[2]=0;y1[2]=280.500000;ey1[2]=9.322911;
x1[3]=15.000000;ex1[3]=0;y1[3]=89.142857;ey1[3]=9.760395;
x1[4]=20.000000;ex1[4]=0;y1[4]=10.375000;ey1[4]=3.425547;
x1[5]=25.000000;ex1[5]=0;y1[5]=3.222222;ey1[5]=1.930905;
x1[6]=30.000000;ex1[6]=0;y1[6]=1.100000;ey1[6]=0.943398;
x1[7]=35.000000;ex1[7]=0;y1[7]=0.250000;ey1[7]=0.433013;
x1[8]=40.000000;ex1[8]=0;y1[8]=0.266667;ey1[8]=0.573488;
x1[9]=50.000000;ex1[9]=0;y1[9]=0.272727;ey1[9]=0.616575;
x1[10]=60.000000;ex1[10]=0;y1[10]=0.000000;ey1[10]=0.000000;
x1[11]=-5.000000;ex1[11]=0;y1[11]=290.800000;ey1[11]=4.791659;
x1[12]=-10.000000;ex1[12]=0;y1[12]=101.142857;ey1[12]=8.708429;
x1[13]=-15.000000;ex1[13]=0;y1[13]=13.000000;ey1[13]=2.984810;
x1[14]=-20.000000;ex1[14]=0;y1[14]=3.733333;ey1[14]=2.048306;
x1[15]=-25.000000;ex1[15]=0;y1[15]=1.200000;ey1[15]=0.979796;
x1[16]=-30.000000;ex1[16]=0;y1[16]=0.777778;ey1[16]=1.030402;
x1[17]=-35.000000;ex1[17]=0;y1[17]=0.100000;ey1[17]=0.300000;
x1[18]=-40.000000;ex1[18]=0;y1[18]=0.000000;ey1[18]=0.000000;
x1[19]=-50.000000;ex1[19]=0;y1[19]=0.200000;ey1[19]=0.600000;
x1[20]=-60.000000;ex1[20]=0;y1[20]=0.000000;ey1[20]=0.000000;

TGraphErrors *gr1 = new TGraphErrors(nrighe,x1,y1,ex1,ey1);
TF1 *ruth = new TF1("ruth", "[0]/((sin(x+[1]))^4)",-60,60);
gr1->Fit("ruth");
gr1->SetMarkerStyle(21);
gr1->SetTitle("Rutherford scattering");
gr1->GetXaxis()->SetTitle("#theta");
gr1->GetXaxis()->CenterTitle();
gr1->GetYaxis()->SetTitle("rate");
gr1->GetYaxis()->CenterTitle();
gr1->Draw("AP");

return 0;

}