Mathc complexes/04b

Application

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

	c00a.c

/* ------------------------------------ */
/*  Save as :   c00a.c                  */
/* ------------------------------------ */
#include "w_a.h"
/* ------------------------------------ */
#define RC        RC5
#define FACTOR_I  10
/* ------------------------------------ */
/* ------------------------------------ */
void fun(void)
{
double **A          = r_mRZ(i_mZ(RC,RC),999.);
double **EigsValue  =       i_mZ(RC,C1);

  clrscrn();
  printf(" Copy/Past into the octave windows \n\n");
  p_Octave_mZ(A,"a",P0,P0);  
  printf(" EigenValues  = eigs (a,%d)\n\n\n",RC);
  stop();

  clrscrn();
  printf(" EigsValue :");
  eigs_mZ(A,EigsValue);  
  p_mZ(EigsValue, S15,P9, S15,P9, C5); 
    
  eyeI_mZ(A,FACTOR_I);
         
  printf(" EigsValue(FACTOR_I) :");
  eigs_mZ(A,EigsValue);  
  p_mZ(EigsValue, S15,P9, S15,P9, C5); 
       
  f_mZ(A);
  f_mZ(EigsValue);
}
/* ------------------------------------ */
int main(void)
{
time_t t;

  srand(time(&t));

do
{
  fun();
  
} while(stop_w());

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

Exemple de sortie écran : Trouvons les valeurs complexes

 Copy/Past into the octave windows 

 a=[
+661+0*i,-574+0*i,+41+0*i,-245+0*i,-936+0*i;
-170+0*i,+165+0*i,-81+0*i,+175+0*i,-739+0*i;
-101+0*i,-773+0*i,-804+0*i,+742+0*i,-955+0*i;
-143+0*i,-336+0*i,+853+0*i,-948+0*i,-902+0*i;
-444+0*i,-10+0*i,+370+0*i,+86+0*i,-722+0*i]

 EigenValues  = eigs (a,5)


 Press return to continue. 


 EigsValue :
-1678.161191722   +0.000000000i 
 +859.134526218   +0.000000000i 
 -861.294778398   +0.000000000i 
 -345.492436810   +0.000000000i 
 +377.813880713   +0.000000000i 

 EigsValue(FACTOR_I) :
-1678.161191722  +10.000000000i 
 +859.134526218  +10.000000000i 
 -603.393715174 +445.986292551i 
 -603.393500034 -425.986292551i 
 +377.813880713  +10.000000000i 


 Press   return to continue
 Press X return to stop