/* ------------------------------------ */
/*  Save as :   gj_z.c                 */
/* ------------------------------------ */
#include "w_a.h"
/* ------------------------------------ */
void fun(int r)
{
double **A  =      i_mZ(r,r);
double **b  =      i_mZ(r,  C1);
double **Ab = i_AbR0_mZ(r,r,C1);
/*            i_AbR0_mZ(R_Ab,C_A,C_b));*/

 do
  {
   r_mZ(A,999);
   printf(".");
  }while(!det_Z(A).r);

  clrscrn();
  printf(" A :");
  p_mZ(A,5,0,5,0,5);
  printf(" b :");
  p_mZ(r_mZ(b,999.),5,0,5,0,5);
  stop();

  clrscrn();
  printf(" Copy/Past into the octave window.\n\n");
  p_octave_mZ(c_A_b_Ab_mZ(A,b,Ab),"Ab",1,0, 0);
  printf(" rref(Ab,.00000000001)\n\n");
  printf(" gj_mZ(Ab) :");
  p_mZ(gj_mZ(Ab),1,5,1,5,3);

  f_mZ(Ab);
  f_mZ(b);
  f_mZ(A);
}
/* ------------------------------------ */
int main(void)
{
time_t t;

  srand(time(&t));

do
{
  fun(rp_I(4)+1);

} while(stop_w());

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