Mathc matrices/a157

Installer et compiler ces fichiers dans votre répertoire de travail.

	`c00b.c`

/* ------------------------------------ */
/*  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 **EVector    =                 i_mR(RCA,RCA);
double **EVectorT   =                 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,EVector);
   transpose_mR(EVector,EVectorT);
    
/* EValue = EVectorT A EVector        */   
   mul_mR(EVectorT,A,T);
   mul_mR(T,EVector,EValue);  
    
   printf(" LogA:");
   f_eigs_mR(log,EValue, Log_EValue); 

/*  LogA = EVector Log_EValue EVectorT */        
      mul_mR(EVector,Log_EValue,T);
      mul_mR(T,EVectorT,LogA);
        p_mR(LogA, S9,P6, C6);
       
   f_mR(A);
   f_mR(EValue);   
       
   f_mR(Log_EValue);
   f_mR(LogA); 
      
   f_mR(EVector);
   f_mR(EVectorT);

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

  srand(time(&t));

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

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

Exemple de sortie écran :

 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