Mathc matrices/c11a3

Un livre de Wikilivres.
Aller à la navigation Aller à la recherche


Sommaire


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


Crystal Clear mimetype source c.png c0a3.c
'
/* ------------------------------------ */
/*  Save as :  c0a3.c               */
/* ------------------------------------ */
#include "v_a.h"
/* ------------------------------------ */
void fun(void)
{
double ab[R3*C3]={
+966,+98,+864,
+98,+692,+214,
+864,+214,+376
};

int r = RC3;
	
double **A         = ca_A_mR(ab,i_mR(r,r));
double **EigsVector=            i_mR(r,r);
double **EigsValue =            i_mR(r,C1);

  clrscrn();
  printf(" Copy/Past into the octave windows \n\n");
  p_Octave_mR(A,"a",P0);
  printf(" [V, E] = eigs (a,10) \n\n");
  stop();

  clrscrn();   
  eigs_V_uv_mR(A,EigsVector); 
  printf(" EigsVector:");
  p_mR(EigsVector,S5,P6,C6);
   
  eigs_value_uv_mR(A,EigsValue);
  printf(" EigsValue :");
  p_mR(EigsValue,S13,P6,C1);   
  stop();
     
  f_mR(A);
  f_mR(EigsVector);
  f_mR(EigsValue);
}
/* ------------------------------------ */
int main(void)
{

  fun();

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


Déclaration des fichiers h.



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

 a=[
+966,+98,+864;
+98,+692,+214;
+864,+214,+376]

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

 Press return to continue. 


/* ------------------------------------ */
 EignVector:
+0.786537 -0.247455 -0.565796 
+0.214554 +0.968632 -0.125376 
+0.579073 -0.022781 +0.814958 

 EigsValue :
 +1628.835762 
  +661.931131 
  -256.766893 

 Press return to continue.