Linear Algebra and the C Language/a0hh

/* ------------------------------------ */
/*  Save as :   c00c.c                  */
/* ------------------------------------ */
#include "v_a.h"
/* ------------------------------------ */
#define   RC  RC6
/* ------------------------------------ */
int main(void)
{
/* Toeplitz Matrix 
   
        V
 
   U    1 5 6 7
        2 
        3         
        4
     
  */
  
double u[R1*RC]={ 1,20,0,0,0,0};
double v[RC*C1]={ 1,
	              30,
	              0,
	              0,
	              0,
	              0};
	              	                 
double **V   = ca_A_mR(v,i_mR(RC,C1));	
double **U   = ca_A_mR(u,i_mR(R1,RC));
double **A   =           i_mR(RC,RC);

  clrscrn();
  
  rToeplitz_mR(U,V,A);
  
  printf(" A: Tridiagonal matrix");
  p_mR(A, S4,P0,C10);
  stop();

  f_mR(U);
  f_mR(V);   
  f_mR(A);
  
  return 0;
}
/* ------------------------------------ */
/* ------------------------------------ */

                                                                                       
 A: Tridiagonal matrix
  +1  +20   +0   +0   +0   +0 
 +30   +1  +20   +0   +0   +0 
  +0  +30   +1  +20   +0   +0 
  +0   +0  +30   +1  +20   +0 
  +0   +0   +0  +30   +1  +20 
  +0   +0   +0   +0  +30   +1 

 Press return to continue.