Mathc matrices/a86

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

	`c00a.c`

/* ------------------------------------ */
/*  Save as:  c00a.c                    */
/* ------------------------------------ */
#include "v_a.h"
/* ------------------------------------ */
#define   RCA        RC5
/* ------------------------------------ */
/* ------------------------------------ */
void fun(void)
{
double **A            = rEsymmetric_mR( i_mR(RCA,RCA),555.,1E-3);
double **coshA        =                 i_mR(RCA,RCA);
double **coshA_EValue =                 i_mR(RCA,RCA);

double **EVector      =                 i_mR(RCA,RCA);
double **EVectorT     =                 i_mR(RCA,RCA);

double **EValue       =                 i_mR(RCA,RCA);
double **cosh_EValue  =                 i_mR(RCA,RCA);
 
double **T1           =                 i_mR(RCA,RCA);
  
  clrscrn();
  printf(" A:");
  p_mR(A,S10,P4,C6); 

/* ================================================================== */ 
  eigs_V_mR(A,EVector); 
  transpose_mR(EVector,EVectorT);
      
/* EValue = EVectorT A EVector */
  mul_mR(EVectorT,A,T1);
  mul_mR(T1,EVector,EValue);     
  
  printf(" cosh(A):");
  f_eigs_mR(cosh,EValue,cosh_EValue); 

/*  coshA = EVector cosh_EValue EVectorT */    
     mul_mR(EVector,cosh_EValue,T1);
     mul_mR(T1,EVectorT,coshA);
     p_mR(coshA,S10,P4,C6); 

/* ================================================================== */ 
  eigs_V_mR(coshA,EVector); 
  transpose_mR(EVector,EVectorT);  
  
/* coshA_EValue = EVectorT coshA EVector */
  mul_mR(EVectorT,coshA,T1);
  mul_mR(T1,EVector,coshA_EValue);  

  printf(" acosh(cosh(A)) = A:");
  f_eigs_mR(acosh,coshA_EValue,EValue);

/*  A = EVector EValue EVectorT */    
     mul_mR(EVector,EValue,T1);
     mul_mR(T1,EVectorT,A);  
       p_mR(A,S10,P4,C6);

  f_mR(A);
  f_mR(coshA);
  f_mR(coshA_EValue);
  
  f_mR(EVector);
  f_mR(EVectorT);
  
  f_mR(EValue);  
  f_mR(cosh_EValue);    
      
  f_mR(T1);
}
/* ------------------------------------ */
int main(void)
{
time_t t;

  srand(time(&t));

do
{
  fun();

} while(stop_w());

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

Exemple de sortie écran :

 A :
   +0.0522    -0.0059    +0.0828    +0.0150    -0.0654 
   -0.0059    +0.5696    -0.0106    -0.2065    -0.2950 
   +0.0828    -0.0106    +0.2041    +0.0014    -0.2069 
   +0.0150    -0.2065    +0.0014    +0.3984    -0.1084 
   -0.0654    -0.2950    -0.2069    -0.1084    +0.6454 

 cosh(A):
   +1.0073    +0.0069    +0.0181    +0.0078    -0.0333 
   +0.0069    +1.2420    +0.0296    -0.0880    -0.1813 
   +0.0181    +0.0296    +1.0476    +0.0137    -0.0950 
   +0.0078    -0.0880    +0.0137    +1.1100    -0.0264 
   -0.0333    -0.1813    -0.0950    -0.0264    +1.3012 

 acosh(cosh(A)) = A :
   +0.0522    -0.0059    +0.0828    +0.0150    -0.0654 
   -0.0059    +0.5696    -0.0106    -0.2065    -0.2950 
   +0.0828    -0.0106    +0.2041    +0.0014    -0.2069 
   +0.0150    -0.2065    +0.0014    +0.3984    -0.1084 
   -0.0654    -0.2950    -0.2069    -0.1084    +0.6454 


 Press   return to continue
 Press X return to stop