Mathc matrices/c15a
Application 
Installer et compiler ces fichiers dans votre répertoire de travail.
/* ------------------------------------ */
/* Save as : gj_r1.c */
void fun(int r)
double **A = i_mR(r,r);
double **b = i_mR(r,C1);
double **Ab = i_Abr_Ac_bc_mR(r,r,C1);
/* i_Abr_Ac_bc_mR(R_Ab,C_A,C_b)); */
 do
 r_mR(A,999.);
 printf(".");
 }while(!det_R(A));
 clrscrn();
 printf(" A :");
 p_mR(A,S8,P0,C7);
 printf(" b :");
 r_mR(b,999.);
 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(RC5)+R1);
} while(stop_w());
 return 0;
/* ------------------------------------ */
Nous utilisons la fonction gj_TP_mR(Ab); pour résoudre ce système.

Exemple de sortie écran :
 A :
 -325 -601 -801 -221 +360 
 -855 +832 +226 -815 -233 
 +724 -933 -49 -231 -965 
 +290 +238 -905 +772 -327 
 +636 -519 +662 -7 -375 
 b :
 +524 
 -143 
 -787 
 +446 
 -813 
 Ab :
 -325 -601 -801 -221 +360 +524 
 -855 +832 +226 -815 -233 -143 
 +724 -933 -49 -231 -965 -787 
 +290 +238 -905 +772 -327 +446 
 +636 -519 +662 -7 -375 -813 
 Press return to continue. 
 Copy/Past into the octave window.
 Ab=[
-325,-601,-801,-221,+360,+524;
-855,+832,+226,-815,-233,-143;
+724,-933,-49,-231,-965,-787;
+290,+238,-905,+772,-327,+446;
+636,-519,+662,-7,-375,-813]
 rref(Ab.00000000001)
 gj_TP_mR(Ab) :
 +1.00000 -0.00000 -0.00000 +0.00000 -0.00000 +2.17421 
 +0.00000 +1.00000 +0.00000 -0.00000 +0.00000 +1.33393 
 +0.00000 +0.00000 +1.00000 -0.00000 +0.00000 -1.37788 
 +0.00000 +0.00000 +0.00000 +1.00000 +0.00000 -1.58506 
 +0.00000 +0.00000 +0.00000 +0.00000 +1.00000 +1.60646 
 Press return to continue
 Press X to stop

Mathc matrices/c15b
Application 
Installer et compiler ces fichiers dans votre répertoire de travail.
/* ------------------------------------ */
/* Save as : gj_r2.c */
void fun(int r)
double **A = i_mR(r,r);
double **b = i_mR(r,C5);
double **Ab = i_Abr_Ac_bc_mR(r,r,C5);
/* i_Abr_Ac_bc_mR(R_Ab,C_A,C_b)); */
 do
 r_mR(A,999.);
 printf(".");
 }while(!det_R(A));
 clrscrn();
 printf(" A :");
 p_mR(A,S8,P0,C6);
 printf(" b1 b2 ... bn :\n");
 r_mR(b,999.);
 p_mR(b,S8,P0,C6);
 stop();

clrscrn(); 
 printf(" Ab :");
 c_A_b_Ab_mR(A,b,Ab);
 p_mR(Ab,S8,P0,C6);
 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,C6);
 f_mR(Ab);
 f_mR(b);
 f_mR(A);
/* ------------------------------------ */
int main(void)
time_t t;
 srand(time(&t));
do
 fun(rp_I(RC5)+R1);
} while(stop_w());
 return 0;
/* ------------------------------------ */
Nous utilisons la fonction gj_TP_mR(Ab); pour résoudre ce système.

Exemple de sortie écran :
 A :
 +788 -473 -703 +478 -281 
 +412 +754 +344 -377 -297 
 -515 -621 -443 +388 -755 
 +2 -359 -789 +882 +958 
 -551 +376 +138 +280 -117 
 b1 b2 ... bn :
 -147 -895 -377 -363 -199 
 +934 +152 +390 +356 +348 
 -389 +486 +820 -733 -863 
 -117 +966 +484 -559 +354 
 -957 +560 -431 +746 +162 
 Press return to continue. 
 Ab :
 +788 -473 -703 +478 -281 -147 
 +412 +754 +344 -377 -297 +934 
 -515 -621 -443 +388 -755 -389 
 +2 -359 -789 +882 +958 -117 
 -551 +376 +138 +280 -117 -957 
 -895 -377 -363 -199 
 +152 +390 +356 +348 
 +486 +820 -733 -863 
 +966 +484 -559 +354 
 +560 -431 +746 +162 
 Press return to continue. 
 Copy/Past into the octave window.
 Ab=[
+788,-473,-703,+478,-281,-147,-895,-377,-363,-199;
+412,+754,+344,-377,-297,+934,+152,+390,+356,+348;
-515,-621,-443,+388,-755,-389,+486,+820,-733,-863;
+2,-359,-789,+882,+958,-117,+966,+484,-559,+354;
-551,+376,+138,+280,-117,-957,+560,-431,+746,+162]
 rref(Ab.00000000001)
 gj_TP_mR(Ab) :
 +1.00000 -0.00000 +0.00000 +0.00000 +0.00000 -1.10815 
 +0.00000 +1.00000 +0.00000 +0.00000 -0.00000 +2.22788 
 -0.00000 -0.00000 +1.00000 -0.00000 +0.00000 -6.30681 
 +0.00000 +0.00000 +0.00000 +1.00000 +0.00000 -5.31055 
 +0.00000 -0.00000 +0.00000 +0.00000 +1.00000 +0.41008 
 -2.37746 -2.59443 +1.06347 +0.34143 
 +2.55587 +2.83566 -0.54926 +0.50154 
 -5.47986 -8.54634 +4.37136 +0.49836 
 -3.22963 -6.15601 +3.25975 +0.53116 
 +0.43136 +0.20223 -0.19249 +0.47818 
 Press return to continue
 Press X to stop

