Linear Algebra and the C Language/a0pg

__NOTOC__ Install and compile this file in your working directory. /* ------------------------------------ */ /* Save as: c00b.c */ /* ------------------------------------ */ #include "v_a.h" /* ------------------------------------ */ #define RCA C2 #define P 20/* Try : 10, 15, 20, 25 */ /* ------------------------------------ */ /* ------------------------------------ */ double f( double x) { return(pow(x,P)); } char feq[] = "x**P"; /* ------------------------------------ */ /* ------------------------------------ */ int main(void) { double a[RCA*RCA]={ +0,+1, -1./2.,+3./2. }; double **A = ca_A_mR(a, i_mR(RCA,RCA)); double **A_P = pow_mR(P, A, i_mR(RCA,RCA)); double **Ev = eigs_V_mR(A, i_mR(RCA,RCA)); double **InvEv = invgj_mR(Ev, i_mR(RCA,RCA)); /* Eval = InvEv A Ev */ double **T = mul_mR(InvEv, A, i_mR(RCA,RCA)); double **Eval = mul_mR(T, Ev, i_mR(RCA,RCA)); double **f_Eval = i_mR(RCA,RCA); clrscrn(); printf(" Copy/Paste into the octave windows \n"); p_Octave_mR(A,"A",P9); printf(" A^%d\n\n",P); printf(" A^%d:",P); pE_mR(A_P,S10,P4,C6); printf(" A^%d: Ev Eval**%d InvEv",P,P); f_eigs_mR(f,Eval,f_Eval); mul_mR(Ev,f_Eval,T); mul_mR(T,InvEv,A_P); pE_mR(A_P,S10,P4,C6); stop(); f_mR(A); f_mR(A_P); f_mR(Ev); f_mR(InvEv); f_mR(Eval); f_mR(f_Eval); f_mR(T); return 0; } /* ------------------------------------ */ /* ------------------------------------ */ Screen output example: Copy/Paste into the octave windows A=[ +0.000000000,+1.000000000; -0.500000000,+1.500000000] A^20 A^20: -1.0000e+00 +2.0000e+00 -1.0000e+00 +2.0000e+00 A^20: Ev Eval**20 InvEv -1.0000e+00 +2.0000e+00 -1.0000e+00 +2.0000e+00 Press return to continue.