Mathc complexes/058

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

	`c00c.c`

/* ------------------------------------ */
/*  Save as :   c00c.c                  */
/* ------------------------------------ */
#include "w_a.h"
/* ------------------------------------ */
#define   RA R5
#define   CA C5
#define   Cb C1
/* ------------------------------------ */
int main(void)
{
double ab[RA*(CA+Cb)]={
// x2    x4    x6    x7             b
   -1,   +0,   +0,   +0,   +0,    (+20-50),
   +1,   +1,   +0,   +0,   +0,    (+60),  
   +0,   +0,   +1,   -1,   +0,    (-60),
   +0,   +0,   +0,   +1,   +0,    (+90-20),
   +0,   -1,   -1,   +0,   +0,    (+50-90) 
};

double **Ab =  ca_A_mRZ(ab,i_Abr_Ac_bc_mZ(RA,CA,Cb));
double **A  = c_Ab_A_mZ(Ab,i_mZ(RA,CA));
double **b  = c_Ab_b_mZ(Ab,i_mZ(RA,Cb));

  clrscrn();
  printf(" A :");
  p_mRZ(A,S5,P0,C7);
  printf(" b :");
  p_mRZ(b,S5,P0,C7);
  printf(" Ab :");
  p_mRZ(Ab,S5,P0,C7);
  getchar();

  clrscrn();
  printf(" Copy/Paste into the octave window.\n\n");
  p_Octave_mZ(Ab,"Ab",P0,P0);
  printf("\n rref(Ab,.00000000001)\n\n");
  printf("  gj_mZ(Ab) :\n\n"
      "   x2    x4    x6    x7 ");
  gj_mZ(Ab);
  p_mRZ(Ab,S5,P0,C7);
  stop();

  f_mZ(Ab);
  f_mZ(b);
  f_mZ(A);
  
  return 0;
}
/* ------------------------------------ */
/* ------------------------------------ */

Exemple de sortie écran :

 A :
   -1    +0    +0    +0    +0 
   +1    +1    +0    +0    +0 
   +0    +0    +1    -1    +0 
   +0    +0    +0    +1    +0 
   +0    -1    -1    +0    +0 

 b :
  -30 
  +60 
  -60 
  +70 
  -40 

 Ab :
   -1    +0    +0    +0    +0   -30 
   +1    +1    +0    +0    +0   +60 
   +0    +0    +1    -1    +0   -60 
   +0    +0    +0    +1    +0   +70 
   +0    -1    -1    +0    +0   -40 


 Copy/Paste into the octave window.

 Ab=[
-1+0*i,+0+0*i,+0+0*i,+0+0*i,+0+0*i,-30+0*i;
+1+0*i,+1+0*i,+0+0*i,+0+0*i,+0+0*i,+60+0*i;
+0+0*i,+0+0*i,+1+0*i,-1+0*i,+0+0*i,-60+0*i;
+0+0*i,+0+0*i,+0+0*i,+1+0*i,+0+0*i,+70+0*i;
+0+0*i,-1+0*i,-1+0*i,+0+0*i,+0+0*i,-40+0*i]


 rref(Ab,.00000000001)

  gj_mZ(Ab) :

   x2    x4    x6    x7 
   +1    +0    +0    +0    +0   +30 
   +0    +1    +0    +0    +0   +30 
   +0    +0    +1    +0    +0   +10 
   +0    +0    +0    +1    +0   +70 
   +0    +0    +0    +0    +0    +0 

 Press return to continue.