Linear Algebra and the C Language/a0mj

/* ------------------------------------ */
/*  Save as :   c00d.c                  */
/* ------------------------------------ */
#include "v_a.h"
/* ------------------------------------ */
void fun(int r,int rn)
{
double **A          =         r_mR(            i_mR(r+rn,r),9);
double **A_T        = transpose_mR(A,          i_mR(r,r+rn));
double **A_TA       =       mul_mR(A_T,A,      i_mR(r,r)); 

double **SvdValue_2 =      eigs_mR(A_TA,       i_mR(r,C1));
double **SvdValue   =  sqrt_svd_mR(SvdValue_2, i_mR(r,C1));

  clrscrn();
  printf(" A :");
  p_mR(A,S5,P0,C6);  
   
  printf(" Copy/Paste into the octave windows \n\n\n");
  p_Octave_mR(A,"a",P0);
  printf(" SvdValue = svd (a,10)\n\n\n");
  stop();
  
  clrscrn();
  printf(" I calculate the square root of\n"
         " each eigenvalue and obtain the singular values.\n\n"
         " sqrt_svd_mR(SvdValue**2, SvdValue);\n\n" 
         " SvdValue**2:");
  p_mR(SvdValue_2,S13,P4,C1);

  printf(" SvdValue:");
  p_mR(SvdValue,S13,P4,C1);  
        
  f_mR(A);
  f_mR(A_T);
  f_mR(A_TA);
  f_mR(SvdValue_2);
  f_mR(SvdValue);
}
/* ------------------------------------ */
int main(void)
{
time_t t;

  srand(time(&t));

do
{
  fun(rp_I(R4),rp0_I(R4));

} while(stop_w());

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

                                                                                       
 A :
   +5    +5    -3    +7 
   -1    +2    -6    -2 
   -4    -5    -6    -6 
   +2    +9    -2    +7 
   +4    -1    +9    -6 

 Copy/Paste into the octave windows 


 a=[
+5,+5,-3,+7;
-1,+2,-6,-2;
-4,-5,-6,-6;
+2,+9,-2,+7;
+4,-1,+9,-6]

 SvdValue = svd (a,10)


 Press return to continue. 


 I calculate the square root of
 each eigenvalue and obtain the singular values.

 sqrt_svd_mR(SvdValue**2, SvdValue);

 SvdValue**2:
    +316.5309 
    +183.3032 
     +25.5570 
     +12.6089 

 SvdValue:
     +17.7913 
     +13.5390 
      +5.0554 
      +3.5509 


 Press   return to continue
 Press X return to stop    

 Press X return to stop