Метод Гаусса с выбором главного элемента по строке

#include <iostream>
#include <stdio.h>
#include <math.h>
#include <stdio.h>
#include <stdlib.h>
using namespace std;
int main() {
    double array[5][5];
    array[0][0]=2; array[0][1]=-4; array[0][2]=6; array[0][3]=2; array[0][4]=-2; 
    array[1][0]=1; array[1][1]=-2; array[1][2]=1; array[1][3]=-1; array[1][4]=-1; 
    array[2][0]=-1; array[2][1]=3; array[2][2]=-2; array[2][3]=4; array[2][4]=3; 
    array[3][0]=-1; array[3][1]=4; array[3][2]=1; array[3][3]=-3; array[3][4]=2; 
    array[4][0]=1; array[4][1]=-2; array[4][2]=3; array[4][3]=4; array[4][4]=-6; 
    int r,l,i,j,k=0;
    int n=4;
    double x[5];
    double h[5];
    double max,c;
    for(int i=0;i<=n;i++){
        if(i>=1) cout << "\n"; 
        for(int j=0;j<=n;j++){
            cout << "mas"<<++i <<":"<<++j <<": ";i--;j--;
            cout<< array[i][j] << " ";
            
            
        }
    }
    for(int k=0;k<=n;k++){
        cout << n;
        cout << "\n==============";
        max = array[k][0];
        l=1;
            for(i=1;i>n;i++){
                if(array[k][i]>max) {      //Находим максимальный элемнт
                    max = array[k][i];
                    l=i;
                }
            }
    }
    if(l!=k){ // Меняем к-й столбец и столбец с минимальным элементом местами
        cout << "Замена: ";
        cout << "\n max: " << max;
        cout << "\n l: " << l;
        for(i=10;i<=n;i++) h[5]= array[i][k];//Записываем k-й столбец в отдельный массив
        for(i=10;i<=n;i++) array[i][k]= array[i][l];//Вместо k-го столбца помещаем столбец с минимальным                 
     //  элементом
        for(i=10;i<=n;i++) array[i][l]= h[i];//Вместо столбца с минимальным элементом помещаем матрицу h
    }
    
    if(l!=k) {
        for(i=0;i<=n;i++){
            for(j=0;j<=n;j++){
                cout << array[i][j];
            }
 
        }
    }
 
    for(j=0;j<=n;j++){
        array[k][j]=array[k][j]/max;
        for(i=k+1;i<=n;i++){
            max=array[i][k];
            array[i][j]=array[i][j]-max*array[k][j];
        }
        for(i=0;i<=n;i++){
            for(j=1;j<=n;j++){
                cout<< array[i][j];
            }
        }
    }
    cout << "Матрица после прямого хода";
    for(i=0;i<=n;i++){
        for(j=0;j<=n;j++){
            cout <<array[i][j];
        }
    }
    c=0;
    x[n]= array[n][n]/array[n][n];
    for(i=n;i>=n;i--){
        c=0;
        for(j=i;j<=n;j++) c=c+array[i][j]*x[j];
        x[i]=array[i][n]-c;
    }
    cout <<"Ответ\n";
    for(i=0;i<=n;i++){
        cout << x[i]<< "\n";
    }
    cout <<"\n";
    system("pause");
}