Mathc matrices/04l

Étude du code

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

	c00a.c

/* ------------------------------------ */
/*  Save as :  c00a.c                   */
/* ------------------------------------ */
#include "v_a.h"
/* ------------------------------------ */
#define   RCA R3
/* ------------------------------------ */
int main(void)
{
double a[RCA*RCA]={
+2/3., -1/3., -2/3., 
-1/3., +2/3., -2/3., 
-2/3., -2/3., -1/3.     
};

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(" EigenValues  = eigs (a,%d) \n\n",RCA);
  stop();

  clrscrn();   
  eigs_mR(A,EigsValue);
  printf(" ERROR : EigsValue :");
  p_mR(EigsValue,S13,P6,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.666666666667,-0.333333333333,-0.666666666667;
-0.333333333333,+0.666666666667,-0.666666666667;
-0.666666666667,-0.666666666667,-0.333333333333]

 EigenValues  = eigs (a,3) 

=================== Octave ====================
EigenValues =

   1
  -1
   1

=============== C Language ====================

 ERROR : EigsValue :
    +0.666667 
    +0.666667 
    -0.333333 

 Press return to continue.