< Mathc matrices
Installer et compiler ces fichiers dans votre répertoire de travail.
 |
gj_r.c' |
|---|
/* ------------------------------------ */
/* Save as : gj_r.c */
/* ------------------------------------ */
#include "v_a.h"
/* ------------------------------------ */
void fun(int r)
{
double **A = i_mR(r,r);
double **b = i_mR(r,C1);
double **Ab = i_AbR0_mR(r,r,C1);
/* i_AbR0_mR(R_Ab,C_A,C_b)); */
do
{
r_mR(A,999.);
printf(".");
}while(!det_R(A));
clrscrn();
printf(" A :");
p_mR(A,8,0,7);
printf(" b :");
p_mR(r_mR(b,999.),8,0,7);
printf(" Ab :");
p_mR(c_A_b_Ab_mR(A,b,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;
}
/* ------------------------------------ */
/* ------------------------------------ */
Déclaration des fichiers h.
Exemple de sortie écran :
------------------------------------
A :
-223 +986 +854 -131 -853
+340 -395 +182 -471 -129
-617 -589 +264 -925 +328
+730 +200 -23 -969 -611
-943 +710 +992 +522 -237
b :
+870
+90
-555
-953
+76
Ab :
-223 +986 +854 -131 -853 +870
+340 -395 +182 -471 -129 +90
-617 -589 +264 -925 +328 -555
+730 +200 -23 -969 -611 -953
-943 +710 +992 +522 -237 +76
------------------------------------
Copy/Past into the octave window.
Ab=[
-223,+986,+854,-131,-853,+870;
+340,-395,+182,-471,-129,+90;
-617,-589,+264,-925,+328,-555;
+730,+200,-23,-969,-611,-953;
-943,+710,+992,+522,-237,+76]
rref(Ab,.00000000001)
gj_mR(Ab) :
+1.00000 +0.00000 +0.00000 +0.00000 +0.00000 -9.73243
-0.00000 +1.00000 -0.00000 -0.00000 +0.00000 -11.16218
+0.00000 +0.00000 +1.00000 +0.00000 +0.00000 -8.46177
+0.00000 +0.00000 +0.00000 +1.00000 +0.00000 +4.50061
-0.00000 -0.00000 -0.00000 -0.00000 +1.00000 -20.54104
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.