Matriz Inversa usando determinante (nxn)

[hezelmet]

Saludos!

Es la segunda vez que escribo en este foro... la primera vez me sirvio de muchisima ayuda, la cual les agradezco infinitamente!!

Esta vez necesito realizar un programa que me entregue la matriz inversa de una matriz n x n...

No me exigen ningun metodo en especifico... decidi hacerlo por utilizando determinantes ya que se me hace mas sencillo (creo) que el metodo de Gauss-Jordan...

En fin, logre hacer el programa para una matriz 3 x 3... Necesito que, aunque sea... me den ideas de como puedo lograr que el programa funcione para cualquier valor de n... 4x4, 5x5, 6x6... nxn...

Aqui lo que tengo:

===============================================

#include<iostream.h>
#define n 3
#define filas 3
#define columnas 3


int main()

{ int a,b;
double B[n][n],A[filas][columnas],x;

for(a=0;a<=filas-1;a++)
{
for(b=0;b<=columnas-1;b++)

{

cout<<"A ["<<a<<"]["<<b<<"] = ";
cin>>A[a][b];

}
}

for(a=0;a<=filas-1;a++)

{
for(b=0;b<=columnas-1;b++)

{

x=(A[0][0]*((A[1][1]*A[2][2])-(A[2][1]*A[1][2])))-(A[0][1]*((A[1][0]*A[2][2])-(A[2][0]*A[1][2])))+(A[0][2]*((A[1][0]*A[2][1])-(A[2][0]*A[1][1])));
}}

cout<<"Determinate = "<<x<<endl; //x es el determinante

if (x!=0)
{

for(a=0;a<=filas-1;a++)
{
for(b=0;b<=columnas-1;b++)
{

B[0][0]=((A[1][1]*A[2][2])-(A[2][1]*A[1][2]))*(1/x);

B[0][1]=(-((A[0][1]*A[2][2])-(A[2][1]*A[0][2])))*(1/x);

B[0][2]=((A[0][1]*A[1][2])-(A[0][2]*A[1][1]))*(1/x);

B[1][0]=(-((A[1][0]*A[2][2]-A[2][0]*A[1][2])))*(1/x);

B[1][1]=((A[0][0]*A[2][2])-(A[2][0]*A[0][2]))*(1/x);

B[1][2]=(-((A[0][0]*A[1][2]-A[1][0]*A[0][2])))*(1/x);

B[2][0]=((A[1][0]*A[2][1])-(A[1][1]*A[2][0]))*(1/x);

B[2][1]=(-((A[0][0]*A[2][1]-A[2][0]*A[0][1])))*(1/x);

B[2][2]=((A[0][0]*A[1][1])-(A[1][0]*A[0][1]))*(1/x); cout<<"A^-1 ["<<a<<"]["<<b<<"] = "<<B[a][b]<<endl;
}
}
}
else
{
cout<<"Matriz no tiene inversa"<<endl;
}
return 0;
}

===============================================

Cualquier ayuda sera agradecida!!

[hezelmet]

Investigue y me acabo de dar cuenta que para matrices mayores a 3x3 hacerlo por Determinante es casi imposible :/


Habra algun otro metodo para relaizar este ejercicio?

[digies]

En C estándar y con el método de Gauss-Jordan lo podrías hacer así:

[hezelmet]

Muchisimas gracias digies!
Ya lo converti a C++

Ahora solo tengo que averiguar como determinar si es o no inversible la matriz para que el programa asi me lo indique...

Saludos y gracias nuevamente!