#include<stdio.h>
#include<stdlib.h>
#include<time.h>
float terme_produit(int i,int j,int n,float**a,float**b)
{ float x=0;
  for(int k=0;k<n;k++) x+=a[i][k]*b[k][j];
  return x;
}
float**produit_mat(int n,float**a,float**b)
{ if(!n) return 0;
  float**c=malloc(  n*sizeof(*c));
  float *t=malloc(n*n*sizeof(*t));
  for(int i=0;i<n;i++) c[i]=t+n*i;
  for(int i=0;i<n;i++) for(int j=0;j<n;j++)
   c[i][j]=terme_produit(i,j,n,a,b);
  return c;
}
float**initiale(int n)
{ if(!n) return 0;
  float**c=malloc(  n*sizeof(*c));
  float *t=malloc(n*n*sizeof(*t));
  for(int i=0;i<n;i++) c[i]=t+n*i;
  for(int i=0;i<n;i++) for(int j=0;j<n;j++)
   c[i][j]=rand()%201/10.-10.;
  return c;
}
void aff_mat(int n,float**a)
{ for(int i=0;printf("\n"),i<n;i++)
  for(int j=0;j<n;j++) printf("%7.2f",a[i][j]);
}
void libere_mat(float**a) {if(a) free(*a), free(a);}
float det(float**a)
{ return a[0][0]*a[1][1]*a[2][2]
        +a[0][1]*a[1][2]*a[2][0]
        +a[0][2]*a[1][0]*a[2][1]
        -a[0][0]*a[1][2]*a[2][1]
        -a[0][1]*a[1][0]*a[2][2]
        -a[0][2]*a[1][1]*a[2][0];
}
int main()
{ srand(time(0));
  int n=3;
  float **a=initiale(n), **b=initiale(n), **c=produit_mat(n,a,b);
  aff_mat(n,a);
  aff_mat(n,b);
  aff_mat(n,c);
  float da=det(a), db=det(b), dc=det(c), dd=da*db;
  printf("%f*%f=%f=%f diff=%e\n",da,db,dc,dd,dc-dd);
  libere_mat(a);
  libere_mat(b);
  libere_mat(c);
  return 0;
}