void line(int x0, int y0, int x1, int y1)     /* first octant: 0 <= dy <= dx */
{
    int dx = x1 - x0, dy = y1 - y0;
    int e = 2 * dy - dx;
    int x, y = y0;

    for (x = x0; x <= x1; x++) {
        plot(x, y);
        if (e > 0) {
            y++;
            e -= 2 * dx;
        }
        e += 2 * dy;
    }
}