< Mathc matrices
C'est le même fichier que le précédent mais cette fois ci la matrice Ab à trois lignes identiques.
gj2_rsing.c ' |
---|
/* ------------------------------------ */
/* Save as : gj2_rsing.c */
/* ------------------------------------ */
#include "v_a.h"
/* ------------------------------------ */
void fun(int r)
{
/* i_AbR0_mR(R_Ab,C_A,C_b)); */
double **Ab = r_singular_Ab_mR(i_AbR0_mR( r, r, C1),999);
double **A = c_Ab_A_mR(Ab,i_mR(r,r));
double **b = c_Ab_b_mR(Ab,i_mR(r,C1));
clrscrn();
printf(" A :");
p_mR(A,8,0,7);
printf(" b :");
p_mR(b,8,0,7); /* I do not change the value of b */
printf(" Ab :");
p_mR(Ab,8,0,7);
getchar();
clrscrn();
printf(" Copy/Past into the octave window.\n\n");
p_octave_mR(Ab,"Ab",1,0);
printf("\n rref(Ab,.00000000001)\n\n");
printf(" gj_mR(Ab) :");
p_mR(gj_mR(Ab),10,5,7);
f_mR(Ab);
f_mR(b);
f_mR(A);
}
/* ------------------------------------ */
int main(void)
{
time_t t;
srand(time(&t));
do
{
fun(rp_I(5)+1);
} while(stop_w());
return 0;
}
/* ------------------------------------ */
/* ------------------------------------ */
C'est le même fichier que le précédent mais cette fois ci la matrice Ab à trois lignes identiques.
Exemple de sortie écran :
------------------------------------
A :
+316 +886 +222 -371 +190
+316 +886 +222 -371 +190
+916 -503 -839 +438 -275
-39 +850 +882 +564 +830
+316 +886 +222 -371 +190
b :
-369
-369
-155
-643
-369
Ab :
+316 +886 +222 -371 +190 -369
+316 +886 +222 -371 +190 -369
+916 -503 -839 +438 -275 -155
-39 +850 +882 +564 +830 -643
+316 +886 +222 -371 +190 -369
------------------------------------
Copy/Past into the octave window.
Ab=[
+316,+886,+222,-371,+190,-369;
+316,+886,+222,-371,+190,-369;
+916,-503,-839,+438,-275,-155;
-39,+850,+882,+564,+830,-643;
+316,+886,+222,-371,+190,-369]
rref(Ab,.00000000001)
gj_mR(Ab) :
+1.00000 +0.00000 +0.00000 +1.65269 +0.71587 -0.92022
+0.00000 +1.00000 +0.00000 -1.56451 -0.37520 +0.13789
+0.00000 +0.00000 +1.00000 +2.22028 +1.33428 -0.90260
+0.00000 +0.00000 +0.00000 +0.00000 +0.00000 +0.00000
+0.00000 +0.00000 +0.00000 +0.00000 +0.00000 +0.00000
Press return to continue
Press X to stop
Cet article est issu de Wikibooks. Le texte est sous licence Creative Commons - Attribution - Partage dans les Mêmes. Des conditions supplémentaires peuvent s'appliquer aux fichiers multimédias.