#include <stdio.h>
#include <limits.h>
#include <math.h>  

void eigvals(double *mat, double *evals, int n) ;
void eigvecs(double *mat, double *evals, double *evecs, int n) ;
void inverse(double *mat, int n);