Mathc matrices/a124

Application

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

	c00a.c

/* ------------------------------------ */
/*  Save as :   c00a.c                  */
/* ------------------------------------ */
#include "v_a.h"
/* ------------------------------------ */
#define   RA R3
#define   CA C3
#define   Cb C1
/* ------------------------------------ */
int main(void)
{
double ab[RA*(CA+Cb)]={
  +4.000,   -2.000,   -7.000, -165.000, 
  +9.000,   -4.000,   -2.000,  -14.000, 
  +1.000,   +6.000,   +9.000, +403.000      
};

double **Ab = ca_A_mR(ab,i_Abr_Ac_bc_mR(RA,CA,Cb));
double **A  = c_Ab_A_mR(Ab,i_mR(RA,CA));
double **b  = c_Ab_b_mR(Ab,i_mR(RA,Cb));

  clrscrn();
  printf(" You have this nonlinear system of equations :\n\n");
  printf(" a x**4 + b y**3 + c z**2 = d \n");
  printf(" e x**4 + f y**3 + g z**2 = h \n");
  printf(" i x**4 + j y**3 + k z**2 = l \n\n");
  printf(" We write X = x**4; Y = y**3; Z = z**2;\n\n"); 
  printf(" a X + b Y + c Z = d\n");
  printf(" e X + f Y + g Z = h\n");
  printf(" i X + j Y + k Z = l\n\n");  
  printf(" Ab :\n");  
  printf(" a  b  c  d\n");
  printf(" e  f  g  h\n");
  printf(" i  j  k  l\n\n\n"); 
  stop();  
  
  clrscrn();
  printf(" With  A :");
  p_mR(A,S8,P4,C7);
  printf(" And b :");
  p_mR(b,S8,P4,C7);
  printf(" Ab :");
  p_mR(Ab,S8,P4,C7);
  getchar();

  clrscrn();
  printf(" Copy/Past into the octave window.\n\n");
  p_Octave_mR(Ab,"Ab",P4);
  printf("\n rref(Ab,.00000000001)\n\n");
  printf("  gj_TP_mR(Ab) :");
  gj_TP_mR(Ab);
  p_mR(Ab,S10,P4,C7);
  printf(" You have  x**4 = 16, y**3 = 27, z**2 = 25 \n\n");
  printf(" These give   x = 2,     y =  3,    z = 5 \n\n\n");
  stop();

  f_mR(Ab);
  f_mR(b);
  f_mR(A);
  
  return 0;
}
/* ------------------------------------ */
/* ------------------------------------ */

Exemple de sortie écran :

 ------------------------------------                                                  
 You have this nonlinear system of equations :

 a x**4 + b y**3 + c z**2 = d 
 e x**4 + f y**3 + g z**2 = h 
 i x**4 + j y**3 + k z**2 = l 

 We write X = x**4; Y = y**3; Z = z**2;

 a X + b Y + c Z = d
 e X + f Y + g Z = h
 i X + j Y + k Z = l

 Ab :
 a  b  c  d
 e  f  g  h
 i  j  k  l


 Press return to continue. 


 ------------------------------------ 
 With  A :
 +4.0000  -2.0000  -7.0000 
 +9.0000  -4.0000  -2.0000 
 +1.0000  +6.0000  +9.0000 

 And b :
-165.0000 
-14.0000 
+403.0000 

 Ab :
 +4.0000  -2.0000  -7.0000 -165.0000 
 +9.0000  -4.0000  -2.0000 -14.0000 
 +1.0000  +6.0000  +9.0000 +403.0000 


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

 Ab=[
+4.0000,-2.0000,-7.0000,-165.0000;
+9.0000,-4.0000,-2.0000,-14.0000;
+1.0000,+6.0000,+9.0000,+403.0000]


 rref(Ab,.00000000001)

  gj_TP_mR(Ab) :
   +1.0000    +0.0000    +0.0000   +16.0000 
   +0.0000    +1.0000    +0.0000   +27.0000 
   +0.0000    +0.0000    +1.0000   +25.0000 

 You have  x**4 = 16, y**3 = 27, z**2 = 25 

 These give   x = 2,     y =  3,    z = 5 


 Press return to continue.