Linear Algebra and the C Language/a0cd

/* ------------------------------------ */
/*  Save as :   c00b.c                  */
/* ------------------------------------ */
#include "v_a.h"
/* ------------------------------------ */
#define RCA      RC5
/* ------------------------------------ */
/* ------------------------------------ */
void fun(void)
{
double **A          = rEsymmetric_mR(i_mR(RCA,RCA),999,+1.E-3);
double **EValue     =                i_mR(RCA,RCA); 

double **Log_EValue = i_mR(RCA,RCA);
double **LogA       = i_mR(RCA,RCA);  

double **V          = i_mR(RCA,RCA);
double **V_T        = i_mR(RCA,RCA);

double **T          = i_mR(RCA,RCA);

   clrscrn();     
   printf(" Copy/Paste into the octave window.       \n\n");
   p_Octave_mR(A,"A", P6);
   printf(" logm (A)\n\n\n");
   
      eigs_V_mR(A,V);
   transpose_mR(V,V_T);
    
/* EValue = V_T * A * V        */   
   mul_mR(V_T,A,T);
   mul_mR(T,V,EValue);  
    
   printf(" LogA :");
   f_eigs_mR(log,EValue, Log_EValue);     
      mul_mR(V,Log_EValue,T);
      mul_mR(T,V_T,LogA);
        p_mR(    LogA, S9,P6, C6);
       
   f_mR(A);
   f_mR(EValue);   
       
   f_mR(Log_EValue);
   f_mR(LogA); 
      
   f_mR(V);
   f_mR(V_T);

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

  srand(time(&t));

do
{
    fun();
    
} while(stop_w());

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

                                                                                       
 Copy/Paste into the octave window.       

 A=[
+1.357734,-0.123337,-1.397583,-0.523404,+0.296125;
-0.123337,+0.350535,-0.197543,+0.337885,-0.144267;
-1.397583,-0.197543,+1.888885,+0.122765,+0.094678;
-0.523404,+0.337885,+0.122765,+0.642718,-0.602939;
+0.296125,-0.144267,+0.094678,-0.602939,+1.464218]

 logm (A)


 LogA :
-1.886935 -0.574739 -2.281263 -1.354664 +0.208755 
-0.574739 -1.888205 -0.716984 +0.981584 +0.179825 
-2.281263 -0.716984 -0.692397 -0.524562 +0.208617 
-1.354664 +0.981584 -0.524562 -1.870239 -0.793007 
+0.208755 +0.179825 +0.208617 -0.793007 +0.154821 


 Press   return to continue
 Press X return to stop