En el siguiente link:

http://blogel.espaciolinux.info/wp-content/archivos/zeka/GAUSSJORDAN.zip

tenemos un programa escrito en C++ (*.zip) que resuelve sistemas lineales utilizando el método de Gauss-Jordan (subrutina gaussjordan.h). Se emplea en el siguiente sistema linealmente independiente (el programa no verifica dependencia lineal):

4x + 2y + 3z -v + w = 16
x + y + z + v + w = 3
x -y + z -v + w = 1
x + y + 2z + v + w = 6
x + y -z + v + 2w = -5

El número de incógnitas y la matriz extendida del sistema los introduce desde datos/datos.dat y los resultados los imprime, con formato, en pantalla y en el archivo salida/salida.dat. Después de desempaquetado el programa, nos ubicamos a través de una cónsola en la carpeta GAUSSJORDAN y compilamos con:

g++ gaussjordan.c++ -o gaussjordan

Su ejecución (./gaussjordan [Enter]) produce está salida en pantalla:

y una similar en datos/datos.dat. Los valores de las incógnitas son 1, 2, 3, -1, -2 para x, y, z, v, w; respectivamente.