Mathc matrices/a222

Application

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

	c00a.c

/* ------------------------------------ */
/*  Save as :   c00a.c                  */
/* ------------------------------------ */
#include "v_a.h"
/* ------------------------------------ */
#define   RA R8
#define   CA C8
#define   Cb C1
/* ------------------------------------ */
int main(void)
{
double ab[RA*(CA+Cb)]={
//   I1    I2    I3    I4    I5    I6    
     -1,   +1,   +1,   +0,   +0,   +0,   +0, +0,   +0,
     +0,   +0,   -1,   +1,   -1,   +0,   +0, +0,   +0,
     +0,   +0,   +0,   -1,   +1,   +1,   +0, +0,   +0,
     +1,   -1,   +0,   +0,   +0,   -1,   +0, +0,   +0,
    +15,  +60,   +0,   +0,   +0,   +0,   +0, +0,  +90,
     +0,  -60,  +15,  +15,   +0,  +15,   +0, +0,   +0,
     +0,   +0,   +0,  -15,  -60,   +0,   +0, +0,  -90,
    +15,   +0,  +15,   +0,  -60,  +15,   +0, +0,   +0
};    

double **Ab = ca_A_mR(ab,i_Abr_Ac_bc_mR(RA,CA,Cb));
double **A  = c_Ab_A_mR(Ab,i_mR(RA,CA));
double **b  = c_Ab_b_mR(Ab,i_mR(RA,Cb));

  clrscrn();
  printf(" Ab :\n"
         "  I1     I2     I3     I4     I5     I6 ");
  p_mR(Ab,S6,P2,C9);
  getchar();

  clrscrn();
  printf(" Copy/Past into the octave window.\n\n");
  p_Octave_mR(Ab,"Ab",P0);
  printf("\n rref(Ab,.00000000001)\n\n");
  getchar();

  clrscrn();
  printf("  gj_TP_mR(Ab) :\n\n"
         "  I1     I2     I3     I4     I5     I6 ");
  gj_TP_mR(Ab);
  p_mR(Ab,S6,P2,C9);
  stop();

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

Exemple de sortie écran :

 ------------------------------------ 
 Ab :
  I1     I2     I3     I4     I5     I6 
 -1.00  +1.00  +1.00  +0.00  +0.00  +0.00  +0.00  +0.00  +0.00 
 +0.00  +0.00  -1.00  +1.00  -1.00  +0.00  +0.00  +0.00  +0.00 
 +0.00  +0.00  +0.00  -1.00  +1.00  +1.00  +0.00  +0.00  +0.00 
 +1.00  -1.00  +0.00  +0.00  +0.00  -1.00  +0.00  +0.00  +0.00 
+15.00 +60.00  +0.00  +0.00  +0.00  +0.00  +0.00  +0.00 +90.00 
 +0.00 -60.00 +15.00 +15.00  +0.00 +15.00  +0.00  +0.00  +0.00 
 +0.00  +0.00  +0.00 -15.00 -60.00  +0.00  +0.00  +0.00 -90.00 
+15.00  +0.00 +15.00  +0.00 -60.00 +15.00  +0.00  +0.00  +0.00 


 ------------------------------------ 
 Copy/Past into the octave window.

 Ab=[
-1,+1,+1,+0,+0,+0,+0,+0,+0;
+0,+0,-1,+1,-1,+0,+0,+0,+0;
+0,+0,+0,-1,+1,+1,+0,+0,+0;
+1,-1,+0,+0,+0,-1,+0,+0,+0;
+15,+60,+0,+0,+0,+0,+0,+0,+90;
+0,-60,+15,+15,+0,+15,+0,+0,+0;
+0,+0,+0,-15,-60,+0,+0,+0,-90;
+15,+0,+15,+0,-60,+15,+0,+0,+0]


 rref(Ab,.00000000001)


 ------------------------------------ 
  gj_TP_mR(Ab) :

  I1     I2     I3     I4     I5     I6 
 +1.00  -0.00  -0.00  +0.00  -0.00  +0.00  -0.00  -0.00  +2.00 
 +0.00  +1.00  +0.00  +0.00  +0.00  +0.00  +0.00  +0.00  +1.00 
 +0.00  +0.00  +1.00  +0.00  +0.00  +0.00  +0.00  +0.00  +1.00 
 +0.00  +0.00  +0.00  +1.00  +0.00  +0.00  +0.00  +0.00  +2.00 
 +0.00  +0.00  +0.00  +0.00  +1.00  +0.00  +0.00  +0.00  +1.00 
 +0.00  +0.00  +0.00  +0.00  +0.00  +1.00  +0.00  +0.00  +1.00 
 +0.00  +0.00  +0.00  +0.00  +0.00  +0.00  +0.00  +0.00  +0.00 
 +0.00  +0.00  +0.00  +0.00  +0.00  +0.00  +0.00  +0.00  +0.00 

 Press return to continue.