Linear Algebra and the C Language/a05p

/* ------------------------------------ */
/*  Save as :   c00a.c                  */
/* ------------------------------------ */
#include "v_a.h"
/* ------------------------------------ */
void fun(void)
{ 
double **U_T = r_mR(i_mR(R1, C3), 9);  
double **V_T = r_mR(i_mR(R1, C3), 9);   
double **W_T = r_mR(i_mR(R1, C3), 9);
 
double **UxV           = UxV_mR(U_T,V_T,    i_mR(R1,C3));
double **UxW           = UxV_mR(U_T,W_T,    i_mR(R1,C3));
double **UxV_pls_UxW   = add_mR(UxV,UxW,    i_mR(R1,C3));

double **VplusW        = add_mR(V_T,W_T,    i_mR(R1,C3));
double **Ux_VplusW     = UxV_mR(U_T,VplusW, i_mR(R1,C3));

  clrscrn();       
  printf("  u_T:");
  p_mR(U_T, S4,P0,C6);
  printf("  v_T:");
  p_mR(V_T, S4,P0,C6);
  printf("  w_T:");
  p_mR(W_T, S4,P0,C6);

  printf("\n\n\n");
  printf("    u  x (v+w):");    
  p_mR(Ux_VplusW, S5,P0,C6);                       
  printf(" (uxv) + (uxw):");
  p_mR(UxV_pls_UxW, S5,P0,C6);
  
  f_mR(U_T); 
  f_mR(V_T);
  f_mR(W_T);  

  f_mR(UxV); 
  f_mR(UxW);
  f_mR(UxV_pls_UxW); 

  f_mR(VplusW);
  f_mR(Ux_VplusW);  
}
/* ------------------------------------ */
int main(void)
{
time_t t;

  srand(time(&t));
  
do
{  
  fun();

} while(stop_w());

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

  u_T  :
  -2   +7   +8 

  v_T  :
  -5   +6   -5 

  w_T  :
  -6   +1   -9 




    u  x (v+w) :
 -154  -116   +63 

 (uxv) + (uxw) :
 -154  -116   +63 


 Press   return to continue
 Press X return to stop