{
  double rmin=0;
  double rmax=TMath::Pi()*2;
  double r1[1000];
  double r2[1000];
  double r3[1000];
  double theta[1000];
  for(int i=0; i<1000; ++i) {
    theta[i] = i*(rmax-rmin)/1000 + rmin;
    r1[i]=cos(theta[i]);
    r2[i]=sin(theta[i]);
    r3[i]=1;
  }
  TCanvas *Cpol = new TCanvas("Cpol1", "TGraphPolar Example: cirles...", 600,600) ;
  TGraphPolar *gr1 = new TGraphPolar (1000, theta, r1) ;
  TGraphPolar *gr2 = new TGraphPolar (1000, theta, r2) ;
  TGraphPolar *gr3 = new TGraphPolar (1000, theta, r3) ;

  gr1->SetTitle("") ;
  gr1->SetLineColor(2) ;
  gr1->SetLineWidth(3) ;
  gr1->Draw ("AOL");

  gr2->SetTitle("") ;
  gr2->SetLineColor(3) ;
  gr2->SetLineWidth(3) ;
  gr2->Draw ("AOLsame");

  gr3->SetTitle("three      circles") ;
  gr3->SetLineColor(4) ;
  gr3->SetLineWidth(3) ;
  gr3->Draw ("AOLsame");



  //----------------------------------------
  // Προσθέστε αυτές τις εντολές για να σχεδιαστεί σωστά η ακτίνα (από 0 και πάνω...)
  Cpol->Draw();
  gr1->SetMinRadial(0);
  gr1->SetMaxRadial(1.1);  // <-- βάλτε τη σωστή τιμή του ρ_max
  gr1->SetMinPolar(0);
  gr1->SetMaxPolar(2*TMath::Pi()); 
  gr2->SetMinRadial(0);
  gr2->SetMaxRadial(1.1);  // <-- βάλτε τη σωστή τιμή του ρ_max
  gr2->SetMinPolar(0);
  gr2->SetMaxPolar(2*TMath::Pi()); 
  gr3->SetMinRadial(0);
  gr3->SetMaxRadial(1.1);  // <-- βάλτε τη σωστή τιμή του ρ_max
  gr3->SetMinPolar(0);
  gr3->SetMaxPolar(2*TMath::Pi()); 
  //----------------------------------------

}