/*matin1 matrix inversion routine used throughout the code*/
void matin1(double avec[], 
  long int NDECLARE, 
  long int n1, 
  long int *nerror);