Application


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


c06.c
/* ------------------------------------ */
/*  Save as :   c06.c                   */
/* ------------------------------------ */
#include "v_a.h"
/* ------------------------------------ */
int main(void)
{
double ab[R4*C7]={
    +1,     -3,     +4,     -2,    +5,    +4,  0,
    +2,     -6,     +9,     -1,    +8,    +2,  0,
    +2,     -6,     +9,     -1,    +9,    +7,  0, 
    -1,      3,     -4,      2,    -5,    -4,  0,  
};

double **Ab = ca_A_mR(ab,i_Abr_Ac_bc_mR(R4,C6,C1));
double **A  = c_Ab_A_mR(Ab,i_mR(R4,C6));
double **b  = c_Ab_b_mR(Ab,i_mR(R4,C1));

double **Ab_free = i_Abr_Ac_bc_mR(csize_A_R(Ab),csize_A_R(Ab),C1+C3);

double **b_free  = i_mR(rsize_R(Ab_free),csize_b_R(Ab_free)); 

int r;

  clrscrn();
  printf(" A :");
  p_mR(A,S6,P1,C10);
  printf(" b :");
  p_mR(b,S6,P1,C10);
  printf(" Ab :");
  p_mR(Ab,S6,P1,C10);
  stop();

  clrscrn();
  printf(" Copy/Past into the octave window.\n\n");
  p_Octave_mR(Ab,"Ab",P1);
  printf("\n rref(Ab,.00000000001)\n\n");
  printf(" Ab :  gj_PP_mR(Ab,NO) :");
  gj_PP_mR(Ab,NO);
  p_mR(Ab,S6,P1,C10);
  stop();
  
  clrscrn();  
  put_zeroR_mR(Ab,Ab_free);  
  printf(" Ab_free : put_zeroR_mR(Ab,Ab_free);");  
  p_mR(Ab_free,S6,P1,C10);  
  stop();
  
  clrscrn();  
  put_freeV_mR(Ab_free);
  printf(" Ab_free : put_freeV_mR(Ab_free);");  
  p_mR(Ab_free,S6,P1,C10);  
  stop();
  
  clrscrn();
  printf(" Copy/Past into the octave window.\n\n");
  p_Octave_mR(Ab_free,"Ab_free",P1);
  printf("\n rref(Ab_free,.00000000001)\n\n");
  stop();
  
  clrscrn();  
  r = rsize_R(Ab_free);
  while(r>R1)    
        zero_below_pivot_gj1Ab_mR(Ab_free,r--);
        
  printf(" Ab_free : zero_below_pivot_gj1Ab_mR(Ab_free,r--);");  
  p_mR(Ab_free,S6,P1,C10);  
  stop();	
  
  clrscrn();
  c_Ab_b_mR(Ab_free,b_free);
  printf(" b_free :"); 
  p_mR(b_free,S10,P4,C7);
  printf(" b_free : free variables");
  p_freeV(b_free,S8,P4);
  stop();

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


Imprimer les variables libres.


Exemple de sortie écran :
 ------------------------------------ 
 A :
  +1.0   -3.0   +4.0   -2.0   +5.0   +4.0 
  +2.0   -6.0   +9.0   -1.0   +8.0   +2.0 
  +2.0   -6.0   +9.0   -1.0   +9.0   +7.0 
  -1.0   +3.0   -4.0   +2.0   -5.0   -4.0 

 b :
  +0.0 
  +0.0 
  +0.0 
  +0.0 

 Ab :
  +1.0   -3.0   +4.0   -2.0   +5.0   +4.0   +0.0 
  +2.0   -6.0   +9.0   -1.0   +8.0   +2.0   +0.0 
  +2.0   -6.0   +9.0   -1.0   +9.0   +7.0   +0.0 
  -1.0   +3.0   -4.0   +2.0   -5.0   -4.0   +0.0 

 Press return to continue. 

...



 Copy/Past into the octave window.

 Ab_free=[
+1.0,-3.0,+4.5,-0.5,+4.0,+1.0,+0.0,+0.0,+0.0,+0.0;
+0.0,+1.0,+0.0,+0.0,+0.0,+0.0,+0.0,+1.0,+0.0,+0.0;
-0.0,-0.0,+1.0,+3.0,-2.0,-6.0,-0.0,+0.0,+0.0,+0.0;
+0.0,+0.0,+0.0,+1.0,+0.0,+0.0,+0.0,+0.0,+1.0,+0.0;
+0.0,+0.0,+0.0,+0.0,+1.0,+5.0,+0.0,+0.0,+0.0,+0.0;
+0.0,+0.0,+0.0,+0.0,+0.0,+1.0,+0.0,+0.0,+0.0,+1.0]


 rref(Ab_free,.00000000001)

 Press return to continue. 


 Ab_free : zero_below_pivot_gj1Ab_mR(Ab_free,r--);
  +1.0   +0.0   +0.0   +0.0   +0.0   +0.0   +0.0   +3.0  +14.0  +37.0 
  +0.0   +1.0   +0.0   +0.0   +0.0   +0.0   +0.0   +1.0   +0.0   +0.0 
  +0.0   +0.0   +1.0   +0.0   +0.0   +0.0   +0.0   +0.0   -3.0   -4.0 
  +0.0   +0.0   +0.0   +1.0   +0.0   +0.0   +0.0   +0.0   +1.0   +0.0 
  +0.0   +0.0   +0.0   +0.0   +1.0   +0.0   +0.0   +0.0   +0.0   -5.0 
  +0.0   +0.0   +0.0   +0.0   +0.0   +1.0   +0.0   +0.0   +0.0   +1.0 

 Press return to continue. 


 b_free :
   +0.0000    +3.0000   +14.0000   +37.0000 
   +0.0000    +1.0000    +0.0000    +0.0000 
   +0.0000    +0.0000    -3.0000    -4.0000 
   +0.0000    +0.0000    +1.0000    +0.0000 
   +0.0000    +0.0000    +0.0000    -5.0000 
   +0.0000    +0.0000    +0.0000    +1.0000 

 b_free : free variables
 x1 =   +0.0000  +3.0000*s +14.0000*t +37.0000*u
 x2 =   +0.0000  +1.0000*s  +0.0000*t  +0.0000*u
 x3 =   +0.0000  +0.0000*s  -3.0000*t  -4.0000*u
 x4 =   +0.0000  +0.0000*s  +1.0000*t  +0.0000*u
 x5 =   +0.0000  +0.0000*s  +0.0000*t  -5.0000*u
 x6 =   +0.0000  +0.0000*s  +0.0000*t  +1.0000*u
 Press return to continue.