Estoy pensando en el algoritmo de triángulo de rasterización. (triangle_rasterization_lesson )

Escribí el siguiente 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);
        }
    }
}

Contriangle(100, 100, 100, 400, 400, 100, image, red); Funciona correctamente. Pero si cambio las coordenadas X (xa, ya) y Z (xc, yc) a no llena mi cuadrado.

Contriangle(70, 50, 200, 100, 20, 150, image, red); dibuja un triángulo, pero llenandosale de límites

¿Dónde está el problema?