Алгоритм Брезенхема для генерации окружности C# OpenGl

using System;
using System.Collections.Generic;
using System.ComponentModel;
using System.Data;
using System.Drawing;
using System.Linq;
using System.Text;
using System.Threading.Tasks;
using System.Windows.Forms;
using Tao.OpenGl;
using Tao.Platform.Windows;
using Tao.FreeGlut;
 
namespace lab1
{
    public partial class Form1 : Form
    {
        public Form1()
        {
            InitializeComponent();
            AnT.InitializeContexts();
            
        }
 
        private void Form1_Load(object sender, EventArgs e)
        {
          
            Glut.glutInit();
            Glut.glutInitDisplayMode(Glut.GLUT_RGB | Glut.GLUT_SINGLE);
            Gl.glClearColor(255, 255, 255, 1); 
            Gl.glViewport(0, 0, AnT.Width, AnT.Height);
            Gl.glMatrixMode(Gl.GL_PROJECTION);
            Gl.glLoadIdentity();
            Glu.gluOrtho2D(0.0, AnT.Height, 0.0, AnT.Width);
            Gl.glMatrixMode(Gl.GL_MODELVIEW);
            Gl.glLoadIdentity();
        }
 
        private void button3_Click(object sender, EventArgs e)
        {
            Application.Exit();
        }
 
        private void button1_Click(object sender, EventArgs e)
        {
            int x = System.Convert.ToInt32(CircleX.Text);
            int y = System.Convert.ToInt32(CircleY.Text);
            int r = System.Convert.ToInt32(CircleR.Text);
            pCil.Text = System.Convert.ToString(Math.Round(Math.Sqrt(Math.Pow(x, 2) + Math.Pow(y, 2)), 3));
            tgCil.Text = System.Convert.ToString(Math.Round(Math.Atan(y / x), 3));
            ZCil.Text = "0.000";
            int yd = r;
            int xd = 0;
            int d = 2-2 *r;
            Gl.glPointSize(1);
            Gl.glBegin(Gl.GL_POINTS);
                      
 
            while (xd <= yd)
            {
                Gl.glVertex3d(x + xd, y + yd, 0);
                Gl.glVertex3d(x - xd, y + yd, 0);
                Gl.glVertex3d(x + xd, y - yd, 0);
                Gl.glVertex3d(x - xd, y - yd, 0);
                Gl.glVertex3d(x + yd, y + xd, 0);
                Gl.glVertex3d(x - yd, y + xd, 0);
                Gl.glVertex3d(x + yd, y - xd, 0);
                Gl.glVertex3d(x - yd, y - xd, 0);
              
                if (d < 0)
                    d += (2 * xd) + 3;
                else
                {
                    d += (2 * (xd - yd)) + 5;
                    yd -= 1;
                }
                xd++;
 
            }
 
            Gl.glEnd(); 
            Gl.glFlush();
            AnT.Invalidate();
        }
 
    }
}