Runge-Kutta-Simulation/AppleX/GRAPHICS/LINE.C

int plot(), bplot(), gplot(), rplot(), vplot();

drw(x1,y1,x2,y2)
int x1,y1,x2,y2;
{
	comline(x1,y1,x2,y2,plot);
}

bdrw(x1,y1,x2,y2)
int x1,y1,x2,y2;
{
	comline(x1,y1,x2,y2,bplot);
}

gdrw(x1,y1,x2,y2)
int x1,y1,x2,y2;
{
	comline(x1,y1,x2,y2,gplot);
}

rdrw(x1,y1,x2,y2)
int x1,y1,x2,y2;
{
	comline(x1,y1,x2,y2,rplot);
}

vdrw(x1,y1,x2,y2)
int x1,y1,x2,y2;
{
    comline(x1,y1,x2,y2,vplot);
}


int __oldx=0, __oldy=0;

lineto(x,y)
int x,y;
{
	drw(__oldx,__oldy,x,y);
}

blineto(x,y)
int x,y;
{
	bdrw(__oldx,__oldy,x,y);
}

glineto(x,y)
int x,y;
{
   	gdrw(__oldx,__oldy,x,y);
}

vlineto(x,y)
int x,y;
{
	vdrw(__oldx,__oldy,x,y);
}

rlineto(x,y)
int x,y;
{
    rdrw(__oldx,__oldy,x,y);
}


static int
comline(x1,y1,x2,y2,func)
int (*func)();
{
	register int dx,dy,ix,iy;
	int e,ei,ed,i;

	__oldx=x2; __oldy=y2;
	dx=x2-x1;
	dy=y2-y1;
	ix=1;
	if(dx<0)
		{ ix=-1;
		  dx=-dx;
		}
	iy=1;
	if(dy<0)
		{ iy=-1;
		  dy=-dy;
		}
	if(dy>dx)
		goto ylin;

	ei=2*dy;
	ed=ei-2*dx;
	e=-dx+ei;
	for(i=1;i<=dx;++i){
		(*func)(x1,y1);
		x1=x1+ix;
		if(e<0)
			e=e+ei;
		else{ 
			y1=y1+iy;e=e+ed;
		}
	}
	goto fin;

ylin:

	ei=2*dx;
	ed=ei-2*dy;
	e=-dy+ei;
	for(i=1;i<=dy;++i)
		{ (*func)(x1,y1);
		  y1=y1+iy;
		  if(e<0)
			e=e+ei;
		  else{ x1=x1+ix;e=e+ed;}
		}
fin:		  
;		   
}