Trabajo en C++ agricola

La inversa de una matriz por método de Gauss-Jordan.

Un pequeño programa hecho en C, que utiliza este método. Antes que veamos el código veamos un pocode teoría...

Dada una matriz A, ¿Podremos encontrar otra matriz B tal que A·B=B·A=I?

Esta matriz B existe aunque no siempre, de existir se le llamamatriz inversa de A y se nota A-1. Para que exista la inversa de A, ésta tiene que ser cuadrada pues de lo contrario no se podría hacer el producto por laizquierda y por la derecha, luego cuando hablamos de matrices invertibles estamos hablando de matrices cuadradas.

Condición necesaria y suficiente para queuna matriz sea invertible es que no sea singular, es decir, que su determinante sea no nulo |A| ≠ 0

Cálculo de la matriz inversa

Método deGauss-Jordan

Este método consiste en colocar junto a la matriz de partida (A) la matriz identidad (I) y hacer operaciones por filas, afectando esas operacionestanto a A como a I, con el objeto de transformar la matriz A en la matriz identidad, la matriz resultante de las operaciones sobre I es la inversa de A (A-1).Las operaciones que podemos hacer sobre las filas son:

a) Sustituir una fila por ella multiplicada por una constante, por ejemplo, sustituimos la fila2 por ella multiplicada por 3.

b) Permutar dos filas

c) Sustituir una fila por una combinación lineal de ella y otras.

#include
#include#include
#include

/*********************************** imprime matriz **********************/
void desplegar(int n, double **A) {
int i,j,f,c;
for(i=1; i