Mathc matrices/04c

Application

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

	c00b.c

/* ------------------------------------ */
/*  Save as :  c00b.c                   */
/* ------------------------------------ */
#include "v_a.h"
/* ------------------------------------ */
#define   RCA    RC6
/* ------------------------------------ */
int main(void)
{
double a[RCA*RCA]={
+0.853214156635, +0.223697489745, -0.130729221324, -0.180799128390, -0.123227848855, -0.101171318892, 
+0.223697489745, +0.659091328080, +0.199227650133, +0.275532777834, +0.187795769834, +0.154182239592, 
-0.130729221324, +0.199227650133, +0.883570997610, -0.161021858298, -0.109748190676, -0.090104382248, 
-0.180799128390, +0.275532777834, -0.161021858298, +0.777306012096, -0.151782264254, -0.124614784740, 
-0.123227848855, +0.187795769834, -0.109748190676, -0.151782264254, +0.896549269431, -0.084934103365, 
-0.101171318892, +0.154182239592, -0.090104382248, -0.124614784740, -0.084934103365, +0.930268236148};

double **A      = ca_A_mR(a, i_mR(RCA,RCA));
double **EValue = eigs_mR(A, i_mR(RCA,C1));

  clrscrn();
  printf(" Copy/Paste into the octave windows \n\n");
  p_Octave_mR(A,"a",P12);
  printf(" EigenValues  = eigs (a,%d) \n\n",RCA);
  stop();
  
  clrscrn();    
  printf(" EigenValues:");
  p_mR(EValue, S13,P3,C1); 
  stop();  
   
  f_mR(A);
  f_mR(EValue);
  
  return 0;
}
/* ------------------------------------ */
/* ------------------------------------ */

Exemple de sortie écran :

 Copy/Paste into the octave windows 

 a=[
+0.853214156635,+0.223697489745,-0.130729221324,-0.180799128390,-0.123227848855,-0.101171318892;
+0.223697489745,+0.659091328080,+0.199227650133,+0.275532777834,+0.187795769834,+0.154182239592;
-0.130729221324,+0.199227650133,+0.883570997610,-0.161021858298,-0.109748190676,-0.090104382248;
-0.180799128390,+0.275532777834,-0.161021858298,+0.777306012096,-0.151782264254,-0.124614784740;
-0.123227848855,+0.187795769834,-0.109748190676,-0.151782264254,+0.896549269431,-0.084934103365;
-0.101171318892,+0.154182239592,-0.090104382248,-0.124614784740,-0.084934103365,+0.930268236148]

 EigenValues  = eigs (a,6) 

 Press return to continue. 


 EigenValues:
       +1.000 
       +1.000 
       +1.000 
       +1.000 
       +1.000 
       -0.000 

 Press return to continue.