¿Cómo resolver Gauss jordán en (C)?

Estoy desarrollando un programa que me permite resolver un sistema de ecuaciones 3X3 por el método de Gauss Jordan utilizando lenguaje de programación (C). Estoy atorado en la parte de la solución en el código.

Las variables fueron guardadas en la mat[i][j].

Para resolver el sistema utilizo el siguiente método

(Los primeros 2 for)

Divido la filas entre la celda que al final tomara el valor de 1

2         4       6     |    8

2/2  4/2    6/2      8/2

Solucion:

1        2       3       |    4

Después necesito hacer el pivote de las columnas que siguen para convertir las celdas en 0

(Los segundos for) pero no hace lo hace

Les dejo el codigo para ver si me pueden ayudar a revisarlo

for(int a=1;a<=3;a++){
z=a;
pivote=mat[z][a];
for(int j=1;j<=(var+1);j++){

mat[a][j]=mat[a][j]/pivote;
}
for(int i=a+1;i<=3;i++){
pivote=mat[i][a];
for(int j=a;j<=3;j++){
res[i][j]=mat[i][j]-pivote*mat[a][j];
}
}