Aller au contenu

Mathc matrices/04h

Un livre de Wikilivres.


É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 R4
/* ------------------------------------ */
int main(void)
{
double a[RCA*RCA]={
+9.825339522090, +7.309094034726, -3.954100051573, -1.198212136840, 
+7.309094034726, +5.437252879491, -2.941464672512, -0.891352931064, 
-3.954100051573, -2.941464672512, +1.591284167096, +0.482207323363, 
-1.198212136840, -0.891352931064, +0.482207323363, +0.146123431322   
};

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

  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);
  stop();

  f_mR(A);
  
  return 0;
}
/* ------------------------------------ */
/* ------------------------------------ */


Vérifier avec octave


Exemple de sortie écran :
 Copy/Past into the octave windows 

 a=[
+9.825339522090,+7.309094034726,-3.954100051573,-1.198212136840;
+7.309094034726,+5.437252879491,-2.941464672512,-0.891352931064;
-3.954100051573,-2.941464672512,+1.591284167096,+0.482207323363;
-1.198212136840,-0.891352931064,+0.482207323363,+0.146123431322]

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

 Press return to continue.