Mathc matrices/c15c
Application 
Installer et compiler ces fichiers dans votre répertoire de travail.
/* Save as : c00a.c */
void fun(int r)
double **A = i_mR(r,r);
double **b = i_mR(r,r);
double **Ab = i_Abr_Ac_bc_mR( r, r, r);
/* i_Abr_Ac_bc_mR(R_Ab,C_A,C_b)); */
double **invA = i_mR(r,r);
 do
 r_mR(A,999.);
 printf(".");
 }while(!det_R(A));
 clrscrn();
 printf(" A :");
 p_mR(A,S8,P0,C4);
 printf(" b : ID :");
 eye_mR(b);
 p_mR(b,S8,P0,C4);
 stop();
 clrscrn(); 
 printf(" Ab :");
 c_A_b_Ab_mR(A,b,Ab);
 p_mR(Ab,S8,P0,C4);
 printf(" gj_TP_mR(Ab) :"); 
 gj_TP_mR(Ab);
 pE_mR(Ab,S14,P4,C4);
 stop();
 clrscrn(); 
 printf(" Copy/Past into the octave window.\n\n");
 p_Octave_mR(A,"A",P0);
 printf("format short e\n");
 printf(" inv(A)\n\n");
 printf(" invA : c_Inv_A_mR(Ab,invA);"); 
 c_Inv_A_mR(Ab,invA);
 pE_mR(invA,S14,P4,C4);
 f_mR(invA);
 f_mR(Ab);
 f_mR(b);
 f_mR(A);
/* ------------------------------------ */
int main(void)
time_t t;
 srand(time(&t));
do
 fun(R4);
} while(stop_w());
 return 0;
/* ------------------------------------ */
 Utiliser la fonction gj_TP_mR(Ab); pour inverser une matrice.
 Exemple de sortie écran :
 A :
 -828 +6 +321 +504 
 +555 +300 +507 -830 
 +303 +959 +719 -612 
 -458 -594 -64 -25 
 b : ID :
 +1 +0 +0 +0 
 +0 +1 +0 +0 
 +0 +0 +1 +0 
 +0 +0 +0 +1 
 Press return to continue. 
 ------------------------------------
 Ab :
 -828 +6 +321 +504 
 +555 +300 +507 -830 
 +303 +959 +719 -612 
 -458 -594 -64 -25 
 +1 +0 +0 +0 
 +0 +1 +0 +0 
 +0 +0 +1 +0 
 +0 +0 +0 +1 
 gj_TP_mR(Ab) :
 +1.0000e+00 +0.0000e+00 +0.0000e+00 +0.0000e+00 
 +0.0000e+00 +1.0000e+00 +0.0000e+00 +0.0000e+00 
 +0.0000e+00 +0.0000e+00 +1.0000e+00 +0.0000e+00 
 +0.0000e+00 +0.0000e+00 +0.0000e+00 +1.0000e+00 
 +9.6778e-04 +2.4025e-03 -2.3558e-03 -2.5803e-03 
 -1.1140e-03 -2.3055e-03 +2.1822e-03 +6.6390e-04 
 +2.6791e-03 +3.5312e-03 -2.4923e-03 -2.2132e-03 
 +1.8810e-03 +1.7253e-03 -2.3089e-03 -2.8374e-03 
 Press return to continue. 
 ------------------------------------
 Copy/Past into the octave window.
 A=[
-828,+6,+321,+504;
+555,+300,+507,-830;
+303,+959,+719,-612;
-458,-594,-64,-25]
format short e
 inv(A)
 invA : c_Inv_A_mR(Ab,invA);
 +9.6778e-04 +2.4025e-03 -2.3558e-03 -2.5803e-03 
 -1.1140e-03 -2.3055e-03 +2.1822e-03 +6.6390e-04 
 +2.6791e-03 +3.5312e-03 -2.4923e-03 -2.2132e-03 
 +1.8810e-03 +1.7253e-03 -2.3089e-03 -2.8374e-03 
 Press return to continue
 Press X return to stop

Mathc matrices/c15d
Application 
Installer et compiler ces fichiers dans votre répertoire de travail.
/* ------------------------------------ */
/* Save as : gj_r_a.c */
int main(void)
double ab[RA*(CA+Cb)]={
 +964, -689, +474, +948, -159, 
 +640, -693, -937, +656, +718, 
 -981, -337, -287, +420, +812, 
 -807, -123, +404, -177, -105 
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(" A :");
 p_mR(A,S8,P0,C7);
 printf(" b :");
 p_mR(b,S8,P0,C7);
 printf(" Ab :");
 p_mR(Ab,S8,P0,C7);
 getchar();
 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);
 stop();
 f_mR(Ab);
 f_mR(b);
 f_mR(A);
 return 0;
/* ------------------------------------ */
Nous utilisons la fonction gj_TP_mR(Ab); pour résoudre ce système.