resolver matrices mostrando su procedimiento mediante el método de gauss jordan

//cpp
void gauss_jordan(float *mat, int n) {
	for(int i = 0; i < n; i++) {
		float pivote = mat[i*n + i];
		for(int j = 0; j < n; j++) {
			mat[i*n + j] /= pivote;
		}
		for(int k = 0; k < n; k++) {
			if(k != i) {
				float factor = mat[k*n + i];
				for(int j = 0; j < n; j++) {
					mat[k*n + j] -= factor * mat[i*n + j];
				}
			}
		}
	}
}

int main() {
	int n = 3;
	float *mat = new float[n*n];
	
	mat[0] = 1; mat[3] = 2; mat[6] = -1;