Mathc matrices/Fichiers c : test02d

Application

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

	c01a.c

/* ------------------------------------ */
/*  Save as :   c00b.c                  */
/* ------------------------------------ */
#include "v_a.h"
/* ------------------------------------ */
#define    RCn            RC3
/* ------------------------------------ */
void fun(void)
{
double a[R3*C3] ={ 1,2,3,
                   5,4,1,
                   7,2,5};
                        
double **A   = ca_A_mR(a,i_mR(RCn,RCn));
double **Inv =           i_mR(RCn,RCn);
double **ID  =           i_mR(RCn,RCn);

  clrscrn();
  printf(" Copy/Past into the octave window.\n\n");
  p_Octave_mR(A,"A",P0);
  printf(" inv(A)\n");
  pE_mR(inv_mR(A,Inv), S12,P4, C4);
  stop();
  
  clrscrn();
  printf(" inv(A)*A\n");  
  p_mR(mul_mR(Inv,A,ID), S8,P4, C8);
  printf(" A*inv(A)\n");  
  p_mR(mul_mR(A,Inv,ID), S8,P4, C8);
  stop();

  f_mR(ID);
  f_mR(Inv);
  f_mR(A);
}
/* ------------------------------------ */
int main(void)
{
  fun();

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

Nous résolvons un système d'équations, et nous vérifions les calculs.



Exemple de sortie écran :

 --------------------------------------
 Copy/Past into the octave window.

 A=[
+1,+2,+3;
+5,+4,+1;
+7,+2,+5]

 inv(A)

 -2.5000e-01  +5.5556e-02  +1.3889e-01 
 +2.5000e-01  +2.2222e-01  -1.9444e-01 
 +2.5000e-01  -1.6667e-01  +8.3333e-02 

 Press return to continue. 


 --------------------------------------
 inv(A)*A

 +1.0000  -0.0000  +0.0000 
 +0.0000  +1.0000  +0.0000 
 -0.0000  +0.0000  +1.0000 

 A*inv(A)

 +1.0000  +0.0000  +0.0000 
 +0.0000  +1.0000  +0.0000 
 -0.0000  -0.0000  +1.0000 

 Press return to continue.