/* ------------------------------------ */
/*  Save as :   c01h.c                  */
/* ------------------------------------ */
#include "w_a.h"
/* ------------------------------------ */
void fun(int r)
{
double **u    = r_mZ(i_mZ(r,C1),9.);
double **n_u  =      i_mZ(r,C1);

  clrscrn();
  printf(" u :");
  p_mZ(u,S8,P2,S6,P2,C6);
            
  smul_mZ(1./normuv_Z(u),u,n_u);              
  printf("     n_u = (1/||u||) u :");
  p_mZ(n_u,S8,P2,S6,P2,C6);           

  printf(" ||n_u|| = %.6f \n\n", normuv_Z(n_u));
            
  f_mZ(u);
  f_mZ(n_u);
}
/* ------------------------------------ */
int main(void)
{
time_t t;

  srand(time(&t));

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

} while(stop_w());

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