#include <stdio.h>
#include <stdlib.h>

double** allocate(int n)
{    double **m;
     int i;
     m=(double**)malloc(n*sizeof(double*));
     for (i=0; i<n; i++)
        m[i]=(double*)malloc(n*sizeof(double));
     return m;
}

void deallocate(double** m, int n)
{    int i;
     for (i=0; i<n; i++)
        free(m[i]);
     free(m);
}

void IspisiMatricu(double** matrica, int n)
{    int i,j;
     for (i=0; i<n; i++)
     {   for(j=0; j<n; j++)
           printf("%f\t",matrica[i][j]);
         printf("\n");
     }       
}

void UnesiMatricu(double*** matrica,int *pn, FILE *fp)
{    int i,j;
     
     fscanf(fp,"%d",pn);

     *matrica=allocate(*pn);
     for (i=0; i<*pn; i++)
        for(j=0; j<*pn; j++)
        {   
            fscanf(fp,"%f",&((*matrica)[i][j]));
        }
     printf("Matrica je uneta\n");
     IspisiMatricu(*matrica,*pn);
}


double determinanta(double** matrica, int n)
{    int i;
     double** podmatrica;
     double det=0;
	 int znak;

     if (n==1) return matrica[0][0];
     podmatrica=allocate(n-1);
     znak=1;
     for (i=0; i<n; i++)
     {  int vrsta,kolona;
        for (kolona=0; kolona<i; kolona++)
           for(vrsta=1; vrsta<n; vrsta++)
             podmatrica[vrsta-1][kolona]=matrica[vrsta][kolona];
        for (kolona=i+1; kolona<n; kolona++)
           for(vrsta=1; vrsta<n; vrsta++)
             podmatrica[vrsta-1][kolona-1]=matrica[vrsta][kolona];

        det+=znak*matrica[0][i]*determinanta(podmatrica,n-1);
        znak*=-1;    
     }
     deallocate(podmatrica,n-1);
     return det;
}

main(int argc, char **argv)
{ 
	double **matrica;
    int n;

   FILE *ul;

   if(argc<2)
      {
         fprintf(stderr,"Nedovojan broj argumenata na ulazu!!!\n");
         exit(1);
      }
      /* Otvaranje datoteke za citanje */
      if(!(ul=fopen(argv[1],"rt")))
      {
         fprintf(stderr,"Ulazna datoteka ne postoji\n");
         exit(2);
      }

   UnesiMatricu(&matrica,&n,ul);
   printf("Determinanta je : %d\n",determinanta(matrica,n));
   if(!determinanta(matrica,n))	   
	   printf("Matrica je singularna\n");
   
   return 0;
}