Linear Algebra and the C Language/a033

Install and compile this file in your working directory.

/* ------------------------------------ */
/*  Save as :   c00b.c                  */
/* ------------------------------------ */
#include "v_a.h"
/* ------------------------------------ */
int main(void)
{ 
double u_T[R1*C3] = { 4, 2, 5};                                             
double v_T[R1*C3] = { 3, 4, 1};

double **U_T = ca_A_mR(u_T,   i_mR(R1, C3));
double **V_T = ca_A_mR(v_T,   i_mR(R1, C3));   
 
double **UxV = UxV_mR(U_T,V_T,i_mR(R1,C3));
double **VxU = UxV_mR(V_T,U_T,i_mR(R1,C3));

  clrscrn();       
  printf("  u_T:");
  p_mR(U_T, S3,P0,C6);
  printf("  v_T:");
  p_mR(V_T, S3,P0,C6);
  
  printf("    u x v: ");    
  p_mR(UxV, S5,P0,C6); 
  printf("    v x u: ");    
  p_mR(VxU, S5,P0,C6);
  stop();
  
  f_mR(U_T); 
  f_mR(V_T);
  f_mR(UxV);  
  f_mR(VxU);   

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

u x v == - v x u:.

Screen output example:

  u_T  :
 +4  +2  +5 

  v_T  :
 +3  +4  +1 

    u x v : 
  -18   +11   +10 

    v x u : 
  +18   -11   -10 

 Press return to continue.