Matrice de Hankel

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

c00a.c
/* ------------------------------------ */
/*  Save as :   c00a.c                  */
/* ------------------------------------ */
#include "v_a.h"
/* ------------------------------------ */
#define FACTOR_E        +1.E-0         
/* ------------------------------------ */
void fun(int r)
{
double **U   = r_mR(          i_mR(R1, r),9);
double **V   = r_mR(          i_mR(r ,C1),9);	
double **A   = rHankel_mR(U,V,i_mR(r , r)  );

double **Adjoint  =           i_mR(r,r);
double **Cofactor =           i_mR(r,r);       

  clrscrn();
  printf(" A :");
  p_mR(A,S3,P0,C6);

  printf(" Cofactor :\n");
  p_mR(cofactor_mR(A,Cofactor),S6,P0,C6);
  
  printf(" Adjoint : Cofactor**T");
  p_mR(adjoint_mR(A,Adjoint),S6,P0,C6);

  f_mR(U);
  f_mR(V);      
  f_mR(A);
  
  f_mR(Cofactor);
  f_mR(Adjoint);
}
/* ------------------------------------ */
int main(void)
{
time_t t;

  srand(time(&t));

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

} while(stop_w());

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


Matrice de Hankel.

Exemple de sortie écran :

 A :
 +1  +3  -4  -2  -5 
 +3  -4  -2  -5  +2 
 -4  -2  -5  +2  -2 
 -2  -5  +2  -2  -3 
 -5  +2  -2  -3  +8 

 Cofactor :

  +416  -1734   +544  +2060  +1602 
 -1734  +1491  +1404   +684   -849 
  +544  +1404  +2180  -1712   -108 
 +2060   +684  -1712  +1940  +1416 
 +1602   -849   -108  +1416   -669 

 Adjoint : Cofactor**T
  +416  -1734   +544  +2060  +1602 
 -1734  +1491  +1404   +684   -849 
  +544  +1404  +2180  -1712   -108 
 +2060   +684  -1712  +1940  +1416 
 +1602   -849   -108  +1416   -669 


 Press   return to continue
 Press X return to stop