/* ------------------------------------ */
/*  Save as :   gjfree_a.c                 */
/* ------------------------------------ */
#include "v_a.h"
/* ------------------------------------ */
int main(void)
{
double ab[4*7]={
 1,  3, -2,  0,  2,  0,  0,
 2,  6, -5, -2,  4, -3, -1,
 0,  0,  5, 10,  0, 15,  5,
 2,  6,  0,  8,  4, 18,  6,
};

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

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

  clrscrn();
  printf(" A :");
  p_mR(A,5,0,9);
  printf(" b :");
  p_mR(b,5,0,9);
  stop();

  clrscrn();
  printf(" Ab ");
  p_mR(Ab,5,0,9);
  gj_freevariable_mR(Ab,1,4);
  stop();

  f_mR(Ab);
  f_mR(b);
  f_mR(A);
  return 0;
}

/* ------------------------------------ */
/* ------------------------------------ */