Mathc matrices/a260
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)
{
double **Ab = r_singular_Ab_mR( i_Abr_Ac_bc_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,S8,P0,C7);
printf(" b :");
p_mR(b,S8,P0,C7);
printf(" Ab :");
p_mR(c_A_b_Ab_mR(A,b,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");
stop();
clrscrn();
printf(" gj_mR(Ab) : The system is consistent");
gj_PP_mR(Ab,YES);
p_mR(Ab,S10,P4,C7);
is_system_consistent_mR(Ab);
f_mR(Ab);
f_mR(b);
f_mR(A);
}
/* ------------------------------------ */
int main(void)
{
time_t t;
srand(time(&t));
do
{
fun(rp_I(R5)+R1);
} while(stop_w());
return 0;
}
/* ------------------------------------ */
/* ------------------------------------ */
Exemple de sortie écran :
------------------------------------
A :
-73 +68 -727 -489 +418 -331
-411 -801 +770 +16 -623 +792
+954 -627 -377 +560 +460 -391
-73 +68 -727 -489 +418 -331
-121 +818 +448 +548 -755 -809
-73 +68 -727 -489 +418 -331
b :
+204
+272
-559
+204
-477
+204
Ab :
-73 +68 -727 -489 +418 -331 +204
-411 -801 +770 +16 -623 +792 +272
+954 -627 -377 +560 +460 -391 -559
-73 +68 -727 -489 +418 -331 +204
-121 +818 +448 +548 -755 -809 -477
-73 +68 -727 -489 +418 -331 +204
Press return to continue.
Copy/Past into the octave window.
Ab=[
-73,+68,-727,-489,+418,-331,+204;
-411,-801,+770,+16,-623,+792,+272;
+954,-627,-377,+560,+460,-391,-559;
-73,+68,-727,-489,+418,-331,+204;
-121,+818,+448,+548,-755,-809,-477;
-73,+68,-727,-489,+418,-331,+204]
rref(Ab,.00000000001)
Press return to continue.
gj_mR(Ab) : The system is consistent
+1.00000 +0.00000 +0.00000 +0.00000 +1.70527 +1.88088 +0.21927
+0.00000 +1.00000 +0.00000 +0.00000 +0.21183 -0.06474 -0.04469
+0.00000 +0.00000 +1.00000 +0.00000 +0.35489 +2.01953 +0.44709
+0.00000 +0.00000 +0.00000 +1.00000 -1.60754 -2.61534 -1.12081
+0.00000 +0.00000 +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