Gauss op, C++ (gcc)

Gauss op

//g++  5.4.0

#include <iostream>
#include <iomanip>

using namespace std;

int gauss4x4(float [4][4], float [4], float [4]);
int mostrar_matriz(float [4][4], float [4]);

int gauss3x3(float [4][4], float [4], float [4]);
int mostrar_matriz3(float [4][4], float [4]);

int main()
{
    float matriz[4][4];
    float sol[4];
    float b[4]={1,2,3,4};
    
    
    for(int i=0; i<4;i++){
        for(int j=0; j<4; j++){
            cin>>matriz[i][j];
        }
    }
    mostrar_matriz(matriz,b);
    mostrar_matriz3(matriz,b);
    cout<<"\n";
    //gauss4x4(matriz, b, sol);
    gauss3x3(matriz, b, sol);
    //
}

int gauss4x4(float matriz[4][4], float b[4], float sol[4]){
    float div[4];
    
    for(int k=0;k<4;k++){
        for(int i=0; i<4; i++) div[i]=matriz[i][k];
        for(int i=0; i<4; i++) matriz[k][i]=matriz[k][i]/div[k];
        b[k]=b[k]/div[k];

for(int i=k+1; i<4; i++){
            for(int j=k; j<4; j++){
                matriz[i][j]=matriz[i][j]-(div[i]*matriz[k][j])/matriz[k][k];
            }
            b[i]=b[i]-(div[i]*b[k])/matriz[k][k];
        }

}
    
    for(int i=0; i<4; i++) div[i]=matriz[i][3];
    for(int i=0; i<4; i++) matriz[3][i]=matriz[3][i]/div[3];
    b[3]=b[3]/div[3];
    
    
    //+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
    
    for(int k=3; k>0; k--){
        for(int i=0; i<4; i++) div[i]=matriz[i][k];

for(int i=0; i<k; i++){
            matriz[i][k]=matriz[i][k]-matriz[k][k]*matriz[i][k];
            b[i]=b[i]-b[k]*div[i];    
        }
    }
    mostrar_matriz(matriz,b);
    
    return 0;
}

int mostrar_matriz(float matriz[4][4],float b[4]){
    printf("\n");
    for(int i=0; i<4;i++){
        for(int j=0; j<4; j++){
            cout<<matriz[i][j]<<" ";
        }
        cout<<" "<<b[i];
        printf("\n");
    }
    return 0;
}

int gauss3x3(float matriz[4][4], float b[4], float sol[4]){
    float div[4];
    
    for(int k=0;k<3;k++){
        for(int i=0; i<3; i++) div[i]=matriz[i][k];
        for(int i=0; i<3; i++) matriz[k][i]=matriz[k][i]/div[k];
        b[k]=b[k]/div[k];

for(int i=k+1; i<3; i++){
            for(int j=k; j<3; j++){
                matriz[i][j]=matriz[i][j]-(div[i]*matriz[k][j])/matriz[k][k];
            }
            b[i]=b[i]-(div[i]*b[k])/matriz[k][k];
        }

}
    
    for(int i=0; i<3; i++) div[i]=matriz[i][2];
    for(int i=0; i<3; i++) matriz[2][i]=matriz[2][i]/div[2];
    b[2]=b[2]/div[2];
    
    
    //+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
    
    for(int k=2; k>0; k--){
        for(int i=0; i<3; i++) div[i]=matriz[i][k];

for(int i=0; i<k; i++){
            matriz[i][k]=matriz[i][k]-matriz[k][k]*matriz[i][k];
            b[i]=b[i]-b[k]*div[i];    
        }
    }
    mostrar_matriz3(matriz,b);
    
    return 0;
}

int mostrar_matriz3(float matriz[4][4],float b[4]){
    printf("\n");
    for(int i=0; i<3;i++){
        for(int j=0; j<3; j++){
            cout<<matriz[i][j]<<" ";
        }
        cout<<" "<<b[i];
        printf("\n");
    }
    return 0;
}


run \| edit \| history \| help	0