Linear Algebra and the C Language/a0my

/* ------------------------------------ */
/*  Save as :  c02b.c                   */
/* ------------------------------------ */
#include "v_a.h"
/* ------------------------------------ */
/* ------------------------------------ */
#define   RA R2
#define   CA C2

#define   RX R2
#define   CX C1
/* ------------------------------------ */
/* ------------------------------------ */
int main(void)
{
double a[RA*CA]={
    +3, -2,  
    +2, -2
};

double x_S[RX*CX]={
  +10.00, 
  +13.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  -2.00 
 +2.00  -2.00 

 And x_s:
+10.00 
+13.00 

 T(x_s) = A x_s
  +4.00 
  -6.00 

 Press return to continue.