Mathc complexes/08n

Application

Installer et compiler ce fichier dans votre répertoire de travail.

	c00a.c

/* ------------------------------------ */
/*  Save as :  c00a.c                   */
/* ------------------------------------ */
#include "w_a.h"
/* ------------------------------------ */
#define  ARRAY  A7 
#define  RCA    R4  
/* ------------------------------------ */
void fun(void)
{
double **A[ARRAY];
double **P[ARRAY];
double **InvP[ARRAY];

double **T = i_mZ(RCA,RCA);
int c;

  for(c=A0; c<ARRAY; c++)
     {   
         A[c] =              i_mZ(RCA,RCA); 
         P[c] =   r_mZ(      i_mZ(RCA,RCA),9); 
      InvP[c] = inv_mZ(P[c], i_mZ(RCA,RCA));       
	   } 

  rcsymmetric_mZ(A[A0],9);
  
    for(c=A0; c<(ARRAY-A1); c++)
     {   
        mul_mZ(InvP[c],A[c],T); 
        mul_mZ(T,P[c],A[c+A1]);
	   } 

  clrscrn();
  for(c=A0; c<ARRAY; c++)
     {
	  printf(" A[%d] : \t\t\t det =  ",c);
      p_Z(det_Z(A[c]),  S8,P0, S5,P0);  
      p_mZ(A[c],  S8,P1, S8,P1,C6); 
     } 

  for(c=A0; c<ARRAY; c++)
     { 
	   f_mZ(A[c]);
	   f_mZ(P[c]);		 
	   f_mZ(InvP[c]); 
	   }
	   
  f_mZ(T);       
}
/* ------------------------------------ */
int main(void)
{
time_t t;

  srand(time(&t));
  
  do{
        fun();
        
  }while(stop_w());

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

Des matrices semblable ont le même déterminant.

Exemple de sortie écran :

 A[0] : 			 det =  +908650760   -0i 
  +200.0    +0.0i    -10.0   -56.0i    -61.0  -176.0i    +39.0    -2.0i 
   -10.0   +56.0i   +190.0    +0.0i    +27.0  -102.0i    -48.0   +23.0i 
   -61.0  +176.0i    +27.0  +102.0i   +322.0    +0.0i     +5.0    +1.0i 
   +39.0    +2.0i    -48.0   -23.0i     +5.0    -1.0i   +301.0    +0.0i 

 A[1] : 			 det =  +908650760   -0i 
  +228.4  -139.4i   -160.7  +105.6i   -262.4  +269.9i   +146.4  -240.3i 
  +175.7  +388.9i   +680.4  -170.5i   +432.7  -295.5i   -499.0   +85.6i 
   -12.1  -277.0i   -201.7   +38.2i    -21.7  +205.3i    +93.4   -54.1i 
  +207.9   +67.2i    +22.3  -175.2i    -21.0  -137.7i   +125.8  +104.6i 

 A[2] : 			 det =  +908650760   -0i 
  -307.6  +128.4i   +140.5   +18.0i   -114.9   +70.3i   -840.2  +587.7i 
  -357.0   +88.2i   +444.5  +172.5i   -225.5  +172.0i   -252.4  +556.8i 
   +17.7   -51.9i   -186.2  +103.0i   +117.7  -155.3i   -178.7  +185.4i 
  +202.7    -8.3i     +5.1    +2.5i    +47.9   +20.1i   +758.5  -145.5i 

 A[3] : 			 det =  +908650760   -0i 
  +410.7  -619.0i   +858.0 -1252.3i    +18.4  -538.9i    +53.6 +1599.3i 
  -502.9   -74.5i   -157.7  -504.7i   -240.0  -496.1i   +608.5   -27.7i 
  +485.9  +154.9i   +733.2  +732.1i   +599.4  +214.4i  -1051.3   -64.2i 
   -66.7  -506.1i   +743.6  -726.8i   +278.6  -710.6i   +160.6  +909.3i 

 A[4] : 			 det =  +908650760   +0i 
  +334.4 -1437.4i   +647.8  +384.0i   +191.3  +147.8i   -561.1   -96.2i 
 +1070.1 -1304.9i   +842.0  +766.0i   -172.4  -118.9i   -723.4 -1003.7i 
 +1196.0  +254.5i   -657.6  +449.8i   +192.1  +565.7i   +215.1  +147.5i 
 -1316.1  -827.2i   +872.0  -453.4i    +46.5  -361.0i   -355.5  +105.6i 

 A[5] : 			 det =  +908650760   -0i 
 -2199.1   +42.2i   +117.5 -1027.9i  +2294.9  +748.8i  -2784.0 -1472.6i 
  +729.1  -791.3i   +859.4  -387.1i   +117.8 +1930.7i   +187.9  -556.1i 
  +212.4 -1424.6i   +874.7   +43.9i   -523.2 +1919.8i   +855.9 -1802.1i 
 +1201.7 -2222.0i   +766.8  +500.3i  -1519.5 +1630.8i  +2876.0 -1574.9i 

 A[6] : 			 det =  +908650760   -0i 
 +6095.9 -4686.4i   +727.8 -3646.8i  -6785.6 -3305.0i  -7326.2  +704.4i 
 -7512.6 +6609.6i   -507.5 +4985.7i  +9937.0 +4381.1i  +9797.4 -1010.5i 
  -526.2  -885.7i   -454.8  -173.1i   -426.6  +807.2i     -5.4  +953.4i 
 +4488.4 -1530.5i  +1226.1 -2007.9i  -3015.3 -3584.9i  -4148.8 -1106.5i 


 Press   return to continue
 Press X return to stop