Application


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


c00a.c
/* ------------------------------------ */
/*  Save as :   c00.c                   */
/* ------------------------------------ */
#include "w_a.h"
/* ------------------------------------ */
/* ------------------------------------ */
#define RA              R4
#define CA              C6
#define FACTOR_E        +1.E-2   
/* ------------------------------------ */       
/* ------------------------------------ */
void fun(void)
{
int r = RA;
int c = CA;  

double a[RA*(CA*C2)] ={   
 -180,+815,  +439,-148,  +755,+959,  +130,-705,  +195,-395,  +525,+941, 
 -839, -97,  +402,-586,  +181,+311,  +794,+666,  +769,-689,  +115,+394, 
 -711, +59,  -927,+824,  +471,+237,  +705,+891,  -225,+942,  +660,-286, 
 -671,-765,  -806, +95,  +385,+909,  +589,+516,   +91,-651,  -576,+119   };

double tb[RA*(C1*C2)] ={   
 +462,+974, 
 +388,-988, 
 +672,-844, 
  +69,+433 };
                         
double **A  = ca_A_mZ(a, i_mZ(r,c));
double **b  = ca_A_mZ(tb, i_mZ(r,C1));
double **Ab = i_Abr_Ac_bc_mZ(r,r+C2,C1);

  clrscrn();
  printf(" A :");
  p_mZ(A, S5,P0, S4,P0, C6);
  printf(" b :");
  p_mZ(r_mZ(b,999.), S5,P0, S4,P0, C6);

  printf(" Copy/Past into the octave window.\n\n");
  p_Octave_mZ(c_A_b_Ab_mZ(A,b,Ab),"Ab",P0, P0);
  printf(" rref(Ab,.00000000001)\n\n");
  stop();

  clrscrn();
  printf(" gj_mZ(Ab) :");
  gj1_mZ(Ab);
  p_mZ(Ab, S8,P4, S8,P4, C3);
  stop();

  f_mZ(Ab);
  f_mZ(b);
  f_mZ(A);
}
/* ------------------------------------ */
int main(void)
{
  fun();
  return 0;
}
/* ------------------------------------ */
/* ------------------------------------ */



Exemple de sortie écran :
 ------------------------------------   
 A :
 -180+815i  +439-148i  +755+959i  +130-705i  +195-395i  +525+941i 
 -839 -97i  +402-586i  +181+311i  +794+666i  +769-689i  +115+394i 
 -711 +59i  -927+824i  +471+237i  +705+891i  -225+942i  +660-286i 
 -671-765i  -806 +95i  +385+909i  +589+516i   +91-651i  -576+119i 

 b :
 +479-154i 
 -501+992i 
 -763-591i 
 -843+708i 

 Copy/Past into the octave window.

 Ab=[
-180+815*i,+439-148*i,+755+959*i,+130-705*i,+195-395*i,+525+941*i,+479-154*i;
-839-97*i,+402-586*i,+181+311*i,+794+666*i,+769-689*i,+115+394*i,-501+992*i;
-711+59*i,-927+824*i,+471+237*i,+705+891*i,-225+942*i,+660-286*i,-763-591*i;
-671-765*i,-806+95*i,+385+909*i,+589+516*i,+91-651*i,-576+119*i,-843+708*i]

 rref(Ab,.00000000001)

 Press return to continue. 


------------------------------------  
 gj_mZ(Ab) :
 +1.0000 +0.0000i  -0.0000 +0.0000i  -0.0000 -0.0000i 
 +0.0000 +0.0000i  +1.0000 +0.0000i  +0.0000 -0.0000i 
 +0.0000 +0.0000i  -0.0000 +0.0000i  +1.0000 -0.0000i 
 +0.0000 +0.0000i  -0.0000 +0.0000i  +0.0000 +0.0000i 

 +0.0000 +0.0000i  +0.9166 +0.2968i  +0.2005 -1.2657i 
 +0.0000 +0.0000i  +0.8802 -0.0597i  -0.1948 +0.1861i 
 +0.0000 +0.0000i  +0.0832 -0.3955i  +0.3289 -0.0106i 
 +1.0000 +0.0000i  +0.9234 -0.4792i  -0.3506 -0.8870i 

 -1.2843 -0.6782i 
 -0.4263 +0.3803i 
 -0.5174 +0.2718i 
 -1.0200 +0.8650i 

 Press return to continue.