Mathc matrices/c11bn

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 c0bn.c
'
/* ------------------------------------ */
/*  Save as :  c0bn.c               */
/* ------------------------------------ */
#include "v_a.h"
/* ------------------------------------ */
void fun(int r)
{
double **A         = rsymmetric_mR(i_mR(r,r),999.);

double **EigsVector=               i_mR(r,r);
double **T_EigsVector=             i_mR(r,r);

double **EigsValue =               i_mR(r,C1);
double **eigs_values=              i_mR(r,r);

double **T1=                       i_mR(r,r);

  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,S10,P4,C6);

  printf(" EigsValue : T_EigsVector * A * EigsVector");  
  transpose_mR(EigsVector,T_EigsVector);   
  mul_mR(T_EigsVector,A,T1);
  mul_mR(T1,EigsVector,eigs_values);
  p_mR(eigs_values,S10,P4,C6);  
   
  eigs_value_uv_mR(A,EigsValue);
  printf(" EigsValue :");
  p_mR(EigsValue,S10,P4,C1);   
   
  f_mR(A);
  
  f_mR(EigsVector);
  f_mR(T_EigsVector);
  
  f_mR(EigsValue);  
  f_mR(eigs_values);
       
  f_mR(T1);
}
/* ------------------------------------ */
int main(void)
{
time_t t;

  srand(time(&t));

do
{
  fun(rp_I(2)+R2);

} while(stop_w());

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


Déclaration des fichiers h.



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

 a=[
-833,-347,-373,-301;
-347,+294,-87,+462;
-373,-87,-749,-41;
-301,+462,-41,-197]

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

 Press return to continue. 


 EigsVector:
   +0.7671    -0.2853    +0.3425    -0.4614 
   +0.1531    +0.8114    -0.3057    -0.4742 
   +0.5984    +0.0103    -0.5607    +0.5722 
   +0.1735    +0.5101    +0.6891    +0.4846 

 EigsValue : T_EigsVector * A * EigsVector
-1261.3064    +0.0000    -0.0000    +0.0000 
   +0.0000  +705.3783    +0.0000    +0.0000 
   -0.0000    +0.0000  -518.2077    -0.0000 
   +0.0000    -0.0000    -0.0000  -410.8641 

 EigsValue :
-1261.3064 
 +705.3783 
 -518.2077 
 -410.8641 


 Press return to continue
 Press X      to stop