#include<stdio.h>
#include<stdlib.h>
#include<time.h>
float terme_produit(int i,int j,int n,float a[n][n],float b[n][n])
{ float x=0;
  for(int k=0;k<n;k++) x+=a[i][k]*b[k][j];
  return x;
}
void produit_mat(int n,float a[n][n],float b[n][n],float c[n][n])
{ for(int i=0;i<n;i++) for(int j=0;j<n;j++)
   c[i][j]=terme_produit(i,j,n,a,b);
}
void initiale(int n,float c[n][n])
{ for(int i=0;i<n;i++) for(int j=0;j<n;j++)
   c[i][j]=rand()%201/10.-10.;
}
void aff_mat(int n,float a[n][n])
{ for(int i=0;printf("\n"),i<n;i++)
  for(int j=0;j<n;j++) printf("%7.2f",a[i][j]);
}
float det(float a[3][3])
{ 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[n][n], b[n][n], c[n][n];
  initiale(n,a), initiale(n,b), produit_mat(n,a,b,c);
  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);
  return 0;
}