Linear Algebra and the C Language/a036

__NOTOC__ Install and compile this file in your working directory. /* ------------------------------------ */ /* Save as : c00b.c */ /* ------------------------------------ */ #include "v_a.h" /* ------------------------------------ */ void fun(int r,int c) { double **A = r_mR(i_mR(r,c),9); double **B = r_mR(i_mR(r,c),9); double **AmnsB = i_mR(r,c); clrscrn(); printf(" A : "); p_mR(A,S4,P0,C6); printf(" B : "); p_mR(B,S4,P0,C6); printf(" A - B : "); p_mR(sub_mR(A,B,AmnsB),S4,P0,C6); stop(); clrscrn(); printf(" Copy/Paste into the octave window.\n\n"); p_Octave_mR(A,"A", P0); p_Octave_mR(B,"B", P0); printf("A - B\n"); p_mR(AmnsB,S4,P0,C6); f_mR(A); f_mR(B); f_mR(AmnsB); } /* ------------------------------------ */ int main(void) { time_t t; srand(time(&t)); do fun(rp_I(R4),rp_I(C4)+C2); while(stop_w()); return 0; } /* ------------------------------------ */ /* ------------------------------------ */ Screen output example: A : -3 +2 +2 -3 +4 +6 +5 -7 +6 +7 -9 -5 +3 -7 -9 -4 B : +7 -9 +8 +6 +9 -8 +5 -4 +3 -3 +1 -5 +7 +5 -8 +6 A - B : -10 +11 -6 -9 -5 +14 +0 -3 +3 +10 -10 +0 -4 -12 -1 -10 Press return to continue. Copy/Paste into the octave window. A=[ -3,+2,+2,-3; +4,+6,+5,-7; +6,+7,-9,-5; +3,-7,-9,-4] B=[ +7,-9,+8,+6; +9,-8,+5,-4; +3,-3,+1,-5; +7,+5,-8,+6] A - B -10 +11 -6 -9 -5 +14 +0 -3 +3 +10 -10 +0 -4 -12 -1 -10 Press return to continue Press X return to stop