/* Simple matrix manipulation routines              */
/* Jan 25 2003       Example 10                      */
#include<stdio.h>
#include<stdlib.h>
#include<sys/types.h>
#include<unistd.h>
#include<math.h>
#define SIZE 3
/* include a user-constructed header in working directory */
#include"matrix.h"


double a[SIZE][SIZE], c[SIZE][SIZE];
int n,m;
char label[10];
       
int main(){
    /* here is one way to initialize matrices */
    sprintf(label,"%s", "DATA");
    
    Mread(a,"DATA", 3);
/* c=a^{-1} */
Minv(a,c,12);

 for(n=0;n<SIZE;n++){
     for(m=0;m<SIZE;m++){
	 /* format the output fourplaces.fourplaces */
	 printf("%4.4f ",a[n][m]);
     }
     printf("\n");
 } 
 exit(0);

}