Construir un código computacional que permita realizar la resolución de ecuaciones simultáneas por Método de Gauss

#include <iostream>
#include <bits/stdc++.h>
using namespace std;

#define size 3

int gauss(double mat[][size+1])
{
    int i, j, k, c, flag = 1;
    double temp, pivote, numerador;

    for(i = 0; i < size; i++)
    {
        for(j = 0; j < size; j++)
        {
            if(j > i)
            {
                pivote = mat[j][i] / mat[i][i];
                for(k = 0; k < size + 1; k++)
                {
                    mat[j][k] -= pivote * mat[i][k];
                }
            }
        }
    }
    cout << "Aproximaciones: " << endl;
    for(j = size - 1; j >= 0; j--)
    {
        temp = 0;
        for(k = j; k < size; k++)
        {
            temp +=