Usando el lenguaje C para realizar el método de eliminación gaussiano del elemento principal de la columna, encontré problemas
for(i=k+1;i { m[i-1][k]=a[i] [k]/a[k][k]; // Las coordenadas de m aquí deben ser [i-1][k] Si es [i][k], provocará fuera de límites for(j =0;j { temp=a[i][j]; a[i ][j]=temp -m[i-1][k]*a[k][j]; //Lo mismo aquí } } m es una matriz de 2X2 y a son datos de 3X3, es decir, la proporción de a[1][0] y a[0][0] debe almacenarse en m[0][0]. ¡Espero que esto ayude!