Aller au contenu

Mathc complexes/c081a

Un livre de Wikilivres.



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


inv_z.c
/* ------------------------------------ */
/*  Save as :   inv_z.c                 */
/* ------------------------------------ */
#include "w_a.h"
/* ------------------------------------ */
void fun(int r)
{
double **A    = r_mZ(i_mZ(r,r),9.);
double **Inv  =      i_mZ(r,r);

  clrscrn();
  printf(" Copy/Past into the octave window.\n\n");
  p_Octave_mZ(A,"A",P0 ,P0);
  printf(" inv(A)\n\n");

  
  printf(" inv_mZ(A,Inv); :");
  pE_mZ(inv_mZ(A,Inv), S1,P4, S1,P4, C2);

  f_mZ(Inv);
  f_mZ(A);
}
/* ------------------------------------ */
int main(void)
{
time_t t;

  srand(time(&t));

do
{
   fun(rp_I(RC3)+RC1);

} while(stop_w());

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


Vérifier avec octave.



Exemple de sortie écran :
 Copy/Past into the octave window.

 A=[
-7-1*i,-9-3*i,-5+6*i,+6+8*i;
-3+4*i,+8-3*i,-1+6*i,+4-3*i;
+4+4*i,-7-1*i,-9-9*i,+2-5*i;
+8+8*i,-3-7*i,-3+2*i,+4-7*i]

 inv(A)

 inv_mZ(A,Inv); :
-3.2154e-02+2.2574e-02i -1.0356e-01+1.1025e-01i -9.3497e-02+1.5867e-02i 
-3.8531e-02+2.4073e-02i +2.4657e-02+5.6872e-02i -3.5626e-02+1.1407e-02i 
-9.6887e-03-2.2186e-02i -2.2642e-02-1.2822e-02i -3.9658e-02+3.7585e-02i 
+2.3037e-02-1.6958e-03i +4.7620e-02+1.5746e-01i -2.4740e-02+1.0619e-01i 

+1.5634e-01-8.6736e-03i 
+3.0915e-02+6.4401e-03i 
+2.4705e-03-1.5545e-02i 
+8.9167e-02-9.3258e-02i 


 Press return to continue
 Press X      to stop