void k1() //0<k<1 
{ 
 glClear(GL_COLOR_BUFFER_BIT); 
 
 glColor3f(0.0,0.0,0.0); 
 glBegin(GL_POINTS); 
 GLint x1=0,y1=0,x2=400,y2=200; 
 GLint x=x1,y=y1; 
 GLint dx=x2-x1,dy=y2-y1,dT=2*(dy-dx),dS=2*dy; 
 GLint d=2*dy-dx; 
 glVertex2i(x,y); 
 while(x<x2) 
 { 
 x++; 
 if(d<0) 
 d=d+dS; 
 else 
 { 
 y++; 
 d=d+dT; 
 } 
 glVertex2i(x,y); 
 } 
 glEnd(); 
 glFlush(); 
 
}

區(qū)域填充
#include "gl/glut.h"
#include "windows.h"
const int POINTNUM=7; //多邊形點(diǎn)數(shù).
 
/******定義結(jié)構(gòu)體用于活性邊表AET和新邊表NET***********************************/
 typedef struct XET
 {
 float x;
 float dx,ymax;
 XET* next;
 }AET,NET;
 
/******定義點(diǎn)結(jié)構(gòu)體point******************************************************/
 struct point
 {
 float x;
 float y;
 }
 polypoint[POINTNUM]={250,50,550,150,550,400,250,250,100,350,100,100,120,30};//多邊形頂點(diǎn)
 
 void PolyScan()
{
/******計(jì)算最高點(diǎn)的y坐標(biāo)(掃描到此結(jié)束)****************************************/
 int MaxY=0;
 int i;
 for(i=0;i<POINTNUM;i++)
 if(polypoint[i].y>MaxY)
 MaxY=polypoint[i].y;
 
/*******初始化AET表***********************************************************/
 AET *pAET=new AET;
 pAET->next=NULL;
 
/******初始化NET表************************************************************/
 NET *pNET[1024];
 for(i=0;i<=MaxY;i++)
 {
 pNET[i]=new NET;
 pNET[i]->next=NULL;
 }
 glClear(GL_COLOR_BUFFER_BIT); //賦值的窗口顯示. 
 glColor3f(0.0,0.0,0.0);  //設(shè)置直線的顏色紅色
 glBegin(GL_POINTS);
/******掃描并建立NET表*********************************************************/
 for(i=0;i<=MaxY;i++)
 {
 for(int j=0;j<POINTNUM;j++)
 if(polypoint[j].y==i)
 { //一個點(diǎn)跟前面的一個點(diǎn)形成一條線段，跟后面的點(diǎn)也形成線段 
 if(polypoint[(j-1+POINTNUM)%POINTNUM].y>polypoint[j].y)
 {
 NET *p=new NET;
 p->x=polypoint[j].x;
 p->ymax=polypoint[(j-1+POINTNUM)%POINTNUM].y;
 p->dx=(polypoint[(j-1+POINTNUM)%POINTNUM].x-polypoint[j].x)/(polypoint[(j-1+POINTNUM)%POINTNUM].y-polypoint[j].y);
 p->next=pNET[i]->next;
 pNET[i]->next=p;
 
 }
 if(polypoint[(j+1+POINTNUM)%POINTNUM].y>polypoint[j].y)
 {
 NET *p=new NET; 
 p->x=polypoint[j].x;
 p->ymax=polypoint[(j+1+POINTNUM)%POINTNUM].y;
 p->dx=(polypoint[(j+1+POINTNUM)%POINTNUM].x-polypoint[j].x)/(polypoint[(j+1+POINTNUM)%POINTNUM].y-polypoint[j].y);
 p->next=pNET[i]->next;
 pNET[i]->next=p;
 }
 }
 }
/******建立并更新活性邊表AET*****************************************************/
for(i=0;i<=MaxY;i++)
 {
 //計(jì)算新的交點(diǎn)x,更新AET
 NET *p=pAET->next;
 while(p)
 {
 p->x=p->x + p->dx;
 p=p->next;
 }
 //更新后新AET先排序*************************************************************/
 //斷表排序,不再開辟空間
 AET *tq=pAET;
 p=pAET->next;
 tq->next=NULL;
 while(p)
 {
 while(tq->next && p->x >= tq->next->x)
 tq=tq->next;
 NET *s=p->next;
 p->next=tq->next;
 tq->next=p;
 p=s;
 tq=pAET;
 }
 //(改進(jìn)算法)先從AET表中刪除ymax==i的結(jié)點(diǎn)****************************************/
 AET *q=pAET;
 p=q->next;
 while(p)
 {
 if(p->ymax==i)
 {
 q->next=p->next;
 delete p;
 p=q->next;
 }
 else
 {
 q=q->next;
 p=q->next;
 }
 }
 //將NET中的新點(diǎn)加入AET,并用插入法按X值遞增排序**********************************/
 p=pNET[i]->next;
 q=pAET;
 while(p)
 {
 while(q->next && p->x >= q->next->x)
 q=q->next;
 NET *s=p->next;
 p->next=q->next;
 q->next=p;
 p=s;
 q=pAET;
 }
/******配對填充顏色***************************************************************/
 
  p=pAET->next;
 while(p && p->next)
 {
 for(float j=p->x;j<=p->next->x;j++)
 glVertex2i(static_cast<int>(j),i);
 p=p->next->next;//考慮端點(diǎn)情況
 }
 
 
 }
 glEnd();
glFlush(); 
}
void init(void)
{glClearColor(1.0,1.0,1.0,0.0);
//窗口的背景顏色設(shè)置為白色
glMatrixMode(GL_PROJECTION);
gluOrtho2D(0.0,600.0,0.0,450.0);
}
 
void main(int argc,char* argv)
{
 glutInit(&argc,&argv);  //I初始化 GLUT.
 glutInitDisplayMode(GLUT_SINGLE | GLUT_RGB); //設(shè)置顯示模式：單個緩存和使用RGB模型
 glutInitWindowPosition(50,100); //設(shè)置窗口的頂部和左邊位置
 glutInitWindowSize(400,300); //設(shè)置窗口的高度和寬度
 glutCreateWindow("An Example OpenGL Program"); //創(chuàng)建顯示窗口
 
 init();    //調(diào)用初始化過程
 glutDisplayFunc(PolyScan); //圖形的定義傳遞給我window.
 glutMainLoop();   //顯示所有的圖形并等待
 }