中点算法之y=x^2/16

/*************************************************************
pb-图形学题3
中点算法之y=x^2/16;

*************************************************************/

#include <GL/glut.h>
#include<cstdio>
#include<cmath>

void init()
{
glutInitDisplayMode(GLUT_RGB | GLUT_SINGLE);                 //什么单缓存的东西，现在不懂什么意思，就当通用条件
glutInitWindowPosition(0, 0);                                //图片出现的位置
glutInitWindowSize(400, 300);                                //图片的长和宽
glutCreateWindow("pb-图形学题1");                            //图片的名字

glClearColor(0.0, 0.0, 0.0, 0.0);
glMatrixMode(GL_PROJECTION);
gluOrtho2D(-100, 100, -10, 200);

glRectf(-0.5f, -0.5f, 0.5f, 0.5f);
glClear(GL_COLOR_BUFFER_BIT);
glColor3f(1.0, 0.0, 0.0);                                    //线的颜色
}

void setPixel (GLint x,GLint y,GLint nx,GLint ny)               //描点函数，表示从x,y到nx,ny的矩形全部涂色
{
glBegin(GL_LINES);
glVertex2i(x, y);
glVertex2i(nx,ny);
glEnd();
}

void myDisplay()
{
int d=-7;
int x1=0,x2=0;
int y1=0,y2=0;
while (x1<=8)
{
x2=x1+1;
if (d<=0)
{
d+=2*x2+1;
y2=y1;
}
else
{
d+=2*x2-15;
y2=y1+1;
}
setPixel(x1,y1,x2,y2);
setPixel(-x1,y1,-x2,y2);
y1=y2;
x1=x2;
}
d=8;
while (y1<=100)
{
y2=y1+1;
if (d>0)
{
x2=x1+1;
d+=16-(x1*x1-(x1-1)*(x1-1)+1);
}
else
{
x2=x1;
d+=16;
}
setPixel(x1,y1,x2,y2);
setPixel(-x1,y1,-x2,y2);
y1=y2;
x1=x2;
}
glFlush();
}

int main(int argc, char *argv[])
{
glutInit(&argc, argv);
init();                                                     //初始化数据
glutDisplayFunc(&myDisplay);                                //调用函数
glutMainLoop();                                             //开始程序
return 0;
}