Saludos a todos, he buscado intensamente un algoritmo que me convierta a matriz identidad un sistema de ecuaciones, buscando solo he encontrado el algoritmos de Gaus Seidel pero eso solo convierta a ceros por debajo de la diagonal y por arriba de la misma, lo que necesito es algo parecido a esto:
SISTEMA A RESOLVER
-4 -5 -1 0 = -51
2 -3 2 0 = -1
4 -5 -1 1 = -28
-3 -5 1 1 = -26
AL FINAL DE VARIAS ITERACIONES ME DEBE SALIR:
1 0 0 0 = 2
0 1 0 0 = 7
0 0 1 0 = 8
0 0 0 1 = 7
Lo debo programar en cualquier lenguaje de programacion (C,C++,Java,etc,etc).
Gracias de antemano, espero su pronta respuesta.