Programowanie C. wyznacznik macierzy

Question

Witam mam prablem z tym zadaniem

Napisac w jezyku C program obliczajacy wyznacznik macierzy 
o dowolnym rozmiarze. Macierz powinna byc wczytywana z pliku tekstowego 
o nastepujacej strukturze:
N
a11 a12 .. a1n
a21 a22 .. a2n
....................
an1 an2 .. ann

generalnie mam trochę kodu i generalny pomysł, ale prosiłbym o pomoc. jakieś wskazówki lub kod do analizy :)

#include <stdio.h>


double det(double ** matrix, int n)
{
  double** minor;
  double sum=0.0;
  double sign=1.0;
  int i, ...;
  if(n==1) return **matrix;//matrix[0][0]
  if(n==2) return matrix[0][0]*matrix[1][1]-matrix[0][1]*matrix[1][0];
  
  minor=matrix_alloc(n-1);
  for(i=0;i<n;i++)
  {
  /* skopiuj odp. podmacierz do minor */
  sum+=sign*matrix[0][i]*det(minor,n-1);
  sign = -sign;
  };
  matrix_free(minor);
  return sum;
};
int main()
{
	int n,i,j,k;
	double l;	
	FILE *f;
	f=fopen("matrix.txt","r");
	if(f==NULL)
	{
		printf("error");
		return 1;
	}
	
	fscanf(f,"%d",&n);
	

	for(i=0;i<n;i++)
	{
		for(j=0;j<n;j++)
		{
		fscanf(f,"%lf",&l);
		printf("%.2f ",l);	
		}if(n)printf("\n");
	}


	
fclose(f);

return 0;}

 

Answer 1

Co dokłądnie jest nie tak z tym kodem?

Answer 2

Po pierwsze treść pytania jest źle sformułowana. Nie można obliczyć wyznacznika macierzy o dowolnym rozmiarze. Można tylko dla kwadratowych. Nie wiem jak u Ciebie działają funkcje matrix_alloc i matrix_free, chociaż nietrudno się domyśleć(ale bez dokładnego kodu nie można mieć pewności). Nie masz żadnej funkcji, która tworzy minor na podstawie istniejącej macierzy i wyrazu odrzucanego. Poza tym próbowałbym ograniczyć metodę rozwinięć Laplace'a do n=4. Dla n=3 mamy tzw regułę Sarrusa. Chociaż program bez tego działałby poprawnie, ale musiałby wywoływać więcej razy tą funkcję.

I zdecydowanie w tych linikach jest coś nie tak

  /* skopiuj odp. podmacierz do minor */             //Rozumiem, że tworzysz ten minor o którym pisałem?
  sum+=sign*matrix[0][i]*det(minor,n-1);             //I temat minora się skończył, przesyłasz do funkcji jakąś niezainicjalizowaną tablicę
  sign = -sign;                                                           //I zmieniasz znak, w metodzie rozwinięć laplace'a jest (-1)^(n+m), gdzie n i m to numery wiersza i kolumny