Linear Algebra and the C Language/a0ma

/* ------------------------------------ */
/*  Save as :   c04c.c                  */
/* ------------------------------------ */
#include "v_a.h"
/* ------------------------------------ */
#define   RCA    RC4
/* ------------------------------------ */
#define    EV    R1  
/* ------------------------------------ */
int main(void)
{
double a[RCA*RCA]={
+9.825339522090, +7.309094034726, -3.954100051573, -1.198212136840, 
+7.309094034726, +5.437252879491, -2.941464672512, -0.891352931064, 
-3.954100051573, -2.941464672512, +1.591284167096, +0.482207323363, 
-1.198212136840, -0.891352931064, +0.482207323363, +0.146123431322 
};

double evalue[RCA*C1]={
+17, 
 -0,  
 -0, 
 -0      
};

double **A        =     ca_A_mR(a,                    i_mR(RCA,RCA));
double **EValue   =     ca_A_mR(evalue,               i_mR(RCA,C1));

double **Ide      =      eye_mR(                      i_mR(RCA,RCA));
double **sIde     =     smul_mR(EValue[EV][C1],Ide,   i_mR(RCA,RCA));
double **AmnssIde =      sub_mR(A,sIde,               i_mR(RCA,RCA));

double **b        =       m0_mR(i_mR(RCA,C1));
double **Ab       = c_A_b_Ab_mR(AmnssIde,b, i_Abr_Ac_bc_mR(RCA,RCA,C1));

double **Ab_free  =                         i_Abr_Ac_bc_mR(RCA,RCA,RCA);
double **b_free   =                                   i_mR(RCA,RCA);

double **V        =                                   i_mR(RCA,RCA);

  clrscrn();
  printf(" An eigenvector associated with an eigenvalue\n"
         " with single multiplicity  is unique, up to a\n"
         " change in the sign.                     \n\n\n"
         " A:");
  p_mR(A,S10,P5,C10);
  
  
        gj_PP_mR(Ab,NO);  
    put_zeroR_mR(Ab,Ab_free);
    put_freeV_mR(Ab_free); 
        gj_PP_mR(Ab_free,YES);      
       c_Ab_b_mR(Ab_free,b_free);

          c_c_mR(b_free,C2,V,EV);  
                    
    Normalize_mR(V);
           
           
          printf("\n\n\n"
                 " Copy V%d into the last file\n\n"
                 " V%d:",EV,EV);
          P_mR(V,S10,P12,C10);   
  stop();
             
  f_mR(A);
  f_mR(EValue);
  
  f_mR(Ide);
  f_mR(sIde);
  f_mR(AmnssIde);
  
  f_mR(Ab);
  f_mR(b);
  f_mR(Ab_free);    
  f_mR(b_free);
  
  f_mR(V); 

  return 0;
}
/* ------------------------------------ */
/* ------------------------------------ */

                                                                                       
 An eigenvector associated with an eigenvalue
 with single multiplicity  is unique, up to a
 change in the sign.                     


 A:
  +9.82534   +7.30909   -3.95410   -1.19821 
  +7.30909   +5.43725   -2.94146   -0.89135 
  -3.95410   -2.94146   +1.59128   +0.48221 
  -1.19821   -0.89135   +0.48221   +0.14612 




 Copy V1 into the last file

 V1:
-0.760237560476,       -nan,       -nan,       -nan, 
-0.565542575476,       -nan,       -nan,       -nan, 
+0.305949262143,       -nan,       -nan,       -nan, 
+0.092711897619,       -nan,       -nan,       -nan  

 Press return to continue.