void DrawClipCircle(short x, short y, short r)
{
  short Xp=r;
  short Yp=0;
  Short Err=r/2;
 
  do
  {
     DrawClipPix (x+Xp,y+Yp);
     DrawClipPix (x-Xp,y+Yp);
     DrawClipPix (x+Xp,y-Yp);
     DrawClipPix (x-Xp,y-Yp);
 
     DrawClipPix (x+Yp,y+Xp);
     DrawClipPix (x-Yp,y+Xp);
     DrawClipPix (x+Yp,y-Xp);
     DrawClipPix (x-Yp,y-Xp);
 
     Yp++;
     Err += Yp;
     If (Err>0)
     {
       Xp--;
       Err -= Xp;
     }
  } while (Xp<Yp);

//==========================================================================//
//  Utilise les symétries du cercle pour en dessiner un plein
//==========================================================================//

void DessineLignes(void * screen_0, void * screen_1, short color, short x, short y, short x_centre, short y_centre)
{
 GrayFastDrawHLine2B(y_centre+y, y_centre-y, x_centre+x, color, screen_0, screen_1);
 GrayFastDrawHLine2B(y_centre+y, y_centre-y, x_centre-x, color, screen_0, screen_1);
 GrayFastDrawHLine2B(y_centre+x, y_centre-x, x_centre+y, color, screen_0, screen_1);
 GrayFastDrawHLine2B(y_centre+x, y_centre-x, x_centre-y, color, screen_0, screen_1);
}

//==========================================================================//
// algo de Brensenham

void DessineCerclePlein(void * screen_0, void * screen_1, short color, short rayon, short x_centre, short y_centre)
{
 register short x, y, d;

 x = 0;
 y = rayon;
 d = (5-rayon)/4;
 
 DessineLignes(screen_0, screen_1, color, x, y, y_centre, x_centre);
 
 while(y > x)
 {
   if(d < 0)
   {
     d += x+x+3;
   }
   else
   {
     d += x+x-y-y+5;
     y --;
   }
   
   x ++;
   DessineLignes(screen_0, screen_1, color, x, y, y_centre, x_centre);
 }
 
} 

//==========================================================================//

// Main Function
void _main(void)
{
  register short i;
  
  ClrScr();
  GrayOn();
  
  for(i = 0; i < 50; i ++)
  {
    DessineCerclePlein(__L_plane, __D_plane, COLOR_BLACK, i, 80, 50);  
  }
  
  ngetchx();
  
  GrayOff();
}

//==========================================================================//

Cls
Input "X";X%
Input "Y";Y%
Input "Rayon";R%
Cls
I%=X%-R%
While I%<=X%+R%
  J%=Y%-R%
  While J%<=Y%+R%
    Tmp%=(I%-X%)^2+(J%-Y%)^2
    If Tmp%<R%*R%+R% And Tmp%>R%*R%-R%   !Pour un cercle vide
    If Tmp%<=R%^2  !Pour un cercle plein
      Plot I%,J%
    Endif
    Inc J%
  Wend
  Inc I%
Wend