/* ------------------------------------ */
/*  Save as :   c05l.c                   */
/* ------------------------------------ */
#include "v_a.h"
/* ------------------------------------ */
void fun(int r)
{
double **U      = r_mR(i_mR(r,r),9.);
double **Orth   =      i_mR(r,r);
double **Orth_T =      i_mR(r,r);

  clrscrn();
  
  printf(" U :");
  p_mR(U,3,0,6);

  orthuv_mR(U,Orth);
  
  printf(" Orth :");
  p_mR(Orth,3,5,6);
  
  printf(" OrthT * Orth  :");
  transpose_mR(Orth,Orth_T);
  mul_mR(Orth_T,Orth, U);  
  p_mR(U,15,9,6);
    
  f_mR(U);
  f_mR(Orth);
  f_mR(Orth_T);
}
/* ------------------------------------ */
int main(void)
{
time_t t;

  srand(time(&t));

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

} while(stop_w());

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