Linear Algebra and the C Language/a06s

Install and compile this file in your working directory.

/* ------------------------------------ */
/*  Save as :   c00d.c                  */
/* ------------------------------------ */
#include "v_a.h"
/* ------------------------------------ */
void fun(int r)
{
double s =  r_I(9);

double **A  =    r_mR(    i_mR(r,r),9.);
double **sA = smul_mR(s,A,i_mR(r,r)    );

  clrscrn();
  printf(" A[r][r] : r = rows \n\n\n"
         "              ( s**r)\n"
         "    A  :      (%+.0f**%d)*det_R(A)  = %+.0f", s,r, 
                                                    pow(s,r)*det_R(A) );
  p_mR(A,  S4, P0, C5);
  printf("   sA  :             det_R(%+.0f A) = %+.0f ", s,  det_R(sA));
  p_mR(sA, S4, P0, C5);

  f_mR(sA);
  f_mR(A);
}
/* ------------------------------------ */
int main(void)
{
time_t t;

  srand(time(&t));
do
{
  fun(rp_I(RC3)+C2);

} while(stop_w());


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

Screen output example:

 A[r][r] : r = rows 


              ( s**r)
    A  :      (-8**4)*det_R(A)  = -11739136
  -5   -9   +7   -2 
  +9   -8   +8   +4 
  +8   +5   -7   +9 
  -2   +5   -5   +6 

   sA  :             det_R(-8 A) = -11739136 
 +40  +72  -56  +16 
 -72  +64  -64  -32 
 -64  -40  +56  -72 
 +16  -40  +40  -48 


 Press   return to continue
 Press X return to stop