{
  double phimin=0;
  double phimax=TMath::Pi();

  double phi[1000];
  double   r[1000];

  TF1 *fp1 = new TF1("fp1", "cos(x)", phimin, phimax) ;
  for(int i=0; i<1000; ++i) {
    phi[i] = i*(phimax-phimin)/1000 + phimin; // γωνία
    r[i]   = fp1->Eval(phi[i]);               // ακτίνα
  }

  TCanvas *Cpol = new TCanvas("Cpol", "TGraphPolar Example", 600,600);

  TGraphPolar *gr = new TGraphPolar(1000, phi, r);
  gr->SetTitle("cos(#varphi) in polar coordinates: a circle!");
  gr->SetLineColor(kRed);
  gr->SetLineWidth(3);
  gr->Draw ("AOL");

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