/* ------------------------------------ */
/*  Save as :   ve_det3.h                 */
/* ------------------------------------ */
double **inv_mR(
double **A,
double **Inverse
)
{
double s = det_R(A);
int r;

  if(!s)
    {
     printf("\n inv_mR error - A determinant is Zero");
     printf("\n Press return to continue.\n");
     fflush(stdout);
     getchar();
     exit(EXIT_FAILURE);
    }

  adjoint_mR(A,Inverse);

  s = (1.)/s;

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

             mulR_mR(Inverse,s,r);
             
return(Inverse);
}
/* ------------------------------------ */
/* ------------------------------------ */