#include <stdio.h>
#include <stdlib.h>
typedef double R;
R det(R *, int, int);

R m[] = {1, 0, 0, 1};
R m2[] = {2, 3, 5, 4};
R m3[] = {1, 0, 0,  0, 1, 0,  0, 0, 1};
R m4[] = {0, 1, 0,  1, 0, 0,  0, 0, 1};
R m5[] = {
  3.3, 5.6, 4.8, 4.3,
  5.6, 4.2, 6.7, 5.2,
  7.6, 5.3, 6.6, 4.3,
  2.6, 4.8, 4.9, 9.6};
R m6[] = {0, 2, 3, 0, -2, 5, 0, 2, 2};
void Do(R * m, int j, R a) {
  printf("("/*)*/);
  {int k; for(k=0; k<j; ++k) {printf("("/*)*/);
    {int n; for(n=0; n<j; ++n) printf(" %4.2f", *(m+j*k+n));}
    printf(/*(*/")%s", k==j-1?"":"\n");}
  printf(/*(*/")\n");}
 {R z[j*j]; {int i=j; while(i--){int k=j; while(k--) z[j*k+i] = m[j*i+k];}}
  printf("D=%10.4f %10.4f %10.4f\n", det(m, j, j), det(z, j, j), a);}}

int main(){
Do(m, 2, 1); // 1
Do(m2, 2, -7); // -7
Do(m3, 3, 1); // 1
Do(m4, 3, -1); // -1
Do(m5, 4, 206.704); // 206.704
Do(m6, 3, 0); // 0
{R a[25], b[25]; {int j=25; while(j--) a[j] = b[j] = random();}
   a[0] = b[0] = 3;
 Do (a, 5, -6473428757628647e31);
 {int j = 5; while(j--) {R t = b[j]; b[j]=b[j+5]; b[j+5] = -t;}}
 Do (b, 5, -6473428757628647e31);}
return 0;}