Linear Algebra and the C Language/a0ip

__NOTOC__ Install and compile this file in your working directory. /* ------------------------------------ */ /* Save as : c00d.c */ /* ------------------------------------ */ #include "v_a.h" /* ------------------------------------ */ #define RA R3 #define CA C3 #define Cb C5 /* ------------------------------------ */ int main(void) { double at[RA*CA]={ +10, -150, -148, +207, -215, +997, -858, +803, +289 }; double bt[RA*Cb]={ +437, -146, -175, -954, -432, -233, +243, +772, +920, -98, +444, -886, -604, -15, +946 }; /* --------- Copy "at" into A ------- */ double **A = ca_A_mR(at, i_mR(RA,CA)); /* --------- Copy "bt" into b ------- */ double **b = ca_A_mR(bt, i_mR(RA,Cb)); /* --------- Initialize Ab ----------- */ double **Ab = i_Abr_Ac_bc_mR(RA,CA,Cb); /* ------- Copy A and b into Ab ----- */ c_A_b_Ab_mR(A,b,Ab); /* ------------------------------------ */ clrscrn(); printf(" A :"); p_mR(A, S8,P0,C6); printf(" b1 b2 ... bn:"); p_mR(b, S8,P0,C6); printf(" Ab:"); p_mR(Ab, S8,P0,C8); stop(); clrscrn(); printf(" Copy/Paste into the octave window.\n\n"); p_Octave_mR(Ab,"Ab",P0); printf(" rref(Ab,.00000000001)\n\n\n\n"); printf(" gj_TP_mR(Ab):"); gj_TP_mR(Ab); p_mR(Ab, S8,P4,C8); stop(); f_mR(Ab); f_mR(b); f_mR(A); return 0; } /* ------------------------------------ */ /* ------------------------------------ */ Screen output example: A : +10 -150 -148 +207 -215 +997 -858 +803 +289 b1 b2 ... bn: +437 -146 -175 -954 -432 -233 +243 +772 +920 -98 +444 -886 -604 -15 +946 Ab: +10 -150 -148 +437 -146 -175 -954 -432 +207 -215 +997 -233 +243 +772 +920 -98 -858 +803 +289 +444 -886 -604 -15 +946 Press return to continue. Copy/Paste into the octave window. Ab=[ +10,-150,-148,+437,-146,-175,-954,-432; +207,-215,+997,-233,+243,+772,+920,-98; -858,+803,+289,+444,-886,-604,-15,+946] rref(Ab,.00000000001) gj_TP_mR(Ab): +1.0000 +0.0000 -0.0000 -3.3492 +2.0425 +1.5396 +5.7575 +1.5874 -0.0000 +1.0000 -0.0000 -2.9619 +1.0616 +0.6767 +5.7825 +2.8101 +0.0000 +0.0000 +1.0000 -0.1771 +0.0486 +0.6006 +0.9744 +0.1781 Press return to continue.