Mathc matrices/a233

Application

Installer et compiler ces fichiers dans votre répertoire de travail.

	c01a.c

/* ------------------------------------ */
/*  Save as :   c01a.c                  */
/* ------------------------------------ */
#include "v_a.h"
/* ------------------------------------ */
void fun(int r)
{
double **u     =   r_mR(      i_mR(r,C1),9.);
double **v     =   r_mR(      i_mR(r,C1),9.);
double **umnsv = sub_mR(u, v, i_mR(r,C1)   );

  clrscrn();
  printf(" u :");
  p_mR(u,S3,P0,C6);
  printf(" v :");
  p_mR(v,S3,P0,C6); 
  
  printf(" ||u||^2 + ||v||^2 - ||u - v||^2 = %.0f \n\n", 
 
                 pow(norm_R(    u),2.)
              +  pow(norm_R(    v),2.) 
              -  pow(norm_R(umnsv),2.));
      
  printf("                       2 * <u,v> = %.0f \n\n", 2 * dot_R(u,v) );

  f_mR(u);
  f_mR(v);
  f_mR(umnsv);
}
/* ------------------------------------ */
int main(void)
{
time_t t;

  srand(time(&t));

do
{  
  fun(rp_I(R3)+R2);
  
} while(stop_w());

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

Nous vérifions une des propriétés de la norme : norm_R(u);



Exemple de sortie écran :

 ---------------------
 u :
 +6 
 -2 
 +9 
 -3 
 +7 

 v :
 +1 
 -6 
 -1 
 +5 
 +7 

 ||u||^2 + ||v||^2 - ||u - v||^2 = 86 

                       2 * <u,v> = 86 


 Press   return to continue
 Press X return to stop