/* ------------------------------------ */
/*  Save as :   inv_r.c                 */
/* ------------------------------------ */
#include "v_a.h"
/* ------------------------------------ */
void fun(int r)
{
double **A   = i_mR(r,r);
double **Inv = i_mR(r,r);

  clrscrn();
  printf(" A\n");
  p_mR(r_mR(A,9999.),8,0,6);
  stop();

  clrscrn();
  printf(" invgj_mR(A,invA); gauss jordan (AId)\n");
  pE_mR(invgj_mR(A,Inv),12,4,6);
  printf(" inv_mR(A,Inv);  det,...,adjoint,... \n");
  pE_mR(inv_mR(A,Inv),12,4,6);
  stop();

  clrscrn();
  printf(" Copy/Past into the octave window.\n\n");
  p_octave_mR(A,"A",1,0);
  printf(" inv(A)\n");
  pE_mR(Inv,12,4,6);

  f_mR(Inv);
  f_mR(A);
}
/* ------------------------------------ */
int main(void)
{
time_t t;

  srand(time(&t));
  
do
{
   fun(rp_I(5)+C1);

} while(stop_w());

  return 0;
}
/* ------------------------------------ */
/* ------------------------------------ */
/* ------------------------------------ */
/* ------------------------------------ */