Linear Algebra and the C Language/a081
Install and compile this file in your working directory.
/* ------------------------------------ */
/* Save as : c00c.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_PP_mR(Ab,YES) : 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;
}
/* ------------------------------------ */
/* ------------------------------------ */
Screen output example:
A :
-139 +570 +872 -301
+923 +274 +811 +460
+107 +982 +360 -48
+923 +274 +811 +460
b :
-431
+848
+180
+848
Ab :
-139 +570 +872 -301 -431
+923 +274 +811 +460 +848
+107 +982 +360 -48 +180
+923 +274 +811 +460 +848
Press return to continue.
Copy/Past into the octave window.
Ab=[
-139,+570,+872,-301,-431;
+923,+274,+811,+460,+848;
+107,+982,+360,-48,+180;
+923,+274,+811,+460,+848]
rref(Ab,.00000000001)
Press return to continue.
gj_PP_mR(Ab,YES) The system is consistent
+1.0000 +0.0000 +0.0000 +0.6910 +1.2372
+0.0000 +1.0000 +0.0000 -0.0500 +0.2070
+0.0000 +0.0000 +1.0000 -0.2024 -0.4324
+0.0000 +0.0000 +0.0000 +0.0000 +0.0000
Press return to continue
Press X return to stop