Mathc matrices/c101a

Un livre de Wikilivres.
Aller à la navigation Aller à la recherche


Sommaire


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


Crystal Clear mimetype source c.png gj_s1.c
/* ------------------------------------ */
/*  Save as :   gj_s1.c                 */
/* ------------------------------------ */
#include "v_a.h"
/* ------------------------------------ */
void fun(int r)
{
double **Ab = r_singular_Ab_mR(i_Abr_Ac_bc_mR(r,r,C1),999);
double **A  = c_Ab_A_mR(Ab,i_mR(r,r));
double **b  = c_Ab_b_mR(Ab,i_mR(r,C1));


  clrscrn();
  printf(" A :");
  p_mR(A,S8,P0,C7);
  printf(" b :");
  p_mR(b,S8,P0,C7); 
  printf(" Ab :");
  p_mR(c_A_b_Ab_mR(A,b,Ab),S8,P0,C7);
  stop();

  clrscrn();
  printf(" Copy/Past into the octave window.\n\n");
  p_Octave_mR(Ab,"Ab",P0);
  printf("\n rref(Ab,.00000000001)\n\n");
  stop();
  
  clrscrn();
  printf(" gj_mR(Ab) :                       The system is consistent");  
  gj_PP_mR(Ab,YES);
  p_mR(Ab,S10,P5,C7);
  
  is_system_consistent_mR(Ab);

  f_mR(Ab);
  f_mR(b);
  f_mR(A);
}
/* ------------------------------------ */
int main(void)
{
time_t t;

  srand(time(&t));
do
{
  fun(rp_I(R5)+R1);


} while(stop_w());

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



Exemple de sortie écran :
 ------------------------------------ 
 A :
     -73      +68     -727     -489     +418     -331 
    -411     -801     +770      +16     -623     +792 
    +954     -627     -377     +560     +460     -391 
     -73      +68     -727     -489     +418     -331 
    -121     +818     +448     +548     -755     -809 
     -73      +68     -727     -489     +418     -331 

 b :
    +204 
    +272 
    -559 
    +204 
    -477 
    +204 

 Ab :
     -73      +68     -727     -489     +418     -331     +204 
    -411     -801     +770      +16     -623     +792     +272 
    +954     -627     -377     +560     +460     -391     -559 
     -73      +68     -727     -489     +418     -331     +204 
    -121     +818     +448     +548     -755     -809     -477 
     -73      +68     -727     -489     +418     -331     +204 

 Press return to continue. 



 Copy/Past into the octave window.

 Ab=[
-73,+68,-727,-489,+418,-331,+204;
-411,-801,+770,+16,-623,+792,+272;
+954,-627,-377,+560,+460,-391,-559;
-73,+68,-727,-489,+418,-331,+204;
-121,+818,+448,+548,-755,-809,-477;
-73,+68,-727,-489,+418,-331,+204]


 rref(Ab,.00000000001)

 Press return to continue. 



 gj_mR(Ab) :                       The system is consistent
  +1.00000   +0.00000   +0.00000   +0.00000   +1.70527   +1.88088   +0.21927 
  +0.00000   +1.00000   +0.00000   +0.00000   +0.21183   -0.06474   -0.04469 
  +0.00000   +0.00000   +1.00000   +0.00000   +0.35489   +2.01953   +0.44709 
  +0.00000   +0.00000   +0.00000   +1.00000   -1.60754   -2.61534   -1.12081 
  +0.00000   +0.00000   +0.00000   +0.00000   +0.00000   -0.00000   +0.00000 
  +0.00000   +0.00000   +0.00000   +0.00000   +0.00000   -0.00000   +0.00000 


 Press return to continue
 Press X      to stop