Linear Algebra and the C Language/a0aq


Install and compile this file in your working directory.

/* ------------------------------------ */
/*  Save as :   c00d.c                  */
/* ------------------------------------ */
#include "v_a.h"
/* ------------------------------------ */
void fun(int r)
{
double **u  =  i_mR(r,C1);

  clrscrn();
  printf(" u :");
  p_mR(u,S3,P0,C6);
  
  printf("      <u,u> >= 0     [Positivity axiom] \n\n");
  printf("      <u,u>  = 0  if and only if u = 0  \n\n");
  
  printf("      <u,u> = u^t.u   \n");
  printf("      <u,u> = %.0f  \n\n", dot_R(u,u));

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

  srand(time(&t));

do
{
  fun(rp_I(R3)+R2);


} while(stop_w());

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


Screen output example:

                                                                                       
 u :
 +0 
 +0 
 +0 
 +0 
 +0 

      <u,u> >= 0     [Positivity axiom] 

      <u,u>  = 0  if and only if u = 0  

      <u,u> = u^t.u   
      <u,u> = 0  


 Press   return to continue
 Press X return to stop