Linear Algebra and the C Language/a0io

Install and compile this file in your working directory.

/* ------------------------------------ */
/*  Save as :   c00c.c                  */
/* ------------------------------------ */
#include "v_a.h"
/* ------------------------------------ */
#define   RA R4
#define   CA C4
#define   Cb C1
/* ------------------------------------ */
int main(void)
{
double ab[RA*(CA+Cb)]={
    +964,     -689,     +474,     +948,     -159, 
    +640,     -693,     -937,     +656,     +718, 
    -981,     -337,     -287,     +420,     +812, 
    -807,     -123,     +404,     -177,     -105 
};
/* ------ Copy ab into Ab ------------- */
double **Ab = ca_A_mR(ab,   i_Abr_Ac_bc_mR(RA,CA,Cb));
/* ------------------------------------ */

/* ------ Copy A from Ab to A --------- */
double **A  = c_Ab_A_mR(Ab, i_mR(RA,CA));

/* ------ Copy b from Ab to b --------- */
double **b  = c_Ab_b_mR(Ab, i_mR(RA,Cb));

  clrscrn();
  printf(" A:");
  p_mR(A, S8,P0,C7);
  
  printf(" b:");
  p_mR(b, S8,P0,C7);
  
  printf(" Ab:");
  p_mR(Ab, S8,P0,C7);
  stop();

  clrscrn();
  printf(" Copy/Paste into the octave window.\n\n");
  p_Octave_mR(Ab,"Ab",P0);
  printf(" rref(Ab,.00000000001)\n\n\n\n");
  
  printf("  gj_TP_mR(Ab):");
  gj_TP_mR(Ab);
  p_mR(Ab, S10,P4,C7);
  stop();

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

Screen output example:

                                                                                       
 A:
    +964     -689     +474     +948 
    +640     -693     -937     +656 
    -981     -337     -287     +420 
    -807     -123     +404     -177 

 b:
    -159 
    +718 
    +812 
    -105 

 Ab:
    +964     -689     +474     +948     -159 
    +640     -693     -937     +656     +718 
    -981     -337     -287     +420     +812 
    -807     -123     +404     -177     -105 

 Press return to continue. 


 Copy/Paste into the octave window.

 Ab=[
+964,-689,+474,+948,-159;
+640,-693,-937,+656,+718;
-981,-337,-287,+420,+812;
-807,-123,+404,-177,-105]

 rref(Ab,.00000000001)



  gj_TP_mR(Ab):
   +1.0000    +0.0000    +0.0000    +0.0000    -0.3916 
   -0.0000    +1.0000    -0.0000    -0.0000    +0.1846 
   +0.0000    +0.0000    +1.0000    +0.0000    -0.6777 
   +0.0000    +0.0000    +0.0000    +1.0000    +0.7036 

 Press return to continue.