Linear Algebra and the C Language/a0g8

/* ------------------------------------ */
/*  Save as :   c00d.c                  */
/* ------------------------------------ */
#include      "v_a.h"
/* ------------------------------------ */
int main(void)
{
double  a[R3*C3]={ 1, 0,0,
                   9, 1,0,
                   9,-8,1};
                   	
double id[R3*C3]={ 1, 0,0,
                   0, 1,0,
                  -9, 8,1};                
              
double **A  = ca_A_mR(a,    i_mR(R3,C3));
double **ID = ca_A_mR(id,   i_mR(R3,C3));

double **AID = mul_mR(A,ID, i_mR(R3,C3));

  clrscrn();
  printf(" Put zero on the left of the pivot: [R3,C3]\n\n");
  
  printf(" ID:");
  p_mR(ID, S3,P0,C6);
  
  printf(" A:");
  p_mR(A, S3,P0,C6);
  
  printf(" A ID:");
  p_mR(AID, S3,P1,C6);
  stop();
   
  f_mR(A);
  f_mR(ID);
  f_mR(AID);

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

                                                                                       
 Put zero on the left of the pivot: [R3,C3]

 ID:
 +1  +0  +0 
 +0  +1  +0 
 -9  +8  +1 

 A:
 +1  +0  +0 
 +9  +1  +0 
 +9  -8  +1 

 A ID:
+1.0 +0.0 +0.0 
+9.0 +1.0 +0.0 
+0.0 +0.0 +1.0 

 Press return to continue.