/* ------------------------------------ */
/*  Save as :   gjfree.c                 */
/* ------------------------------------ */
#include "v_a.h"
/* ------------------------------------ */
void fun(int r)
{
double **A = r_mR(i_mR(r,r+C2),999.);
double **b = r_mR(i_mR(r,C1),999.);
double **Ab = c_A_b_Ab_mR(A,b,i_AbR0_mR(r,r+C2,C1));
/*                            i_AbR0_mR(R_Ab,C_A,C_b));  */

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

  clrscrn();
  gj_freevariable_mR(Ab,5,5);

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

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


} while(stop_w());

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