Linear Algebra and the C Language/a0mp

/* ------------------------------------ */
/*  Save as :  c00b.c                   */
/* ------------------------------------ */
#include "v_a.h"
/* ------------------------------------ */
/* ------------------------------------ */
#define   RA R4
#define   CA C4

#define   RX R4
#define   CX C1
/* ------------------------------------ */
/* ------------------------------------ */
int main(void)
{
double a[RA*CA]={
    +3, -4, -2, -1,
    +5, -3, -0, -2,      
    +6, -0, -3, -2,
    +3, -3, -1, -3   
};

double x_S[RX*CX]={
  +47.00, 
  +34.00, 
  +48.00, 
  +37.00 
};

double **A    = ca_A_mR(a,     i_mR(RA,CA));
double **x_s  = ca_A_mR(x_S,   i_mR(RX,CX));
double **Ax_s =  mul_mR(A,x_s, i_mR(RA,CX));

  clrscrn();
  printf(" In the Standard basis\n\n"
         " Calculate the linear application T(x_s) = A x_s \n\n"
         " With A :");
  p_mR(A,S6,P2,C7);
  
  printf(" And x_s:");
  p_mR(x_s,S6,P2,C7);
  
  printf(" T(x_s) = A x_s");
  p_mR(Ax_s,S7,P2,C7);
  stop();


  f_mR(A);
  f_mR(x_s);
  f_mR(Ax_s);
  
  return 0;
}
/* ------------------------------------ */
/* ------------------------------------ */

                                                                                       
 In the Standard basis

 Calculate the linear application T(x_s) = A x_s 

 With A :
 +3.00  -4.00  -2.00  -1.00 
 +5.00  -3.00  +0.00  -2.00 
 +6.00  +0.00  -3.00  -2.00 
 +3.00  -3.00  -1.00  -3.00 

 And x_s:
+47.00 
+34.00 
+48.00 
+37.00 

 T(x_s) = A x_s
-128.00 
 +59.00 
 +64.00 
-120.00 

 Press return to continue.