Linear Algebra and the C Language/a0k3

__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 C1 /* ------------------------------------ */ int main(void) { double ab[RA*(CA+Cb)]={ +4.000, -2.000, -7.000, -165.000, +9.000, -4.000, -2.000, -14.000, +1.000, +6.000, +9.000, +403.000 }; double **Ab = ca_A_mR(ab,i_Abr_Ac_bc_mR(RA,CA,Cb)); double **A = c_Ab_A_mR(Ab,i_mR(RA,CA)); double **b = c_Ab_b_mR(Ab,i_mR(RA,Cb)); clrscrn(); printf(" You have this nonlinear system of equations :\n\n"); printf(" a x**4 + b y**3 + c z**2 = d \n"); printf(" e x**4 + f y**3 + g z**2 = h \n"); printf(" i x**4 + j y**3 + k z**2 = l \n\n"); printf(" We write X = x**4; Y = y**3; Z = z**2;\n\n"); printf(" a X + b Y + c Z = d\n"); printf(" e X + f Y + g Z = h\n"); printf(" i X + j Y + k Z = l\n\n"); printf(" Ab :\n"); printf(" a b c d\n"); printf(" e f g h\n"); printf(" i j k l\n\n\n"); stop(); clrscrn(); printf(" With A:"); p_mR(A, S8,P4,C7); printf(" And b:"); p_mR(b, S8,P4,C7); printf(" Ab:"); p_mR(Ab, S8,P4,C7); stop(); clrscrn(); printf(" Copy/Paste into the octave window.\n\n"); p_Octave_mR(Ab,"Ab",P4); printf("\n rref(Ab,.00000000001)\n\n"); printf(" gj_TP_mR(Ab):"); p_mR(gj_TP_mR(Ab), S10,P4,C7); printf(" You have x**4 = 16, y**3 = 27, z**2 = 25 \n\n"); printf(" These give x = 2, y = 3, z = 5 \n\n\n"); stop(); f_mR(Ab); f_mR(b); f_mR(A); return 0; } /* ------------------------------------ */ /* ------------------------------------ */ Screen output example: You have this nonlinear system of equations : a x**4 + b y**3 + c z**2 = d e x**4 + f y**3 + g z**2 = h i x**4 + j y**3 + k z**2 = l We write X = x**4; Y = y**3; Z = z**2; a X + b Y + c Z = d e X + f Y + g Z = h i X + j Y + k Z = l Ab : a b c d e f g h i j k l Press return to continue. With A: +4.0000 -2.0000 -7.0000 +9.0000 -4.0000 -2.0000 +1.0000 +6.0000 +9.0000 And b: -165.0000 -14.0000 +403.0000 Ab: +4.0000 -2.0000 -7.0000 -165.0000 +9.0000 -4.0000 -2.0000 -14.0000 +1.0000 +6.0000 +9.0000 +403.0000 Press return to continue. Copy/Paste into the octave window. Ab=[ +4.0000,-2.0000,-7.0000,-165.0000; +9.0000,-4.0000,-2.0000,-14.0000; +1.0000,+6.0000,+9.0000,+403.0000] rref(Ab,.00000000001) gj_TP_mR(Ab): +1.0000 +0.0000 +0.0000 +16.0000 +0.0000 +1.0000 +0.0000 +27.0000 +0.0000 +0.0000 +1.0000 +25.0000 You have x**4 = 16, y**3 = 27, z**2 = 25 These give x = 2, y = 3, z = 5 Press return to continue.