Mathc complexes/a305

Quelques propriétés des valeurs propres. Matrices non symétrique conjuguée

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

	c02a.c

/* ------------------------------------ */
/*  Save as :  c02a.c                   */
/* ------------------------------------ */
#include "w_a.h"
/* ------------------------------------ */
void fun(int r)
{
double **A         = r_mZ(i_mZ(r,r),9.);
double **EigsValue =      i_mZ(r,C1);

nb_Z a =i_Z(1,0);
nb_Z b =i_Z(1,0);

int  i = 0;

  clrscrn();

  printf(" A :");
  p_mZ(A,S5,P0,S4,P0,C6);
       
  printf(" EigsValue :");
  p_mZ(eigs_mZ(A,EigsValue),S10,P3,S10,P3,C6);   
  
  for(i=R1;i<=r;i++)
     {
      a = i_Z(EigsValue[i][C1],EigsValue[i][C2]);
      b = mul_Z(a,b);
      }
        
  printf(" Product of the eigenvalues of A : "); 
  p_Z(b,S10,P3,S10,P3); printf("\n\n");  
  printf(" Determinant of A                : "); 
  p_Z(det_Z(A),S10,P3,S10,P3); printf("\n\n");
   
  f_mZ(A);
  f_mZ(EigsValue);
}
/* ------------------------------------ */
int main(void)
{
time_t t;

  srand(time(&t));

do
{
 fun(rp_I(R3)+R1);

} while(stop_w());

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

Le produit des valeurs propres de A est égale au déterminant de A :


Exemple de sortie écran :

 A :
   -7  -5i    -2  -8i    +4  -4i    -5  -9i 
   +1  +9i    -5  +3i    -9  -7i    -9  +5i 
   +2  +9i    +9  +5i    -1  +5i    +7  -6i 
   +2  -8i    -6  -9i    +3  -7i    +6  +4i 

 EigsValue :
   -14.511    +2.067i 
    +8.777    +5.334i 
    -5.361    +8.056i 
    +4.095    -8.457i 

 Product of the eigenvalues of A :  -1748.000-13576.000i 

 Determinant of A                :  -1748.000-13576.000i 


 Press   return to continue
 Press X return to stop