/* ------------------------------------ */
/*  Save as :   gj_z_a.c                 */
/* ------------------------------------ */
#include "w_a.h"
/* ------------------------------------ */
void fun(void)
{
double a[R3*(C3*C2)] ={ 1,2,  3,4,  5,6,
                        5,4,  1,3,  6,8,
                        7,2,  5,1,  1,1};

double b0[R3*(C1*C2)] ={ 1,4, 
                         2,5, 
                         3,6};   
                        
double **A = ca_A_mZ(a, i_mZ(R3,C3));                        
double **b = ca_A_mZ(b0,i_mZ(R3,C1));                                             
                                               
double **Ab = c_A_b_Ab_mZ(A,b,i_AbR0_mZ(R3,  C3, C1));	
/*                            i_AbR0_mZ(R_Ab,C_A,C_b));*/
	
  clrscrn();
  printf(" A :");
  p_mZ(A,S4,P0,S3,P0,C5);
  printf(" b :");
  p_mZ(b,S4,P0,S3,P0,C5);
  stop();

  clrscrn();
  printf(" Copy/Past into the octave window.\n\n");
  p_octave_mZ(Ab,"Ab",1,0, 0);
  printf(" rref(Ab,.00000000001)\n\n");
  printf(" gj_mZ(Ab) :");
  p_mZ(gj_mZ(Ab),S5,P5,S9,P5,C3);
  stop();

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

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