Mathc matrices/c17b

Application

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

	gj_r2.c

/* ------------------------------------ */
/*  Save as :   gj_r2.c                 */
/* ------------------------------------ */
#include "v_a.h"
/* ------------------------------------ */
void fun(int r, int c)
{
double **A  = i_mR(r,c);
double **b  = i_mR(r,C1);
double **Ab = i_Abr_Ac_bc_mR(r,c,C1);
/*            i_Abr_Ac_bc_mR(R_Ab,C_A,C_b));  */

  r_mR(A,999.);
  r_mR(b,999.);
  
  clrscrn();
  printf(" A :");
  p_mR(A,S8,P0,C7);
  
  printf(" b :");
  p_mR(b,S8,P0,C7);
  
  printf(" Ab :");
  c_A_b_Ab_mR(A,b,Ab);
  p_mR(Ab,S8,P0,C7);
  stop();

  clrscrn();
  printf(" Copy/Past into the octave window.\n\n");
  p_Octave_mR(Ab,"Ab",P0);
  printf("\n rref(Ab,.00000000001)\n\n");
  printf("  gj_TP_mR(Ab) :");
  GJ_TP_mR(Ab);
  p_mR(Ab,S10,P4,C7);

  f_mR(Ab);
  f_mR(b);
  f_mR(A);
}
/* ------------------------------------ */
int main(void)
{
time_t t;

  srand(time(&t));
do
{
  fun(rp_I(R3)+R1, rp_I(R3)+R4);


} while(stop_w());

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

Nous utilisons la fonction GJ_TP_mR(Ab); pour résoudre ce système.

Exemple de sortie écran :

 ------------------------------------ 
 A :
    +816     +496     +678     -837     +126 
    +490     +624     +832     -285     +632 
    +140     -101     +464     +742     +884 

 b :
    +530 
    -895 
    -313 

 Ab :
    +816     +496     +678     -837     +126     +530 
    +490     +624     +832     -285     +632     -895 
    +140     -101     +464     +742     +884     -313 

 Press return to continue. 



 Copy/Past into the octave window.

 Ab=[
+816,+496,+678,-837,+126,+530;
+490,+624,+832,-285,+632,-895;
+140,-101,+464,+742,+884,-313]


 rref(Ab,.00000000001)

  gj_TP_mR(Ab) :
  +1.00000   +0.00000   +0.00000   -1.49815   -0.96054   +2.98320 
  +0.00000   +1.00000   -0.00000   -1.56188   -0.89873   -1.30001 
  +0.00000   +0.00000   +1.00000   +1.71119   +1.99936   -1.85765 


 Press return to continue
 Press X      to stop