Giải thuật sinh đường tròn Midpoint

Phương trình đường tròn không tường minh:
f(x,y) = x^2+y^2-R^2 =0
Nếu f(x,y) = 0 thì nằm trên đường tròn
f(x,y) > 0 thì nằm bên ngoài đường tròn
f(x,y) < 0thì nằm bên trong đường tròn
Thực hiện giải thuật trên 1/8 đường tròn và lấy đối xứng cho các góc còn lại.
Với M là điểm giữa của AB
Với di là giá trị của đường tròn tại một điểm bất kỳ



Cài đặt
void Mid_circle(int xc, int yc, int Radius, int
color)
{
int x, y, d;
x = 0;
y = Radius;
d = 1- Radius;
while (x <= y)
{
putpixel(xc + x, yc + y, color);
if (d< 0)
d +=2 * x + 3;
else
{
d += 2 * (x-y) + 5;
y--;
}
x++;
}
}