Exemple de sortie écran :
 A :
 -944.0000 +337.0000 -777.0000 -131.0000 
 -944.0000 +337.0000 -777.0000 -131.0000 
 +509.0000 +581.0000 -475.0000 -544.0000 
 +69.0000 +777.0000 -245.0000 -51.0000 
 b :
 -592.0000 
 -592.0000 
 -649.0000 
 +974.0000 
 Ab :
 -944.0000 +337.0000 -777.0000 -131.0000 -592.0000 
 -944.0000 +337.0000 -777.0000 -131.0000 -592.0000 
 +509.0000 +581.0000 -475.0000 -544.0000 -649.0000 
 +69.0000 +777.0000 -245.0000 -51.0000 +974.0000 
 Press return to continue. 
 ---------------------
 Copy/Past into the octave window.
 Ab=[
-944.0000,+337.0000,-777.0000,-131.0000,-592.0000;
-944.0000,+337.0000,-777.0000,-131.0000,-592.0000;
+509.0000,+581.0000,-475.0000,-544.0000,-649.0000;
+69.0000,+777.0000,-245.0000,-51.0000,+974.0000]
 rref(Ab.00000000001)
 Ab : gj_PP_mR(Ab,NO) :
 +1.0000 -0.3570 +0.8231 +0.1388 +0.6271 
 +0.0000 +1.0000 -0.3765 -0.0756 +1.1610 
 -0.0000 -0.0000 +1.0000 +0.9179 +3.0549 
 +0.0000 +0.0000 +0.0000 +0.0000 +0.0000 
 Press return to continue. 
 ---------------------
 Ab_free : put_var_mR(Ab,Ab_free);
 +1.0000 -0.3570 +0.8231 +0.1388 +0.6271 +0.0000 
 +0.0000 +1.0000 -0.3765 -0.0756 +1.1610 +0.0000 
 -0.0000 -0.0000 +1.0000 +0.9179 +3.0549 +0.0000 
 +0.0000 +0.0000 +0.0000 +0.0000 +0.0000 +0.0000 
 Ab_free : put_freeV_mR(Ab_free);
 +1.0000 -0.3570 +0.8231 +0.1388 +0.6271 +0.0000 
 +0.0000 +1.0000 -0.3765 -0.0756 +1.1610 +0.0000 
 -0.0000 -0.0000 +1.0000 +0.9179 +3.0549 +0.0000 
 +0.0000 +0.0000 +0.0000 +1.0000 +0.0000 +1.0000 
 Press return to continue. 
 ---------------------
 Copy/Past into the octave window.
 Ab_free=[
+1.0000,-0.3570,+0.8231,+0.1388,+0.6271,+0.0000;
+0.0000,+1.0000,-0.3765,-0.0756,+1.1610,+0.0000;
-0.0000,-0.0000,+1.0000,+0.9179,+3.0549,+0.0000;
+0.0000,+0.0000,+0.0000,+1.0000,+0.0000,+1.0000]
 rref(Ab_free.00000000001)
 Ab_free : zero_below_pivot_gj1Ab_mR(Ab_free,r--);
 +1.0000 +0.0000 +0.0000 +0.0000 -1.0623 +0.5204 
 +0.0000 +1.0000 +0.0000 +0.0000 +2.3111 -0.2700 
 -0.0000 -0.0000 +1.0000 +0.0000 +3.0549 -0.9179 
 +0.0000 +0.0000 +0.0000 +1.0000 +0.0000 +1.0000 
 Press return to continue. 
 ---------------------
 b_free :
 -1.0623 +0.5204 
 +2.3111 -0.2700 
 +3.0549 -0.9179 
 +0.0000 +1.0000 
 b_free : free variables
 x1 = -1.0623 +0.5204*u
 x2 = +2.3111 -0.2700*u
 x3 = +3.0549 -0.9179*u
 x4 = +0.0000 +1.0000*u
 Press return to continue.

Mathc matrices/a117
Application 
Installer et compiler ces fichiers dans votre répertoire de travail.
/* ------------------------------------ */
/* Save as : c00a.c */
int main(void)
double ab[RA*(CA+C1)]={
 +692.0, +524.0, +86.0, -1156.0, 
 +346.0, +262.0, +43.0, -578.0, 
-1038.0, -786.0, -129.0, +1734.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));
double **Ab_free = i_Abr_Ac_bc_mR(csize_A_R(Ab),csize_A_R(Ab),Cb+CFV);
double **b_free = i_mR(rsize_R(Ab_free),csize_b_R(Ab_free)); 
int r;
 clrscrn();
 printf(" A :");
 p_mR(A,S10,P4,C10);
 printf(" b :");
 p_mR(b,S10,P4,C10);
 printf(" Ab :");
 p_mR(Ab,S10,P4,C10);
 stop();
 clrscrn();
 printf(" Copy/Past into the octave window.\n\n");
 p_Octave_mR(Ab,"Ab",P4);
 printf("\n rref(Ab.00000001)\n\n");
 printf(" Ab : gj_PP_mR(Ab,NO) :");
 gj_PP_mR(Ab,NO);
 p_mR(Ab,S10,P4,C10);
 stop();
 clrscrn(); 
 put_var_mR(Ab,Ab_free); 
 printf(" Ab_free : put_var_mR(Ab,Ab_free);"); 
 p_mR(Ab_free,S10,P4,C10); 
 put_freeV_mR(Ab_free);
 printf(" Ab_free : put_freeV_mR(Ab_free);"); 
 p_mR(Ab_free,S10,P4,C10); 
 stop();
 clrscrn();
 printf(" Copy/Past into the octave window.\n\n");
 p_Octave_mR(Ab_free,"Ab_free",P4);
 printf("\n rref(Ab_free.00000001)\n\n");
 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,S10,P4,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;
