Skip to content

Instantly share code, notes, and snippets.

@bencz
Created March 17, 2013 09:13
Show Gist options
  • Save bencz/5180766 to your computer and use it in GitHub Desktop.
Save bencz/5180766 to your computer and use it in GitHub Desktop.
Programa para calcular o determinante de uma matriz NxN
#include <stdio.h>
#include <stdlib.h>
int main(void)
{
int m = 0;
double **a = 0;
int i = 0, j = 0, k = 0;
double factor = 0;
double temp = 0;
int count = 0;
printf("dimensao => ");
scanf("%d", &m);
a = (double **) calloc(m, sizeof(double *));
for(i = 0; i < m; i++)
{
a[i] = (double *) calloc(m, sizeof(double));
}
printf("\n\nEntre com o conteudo da matriz\n\n");
for(i = 0; i < m; i++)
{
for(j = 0; j < m; j++)
{
printf("A[%d ; %d] => ", i+1, j+1);
scanf("%lf", &a[i][j]);
}
}
// mostra a matriz
printf("\nMatriz digitada:\n");
for(i = 0; i < m; i++)
{
for(j = 0; j < m; j++)
{
printf("%8.3f ", a[i][j]);
}
printf("\n");
}
// faz a transformação em um triangulo...
for(i = 0; i < m - 1; i++)
{
if(a[i][i] == 0)
{
for(k = i; k < m; k++)
{
if(a[k][i] != 0)
{
for(j = 0; j < m; j++)
{
temp = a[i][j];
a[i][j] = a[k][j];
a[k][j] = temp;
}
k = m;
}
}
count++;
}
if(a[i][i] != 0)
{
for(k = i + 1; k < m; k++)
{
factor = -1.0 * a[k][i] / a[i][i];
for(j = i; j < m; j++)
{
a[k][j] = a[k][j] + (factor * a[i][j]);
}
}
}
}
temp = 1.0;
// Calcula o determinante
for(i = 0; i < m; i++)
temp *= a[i][i];
printf("\nDeterminante:\n");
if(count % 2 == 0)
printf("%8.3f \n", temp);
else
printf("%8.3f \n", -1.0 * temp);
return 0;
}
@jbrenorv
Copy link

What mathematical method was used to calculate the determinant?

@bencz
Copy link
Author

bencz commented Sep 29, 2019

Oh... a shot from the past! ahahha
This is the method: https://s-mat-pcs.oulu.fi/~mpa/matreng/eem3_4-2.htm

@jbrenorv
Copy link

Ok, does it work for any matrix?
From any matrix can one arrive at a triangular matrix?
And what is the definition of even higher order triangular matrix?

@bencz
Copy link
Author

bencz commented Sep 29, 2019

Vish, I do not remember ...

@jonesabr
Copy link

The allocated memory was not freed.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment