Mathc matrices/c17b
Installer et compiler ces fichiers dans votre répertoire de travail.
c00a.c |
---|
/* ------------------------------------ */
/* Save as : c00a.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