Gauss v1.1, C++ (gcc)

Gauss v1.1

//g++  5.4.0

#include <iostream>
#include <iomanip>

using namespace std;

int gauss4(float [4][4], float [4], float [4]);
int mostrar_matriz(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);
    cout<<"\n";
    gauss4(matriz, b, sol);
    //
}

int gauss4(float matriz[4][4], float b[4], float sol[4]){
    float div[4];
    
    for(int i=0; i<4; i++) div[i]=matriz[i][0];
    
    for(int i=0; i<4; i++){
        for(int j=1; j<4; j++ ){
            matriz[j][i]=(matriz[j][i]*matriz[0][0])/div[j];
            b[j]=(b[j]*matriz[0][0])/div[j];
        }        
    }
    mostrar_matriz(matriz,b);
    
    for(int i=0; i<4; i++){
        for(int j=1; j<4; j++ ){
            matriz[j][i]=matriz[j][i]-matriz[0][i];
            b[j]=b[j]-b[0];
        }        
    }//printf("\n");
    
    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;
}


run \| edit \| history \| help	0