/* ------------------------------------ */

Exemple de sortie écran :
 A :
 +692.0000 +524.0000 +86.0000 
 +346.0000 +262.0000 +43.0000 
-1038.0000 -786.0000 -129.0000 
 b :
-1156.0000 
 -578.0000 
+1734.0000 
 Ab :
 +692.0000 +524.0000 +86.0000 -1156.0000 
 +346.0000 +262.0000 +43.0000 -578.0000 
-1038.0000 -786.0000 -129.0000 +1734.0000 
 Press return to continue. 
 ---------------------
 Copy/Past into the octave window.
 Ab=[
+692.0000,+524.0000,+86.0000,-1156.0000;
+346.0000,+262.0000,+43.0000,-578.0000;
-1038.0000,-786.0000,-129.0000,+1734.0000]
 rref(Ab.00000001)
 Ab : gj_PP_mR(Ab,NO) :
 +1.0000 +0.7572 +0.1243 -1.6705 
 +0.0000 +0.0000 +0.0000 +0.0000 
 +0.0000 +0.0000 +0.0000 +0.0000 
 Press return to continue. 
 ---------------------
 Ab_free : put_var_mR(Ab,Ab_free);
 +1.0000 +0.7572 +0.1243 -1.6705 +0.0000 +0.0000 
 +0.0000 +0.0000 +0.0000 +0.0000 +0.0000 +0.0000 
 +0.0000 +0.0000 +0.0000 +0.0000 +0.0000 +0.0000 
 Ab_free : put_freeV_mR(Ab_free);
 +1.0000 +0.7572 +0.1243 -1.6705 +0.0000 +0.0000 
 +0.0000 +1.0000 +0.0000 +0.0000 +1.0000 +0.0000 
 +0.0000 +0.0000 +1.0000 +0.0000 +0.0000 +1.0000 
 Press return to continue. 
 ---------------------
 Copy/Past into the octave window.
 Ab_free=[
+1.0000,+0.7572,+0.1243,-1.6705,+0.0000,+0.0000;
+0.0000,+1.0000,+0.0000,+0.0000,+1.0000,+0.0000;
+0.0000,+0.0000,+1.0000,+0.0000,+0.0000,+1.0000]
 rref(Ab_free.00000001)
 Ab_free : zero_below_pivot_gj1Ab_mR(Ab_free,r--);
 +1.0000 +0.0000 +0.0000 -1.6705 -0.7572 -0.1243 
 +0.0000 +1.0000 +0.0000 +0.0000 +1.0000 +0.0000 
 +0.0000 +0.0000 +1.0000 +0.0000 +0.0000 +1.0000 
 Press return to continue. 
 ---------------------
 b_free :
 -1.6705 -0.7572 -0.1243 
 +0.0000 +1.0000 +0.0000 
 +0.0000 +0.0000 +1.0000 
 b_free : free variables
 x1 = -1.6705 -0.7572*v -0.1243*w
 x2 = +0.0000 +1.0000*v +0.0000*w
 x3 = +0.0000 +0.0000*v +1.0000*w
 Press return to continue.

Mathc matrices/a118
Application
Installer et compiler ces fichiers dans votre répertoire de travail.
/* ------------------------------------ */
/* Save as : c00a.c */
int main(void)
double ab[RA*(CA+C1)]={
-1888.0, +674.0, -1554.0, -262.0, -1184.0, 
+2832.0, -1011.0, +2331.0, +393.0, +1776.0, 
 -944.0, +337.0, -777.0, -131.0, -592.0, 
 +69.0, +777.0, -245.0, -51.0, +974.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));
double **Ab_free = i_Abr_Ac_bc_mR(csize_A_R(Ab),csize_A_R(Ab),Cb+CFV);
double **b_free = i_mR(rsize_R(Ab_free),csize_b_R(Ab_free)); 
int r;
 clrscrn();
 printf(" A :");
 p_mR(A,S10,P4,C10);
 printf(" b :");
 p_mR(b,S10,P4,C10);
 printf(" Ab :");
 P_mR(Ab,S10,P4,C10);
 stop();
 clrscrn();
 printf(" Copy/Past into the octave window.\n\n");
 p_Octave_mR(Ab,"Ab",P4);
 printf("\n rref(Ab.00000001)\n\n");
 printf(" Ab : gj_PP_mR(Ab,NO) :");
 gj_PP_mR(Ab,NO);
 p_mR(Ab,S10,P4,C10);
 stop();
 clrscrn(); 
 put_var_mR(Ab,Ab_free); 
 printf(" Ab_free : put_var_mR(Ab,Ab_free);"); 
 p_mR(Ab_free,S10,P4,C10); 
 put_freeV_mR(Ab_free);
 printf(" Ab_free : put_freeV_mR(Ab_free);"); 
 p_mR(Ab_free,S10,P4,C10); 
 stop();
 clrscrn();
 printf(" Copy/Past into the octave window.\n\n");
 p_Octave_mR(Ab_free,"Ab_free",P4);
 printf("\n rref(Ab_free.00000001)\n\n");
 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,S10,P4,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;
/* ------------------------------------ */