Свести матрицу к треугольному виду

[Pein95]

void triangle_view(double mas[][N],int n)
{
	  

	for (int i = 1; i < N; ++i)
		for (int k = i; k < N; ++k)
			for (int j = N-1; j >= 0; --j)
				mas[k][j] -= ceil((mas[k][i-1] / mas[i-1][i-1] * mas[i-1][j]));
}

Етот алгоритм правильный?
просто иногда в матрице появляются результаты типа 1.#inf и -1.#ind

[Ezhik Kihze]

просто иногда в матрице появляются результаты типа 1.#inf и -1.#ind

Просто на нулевые элементы делить не надо...

[Pein95]

for (int i = 1; i < N; ++i)
		for (int k = i; k < N; ++k)
			for (int j = N-1; j >= 0; --j)
				if (mas[k][i-1] !=0 &&  mas[i-1][i-1] !=0  && mas[i-1][j] != 0)
				mas[k][j] -= ceil((mas[k][i-1] / mas[i-1][i-1] * mas[i-1][j]));

с такой поправкой все будет ок?