/* ------------------------------------ */
/*  Save as :   gj1_rsing.c             */
/* ------------------------------------ */
#include "v_a.h"
/* ------------------------------------ */
void fun(int r)
{
double **Ab = r_singular_Ab_mR(i_AbR0_mR(r,r,C1),999);
double **A  = c_Ab_A_mR(Ab,i_mR(r,r));
double **b  = c_Ab_b_mR(Ab,i_mR(r,C1));


  clrscrn();
  printf(" A :");
  p_mR(A,8,0,7);
  printf(" b :");
  p_mR(r_mR(b,999.),8,0,7); /* I change the value of b */
  printf(" Ab :");
  p_mR(c_A_b_Ab_mR(A,b,Ab),8,0,7);
  getchar();

  clrscrn();
  printf(" Copy/Past into the octave window.\n\n");
  p_octave_mR(Ab,"Ab",1,0);
  printf("\n rref(Ab,.00000000001)\n\n");
  printf(" gj_mR(Ab) :");
  p_mR(gj_mR(Ab),10,5,7);

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

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


} while(stop_w());

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