Linear Algebra and the C Language/a0a7

__NOTOC__ Install and compile this file in your working directory. /* ------------------------------------ */ /* Save as : c00f.c */ /* ------------------------------------ */ #include "v_a.h" /* ------------------------------------ */ void fun(int r) { double **u = r_mR(i_mR(r,C1),9.); double **v = r_mR(i_mR(r,C1),9.); double **projuv = i_mR(r,C1); double **umnsprojuv = i_mR(r,C1); clrscrn(); printf(" u :"); p_mR(u,S3,P0,C6); printf(" v :"); p_mR(v,S3,P0,C6); proj_mR(u,v,projuv); printf(" Vector component of u along v \n\n"); printf(" projuv : [ (<u,v> / ||v||^2) * v ]"); p_mR(projuv,S3,P3,C6); stop(); clrscrn(); printf(" Vector component of (u-projuv) is orthogonal to v \n\n"); printf(" u - projuv :"); sub_mR(u,projuv,umnsprojuv); p_mR(umnsprojuv,S3,P3,C6); printf(" <v,u-projuv> = %.6f \n\n", dot_R(v,umnsprojuv)); f_mR(u); f_mR(v); f_mR(projuv); } /* ------------------------------------ */ int main(void) { time_t t; srand(time(&t)); do { fun(rp_I(R3)+R2); } while(stop_w()); return 0; } /* ------------------------------------ */ /* ------------------------------------ */ Screen output example: u : +8 -2 +3 +1 v : +2 -6 -2 +1 Vector component of u along v projuv : [ (<u,v> / ||v||^2) * v ] +1.022 -3.067 -1.022 +0.511 Press return to continue. Vector component of (u-projuv) is orthogonal to v u - projuv : +6.978 +1.067 +4.022 +0.489 <v,u-projuv> = 0.000000 Press return to continue Press X return to stop