двухмерный мерный массив, превратить в треугольник по гауссу - C++ - Ответ 2225483

#include <stdio.h>
#include <stdlib.h>
#include <iostream>
#include <time.h>
using namespace std;
int main(int argc, char *argv[])
{
    const int b=100;
    //Матрица системы
    int  A[b][b];
    double summa[b], A1[b][b], summa2[b];
    int t=0,n,m,i,j,k;
      cout<<"Vvtdite n: ";
        cin>>n;
         cout<<"Vvtdite m: ";
        cin>>m;
            for (i=0; i<n; i++){
for (j=0;j<m;j++) {
printf(" A[%i][%i]:",i+1,j+1);
scanf("%i",&A[i][j]);
}}
             for(i = 0; i < n; ++i) 
        {       summa[i]=0;
                summa2[i]=0;
                for(j = 0; j < m; ++j) 
                        cout<<A[i][j]<<"\t"; 
                cout<<endl;
        }
 
             for(i = 0; i < n; ++i) 
        {    for(j = 0; j < m; ++j) 
                 A1[i][j]=A[i][j];}
        //Число перестановок строк
    int rCnt=0;
    //Ведущий элемент
    double ldn;
    //Номер ведущего элемента
    int index;
        double temp;
 
    //Прямой ход
    for(j=0; j<m; ++j)
    {
        //Выбор ведущего элемента столбца
        ldn=A1[j][j];
        index=j;
        for(i=j; i<n; ++i)
            if(fabs(A1[i][j])>fabs(ldn))
            {
                ldn=A1[i][j];
                index=i;
            }
        //Если ведущий элемент - не на гл. диагонали
        //меняем строки местами
        if(index!=j)
        {
            ++rCnt;
            for(i=0; i<n; ++i)
            {
                temp=A1[j][i];
                A1[j][i]=A1[index][i];
                A1[index][i]=temp;
            }
        }
        
        for (i = 0; i<m; i++)
            A1[j][i] = A1[j][i]/ldn; //делим строку на ведущий элемент
 
        //Обнуление элементов ниже гл. диагонали
        for(i=j+1; i<n; ++i)
        {
               for(k=0; k<m; ++k)
                A1[i][k]= fabs (A1[i][k]-A1[j][k]*A1[i][k]);
        }   }
 
    for(i = 0; i < n; ++i)  
        {
                for(j = 0; j < m; ++j) 
                        cout<<A1[i][j]<<"\t"; 
                cout<<endl; 
        }
                for (i=0; i<n; i++)
        { for (j=0; j<m; j++)
        summa2[i]=summa2[i]+A1[i][j];
        //cout<<summa[i]<<endl;
}
 
                for (i=0; i<n; ++i)
                {  summa[i]=summa[i]/m;
                   summa2[i]=summa2[i]/m;
                    if (summa2[i]<summa[i]) 
                    {
                        t=t+1;
                    }
                }
  system("PAUSE");   //для просмотра результата
  return 0;
}