C++. Двумерный массив

[Trinity13]

Помогите, пожалуйста, решить задачу:
Коэффициенты системы линейных уравнений заданы в виде прямоугольной матрицы. С помощью допустимых преобразований привести систему к треугольному виду. Найти количество строк, среднее арифметическое элементов которых меньше заданной величины.
Усть некоторые соображения по этому поводу, но прога не работает =(
#include <iostream>
#include <ctime>
using namespace std;

int NOD (int a, int b)
{
int c;
c=(a>b)?a:b;
while (!((c%b==0)&&(c%a==0)))
c++;
return c;
}

void main()
{
int n, i, j, k, l, c, sign;
cout<<"Write how many equations?-->";
cin>>n;
int **mas=new int *[n];
for (i=0;i<n;i++)
mas[i]=new int[n];
srand(time(NULL));
for (i=0;i<n;i++)
{
for (j=0;j<n;j++)
{
mas[i][j]=rand()%9+1;
sign=rand()%4;
if (sign==0)
mas[i][j]=mas[i][j]*(-1);
cout<<mas[i][j]<<'\t';
}
cout<<'\n';
}
for (j=n-1;j>0;j--)
{
for(i=n-1;i>0;i--)
{
c=NOD(mas[i][j],mas[i-1][j]);
for (k=0;k<n;k++)
{
mas[i][k]=mas[i][k]*(c%mas[i][k]);
mas[i-1][k]=mas[i-1][k]*(c%mas[i-1][k]);
mas[i][k]=mas[i][k]-mas[i-1][k];
}
}
}
for (i=0;i<n;i++)
{
for (j=0;j<n;j++)
cout<<mas[i][j]<<'\t';
cout<<'\n';
}
}
Заранее спасибо! =)

[Vago]

1. П/п следовало бы, всё-таки, назвать NOK

2. Фокусы с нахождением наименьшего общего кратного вообще-то имеют смысл только если заведомо известно, что система имеет целочисленное решение. Что на практике бывает редко. Но если уж очень хочется...

Верхняя треугольная матрица:

#!/usr/bin/python
# -*- coding: cp1251 -*-

from numpy import array

def NOK( a, b ):

    if a > b:
        c = a
    else:
        c = b

    while ( not ((c % b == 0) and (c % a == 0)) ):
        c += 1

    return c


mas = array( [[ 3, 3, 6], [ -2, 2, 5], [-4, 1, 7]], dtype = int )
print mas
n = len( mas )

for i in range ( 0, n-1 ):
    for j in range ( i+1, n ):
        c = NOK( mas[i][i], mas[j][i] )
        c1 = c / mas[i][i]
        c2 = c / mas[j][i]
        for k in range( i, n ):
            mas[j][k] = mas[i][k] * c1 - mas[j][k] * c2

print mas
#

3. Питоновский

for k in range( i, n )

эквивалентен C-шному

for ( k = i; k < n; k++ )

Остальные вычисления переводятся на C практически автоматически.

4. Понятно, что если a или b окажутся равными нулю, то п/п выдаст бред. Но эту ситуацию уж сами обходите...