/* ------------------------------------ */
/*  Save as :   c05n.c                   */
/* ------------------------------------ */
#include "v_a.h"
/* ------------------------------------ */
void fun(int r,int rn)
{
double **A      = r_mR(i_mR(r+rn,r),9.);
double **Q      =      i_mR(r+rn,r);
double **R      =      i_mR(r,r);

  clrscrn();
  printf(" A :");
  p_mR(A,3,3,6);
  
  QRuv_mR(A,Q,R);
    
  printf(" Q :");
  p_mR(Q,3,5,6);
  
  printf(" R :");
  p_mR(R,10,5,6);
  stop();
  
  clrscrn();
  printf(" A = Q * R  :");  
  mul_mR(Q,R, A);
  p_mR(A,3,3,6);  
  
  f_mR(A);
  f_mR(Q);
  f_mR(R);
}
/* ------------------------------------ */
int main(void)
{
time_t t;

  srand(time(&t));

do
{
  fun(rp_I(2)+R2,rp_I(3));

} while(stop_w());

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