Linear Algebra and the C Language/a0md


Install and compile this file in your working directory.

/* ------------------------------------ */
/*  Save as :   c00b.c                  */
/* ------------------------------------ */
#include "v_a.h"
/* ------------------------------------ */     
#define RCA          RC3  
#define EVALUE       RC3  
/* ------------------------------------ */       
/* ------------------------------------ */
void fun(void)
{                          
double a[RCA*RCA] ={   
+1.345288030, +1.760536858, -1.232763852, 
+1.760536858, +4.139646496, +0.479295438, 
-1.232763852, +0.479295438, +3.515065475    
};
                       
double **A      =   ca_A_mR(a,        i_mR(RCA,RCA));
double **V      = eigs_V_mR(A,        i_mR(RCA,RCA));

double **V3     =    c_c_mR(V,EVALUE, i_mR(RCA, C1),C1);
double **AV3    =    mul_mR(A,V3,     i_mR(RCA, C1)); 

  clrscrn(); 
  printf(" We know: A V3 = EVALUE3 V3             \n\n"
         " Verify if A V3 = 0 with A!= 0  and V3 != 0\n\n"
         " When  EVALUE3 = 0                         \n\n");
  stop();
         
  clrscrn(); 
  printf(" A:");
  p_mR(A, S9,P5, C3);     

  printf(" V:");
  p_mR(V, S9,P5, C4); 
  
  printf(" V%d:",EVALUE);
  p_mR(V3, S9,P5, C3);

  printf(" A V%d:",EVALUE);
  p_mR(AV3, S9,P5, C3);
  stop();
              
  f_mR(A);
  f_mR(V);  

  f_mR(V3); 
  f_mR(AV3); 
}
/* ------------------------------------ */
int main(void)
{
    fun();
    
  return 0;
}
/* ------------------------------------ */
/* ------------------------------------ */

Screen output example:

                                                                                       
 We know: A V3 = EVALUE3 V3             

 Verify if A V3 = 0 with A!= 0  and V3 != 0

 When  EVALUE3 = 0                         

 Press return to continue. 


 A:
 +1.34529  +1.76054  -1.23276 
 +1.76054  +4.13965  +0.47930 
 -1.23276  +0.47930  +3.51507 

 V:
 +0.45779  -0.27269  +0.84621 
 +0.88400  +0.24101  -0.40057 
 -0.09471  +0.93143  +0.35139 

 V3:
 +0.84621 
 -0.40057 
 +0.35139 

 A V3:
 +0.00000 
 -0.00000 
 +0.00000 

 Press return to continue.