1:
2:
3:
4:
5:
6:
7:
8:
9:
10:
11:
12:
13:
14:
15:
16:
17:
18:
19:
20:
21:
22:
23:
24:
25:
26:
27:
28:
29:
30:
31:
32:
33:
34:
35:
36:
37:
38:
39:
40:
41:
42:
43:
44:
45:

procedure gbham(xstart,ystart,xend,yend:Integer;var npix:Integer;var koords:Tposarray);
var x, y, t, dist, xerr, yerr, dx, dy, incx, incy:Integer;
begin
// Entfernung in beiden Dimensionen berechnen
     dx:=xend - xstart;
     dy:=yend - ystart;
// Vorzeichen des Inkrements bestimmen
     if(dx<0)then begin
       incx:=-1;
       dx:=-dx;
     end else if (dx > 0) then incx:=1
     else incx:=0;
     if(dy<0) then begin
       incy:=-1;
       dy:=-dy;
     end else if (dy > 0) then incy:=1
     else incy:=0;
// feststellen, welche Entfernung größer ist
     dist:=dy;
     if(dx>=dy) then dist:=dx;
// Initialisierungen vor Schleifenbeginn
     x:=xstart;
     y:=ystart;
     xerr:=dx;
     yerr:=dy;
     npix:=dist+1;
     Setlength(koords,npix);
// Pixel berechnen
     for t:=0 to dist-1 do begin
       koords[t][0]:=x;
       koords[t][1]:=y;
       xerr:=xerr+dx;
       yerr:=yerr+dy;
       if(xerr>dist) then begin
         xerr:=xerr-dist;
         x:=x+incx;
       end;
       if(yerr>dist) then begin
         yerr:=yerr-dist;
         y:=y+incy;
       end;
     end;
     koords[dist][0]:=xend;
     koords[dist][1]:=yend;
end;					

1:

if(dx>=dy) then dist:=dx;					

1:

dist = (dx > dy)?dx:dy;					

1:

  dist := max(dx, dy);