Mathc matrices/05h

Application

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

	c00a.c

/* ------------------------------------ */
/*  Save as :  c00a.c                   */
/* ------------------------------------ */
#include "v_a.h"
/* ------------------------------------ */
#define   RCA C5
/* ------------------------------------ */
int main(void)
{
double a[RCA*RCA]={
+0.329527604306, -0.365886364507, -0.067047239569, +0.001915635416, -0.287345312440, 
-0.365886364507, +0.800330583940, -0.036588636451, +0.001045389613, -0.156808441932, 
-0.067047239569, -0.036588636451, +0.993295276043, +0.000191563542, -0.028734531244, 
+0.001915635416, +0.001045389613, +0.000191563542, +0.999994526756, +0.000820986607, 
-0.287345312440, -0.156808441932, -0.028734531244, +0.000820986607, +0.876852008954                    
};

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

  clrscrn();
  printf(" Copy/Past into the octave windows \n\n");
  p_Octave_mR(A,"a",P12);
  printf(" [V, E] = eigs (a,%d) \n\n",RCA);
   
  eigs_mR(A,EigsValue);
  printf(" EigsValue :");
  p_mR(EigsValue,S10,P4,C1);
     
  stop();
   
  f_mR(A);
  f_mR(EigsValue);  
     
    return 0;
}
/* ------------------------------------ */
/* ------------------------------------ */

Vérifier avec Octave. 

Exemple de sortie écran :

 Copy/Past into the octave windows 

 a=[
+0.329527604306,-0.365886364507,-0.067047239569,+0.001915635416,-0.28734531244;,
-0.365886364507,+0.800330583940,-0.036588636451,+0.001045389613,-0.15680844193;,
-0.067047239569,-0.036588636451,+0.993295276043,+0.000191563542,-0.02873453124;,
+0.001915635416,+0.001045389613,+0.000191563542,+0.999994526756,+0.00082098660;,
-0.287345312440,-0.156808441932,-0.028734531244,+0.000820986607,+0.87685200895],

 [V, E] = eigs (a,5) 

 EigsValue :
   +1.0000 
   +1.0000 
   +1.0000 
   +1.0000 
   -0.0000 

 Press return to continue.