Determinante de una matriz de n*n

Autor: Ricardo Alberto Rosales Guerre

Este programa obtiene el determinante de una matriz n*n empleando recursividad. Hay que darle como entreada el tamaño de la matriz (n) y los n*n elementos de la matriz (m[i][j]). Devuelve el determinante de la matriz.

#include <iostream.h>
#include <alloc.h>
#include <math.h>

double det(double **m, int b){
	double determinante = 0, aux = 0;
	int c;
	if(b==2)
		return m[0][0]*m[1][1] - m[1][0]*m[0][1];
	else{
		for(int j=0; j<b; j++){
			double **menor = (double **)malloc(sizeof(double)*(b-1));
			for(int h=0; h<(b-1); h++) menor[h] = (double *)malloc(sizeof(double)*(b-1));
			for(int k=1; k<b; k++){
				c = 0;
				for(int l=0; l<b; l++){
					if(l!=j){
						menor[k-1][c] = m[k][l];
						c++;
					}
				}
			}
			aux = pow(-1, 2+j)*m[0][j]*det(menor, b-1);
			determinante += aux;
			for(int q = 0; q<(b-1); q++)
				free(menor[q]);
			free(menor);
		}
		return determinante;
	}
}

int main(){
	int n;  // n = numero de renglones = numero de columnas
	double **m=NULL;
	cout<<\"Elija el numero de renglones o columnas que tiene su matriz: \";
	cin>>n;

// Se crea la matriz de forma dinamica
	m = (double **)malloc(sizeof(double)*n);
	for(int i=0; i<n; i++) m[i]=(double *)malloc(sizeof(double)*n);

// Pide el valor de cada elemento de la matriz
	for(int y=0; y<n; y++)
		for(int w=0; w<n; w++){
			cout<<\"Introduzca el valor de matriz[\"<<y<<\"][\"<<w<<\"]: \";
			cin>>m[y][w];
		}

// Despliego el resultado de la funcion det
	cout<<\"nDeterminante: \"<<det(m, n);

// Libero la memoria utilizada por la matriz
	for(int r=0; r<n; r++) free(m[r]);
	free(m);

return 0;
}

Utilidades

	Leer comentarios (12)
	Escribir comentario
	Puntuación: (1 voto)
	Votar
	Recomendar este código
	Estadísticas

Patrocinados

Hospedaje web y servidores dedicados linux por Ferca Network