Estou pensando no algoritmo do triângulo de rasterização. (triangle_rasterization_lesson )

Anotei o seguinte código:

void triangle(int xa, int ya, int xb, int yb, int xc, int yc, TGAImage &image, TGAColor color)
{
    line(xa, ya, xb, yb, image, color);
    line(xa, ya, xc, yc, image, color);
    line(xb, yb, xc, yc, image, color);
    for (int x = xa; x<=xb; x++)
    {
        for (int y = ya; y<=yb; y++)
        {
            line(xc, yc, x, y, image, white);
        }
    }
}

Comtriangle(100, 100, 100, 400, 400, 100, image, red); funciona corretamente. Mas se eu trocar as coordenadas X (xa, ya) e Z (xc, yc) para não preencher meu quadrado.

Comtriangle(70, 50, 200, 100, 20, 150, image, red); desenha triângulo, mas recheiosai dos limites.

Onde está o problema?