Mathc complexes/c101h

Application

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

	c00a.c

/* ------------------------------------ */
/*  Save as :   c00a.c                  */
/* ------------------------------------ */
#include "w_a.h"
/* ------------------------------------ */
void fun(int r)
{
double **Ab = r_singular_Ab_mZ(    i_Abr_Ac_bc_mZ(r,r,C1),999);
double **A  =        c_Ab_A_mZ(Ab, i_mZ(r,r));
double **b  =             r_mZ(    i_mZ(r,C1),999);

  c_A_b_Ab_mZ(A,b,Ab);  	
	
  clrscrn();
  printf(" A :");
  p_mZ(A, S5,P0, S4,P0, C6);
  printf(" b :");
  p_mZ(b, S5,P0, S4,P0, C6);
  stop();

  clrscrn();
  printf(" Copy/Past into the octave window.\n\n");
  p_Octave_mZ(Ab,"Ab",P0, P0);
  printf(" rref(Ab,.00000000001)\n\n");
  
  printf(" gj_mZ(Ab) :");
  p_mZ(gj_mZ(Ab), S8,P4, S8,P4, C4);
  
  f_mZ(Ab);
  f_mZ(b);
  f_mZ(A);
}
/* ------------------------------------ */
int main(void)
{
time_t t;

  srand(time(&t));

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

} while(stop_w());

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

La matrice A a deux lignes identique mais avec un b différent. Système incompatible.



Exemple de sortie écran :

 ------------------------------------ 
 A :
  +48 +728i  +538 +498i  -941 +322i  +142 +596i 
 +102 -397i  -873 +788i  -441 +780i  -937 +234i 
  +48 +728i  +538 +498i  -941 +322i  +142 +596i 
 +102 -397i  -873 +788i  -441 +780i  -937 +234i 

 b :
 -927 -603i 
 -445 -975i 
 -331 -701i 
 +192 -273i 

 Press return to continue. 



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

 Ab=[
+48+728*i,+538+498*i,-941+322*i,+142+596*i,-927-603*i;
+102-397*i,-873+788*i,-441+780*i,-937+234*i,-445-975*i;
+48+728*i,+538+498*i,-941+322*i,+142+596*i,-331-701*i;
+102-397*i,-873+788*i,-441+780*i,-937+234*i,+192-273*i]

 rref(Ab,.00000000001)

 gj_mZ(Ab) : This(ese) row(s) are incompatible.
 row 3;   row 4;  

 This(ese) equation(s) was already defined
 with an other value(s)        

 The system is inconsistent. 
 I prefer to close the program.
 Press return to continue.