/* ------------------------------------ */
/*  Save as :   gj_r_a.c                 */
/* ------------------------------------ */
#include "v_a.h"
/* ------------------------------------ */
int main(void)
{
double ab[R4*C5]={
    +964,     -689,     +474,     +948,     -159, 
    +640,     -693,     -937,     +656,     +718, 
    -981,     -337,     -287,     +420,     +812, 
    -807,     -123,     +404,     -177,     -105 
};

/*                       i_AbR0_mR(R_Ab,C_A,C_b));  */

double **Ab = ca_A_mR(ab,i_AbR0_mR(R4,C4,C1));
double **A  = c_Ab_A_mR(Ab,i_mR(R4,C4));
double **b  = c_Ab_b_mR(Ab,i_mR(R4,C1));

  clrscrn();
  printf(" A :");
  p_mR(A,8,0,7);
  printf(" b :");
  p_mR(b,8,0,7);
  printf(" Ab :");
  p_mR(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);
  stop();

  f_mR(Ab);
  f_mR(b);
  f_mR(A);
  
  return 0;
}
/* ------------------------------------ */
/* ------------------------------------ */