/* ------------------------------------ */
/*  Save as :   we_det3.h                 */
/* ------------------------------------ */

/* ------------------------------------ */
double **inv_mZ(
double **A,
double **Inverse
)
{
nb_Z z = det_Z(A);
int r;

  if(!(z.r && z.i))
    {
     printf("\n inv_mZ error - A determinant is Zero");
     printf("\n Press return to continue.\n");
     fflush(stdout);
     getchar();
     exit(EXIT_FAILURE);
    }

  adjoint_mZ(A,Inverse);

  z = inv_Z(z);

  for(r=R1; r<A[R_SIZE][C0]; r++)

             mulR_mZ(Inverse,z,r);
             
return(Inverse);
}
/* ------------------------------------ */
/* ------------------------------------ */