Mathc matrices/a261
Installer et compiler ces fichiers dans votre répertoire de travail.
|
c00a.c |
|---|
/* ------------------------------------ */
/* Save as : c00a.c */
/* ------------------------------------ */
#include "v_a.h"
/* ------------------------------------ */
#define ARRAY A7
#define RCA R4
/* ------------------------------------ */
void fun(void)
{
double **A[ARRAY];
double **P[ARRAY];
double **InvP[ARRAY];
double **T = i_mR(RCA,RCA);
int i;
for(i=A0; i<ARRAY; i++)
{
A[i] = i_mR(RCA,RCA);
P[i] = r_mR( i_mR(RCA,RCA),9);
InvP[i] = inv_mR(P[i], i_mR(RCA,RCA));
}
rsymmetric_mR(A[0],99);
for(i=A0; i<(ARRAY-C1); i++)
{
mul_mR(InvP[i],A[i],T);
mul_mR(T,P[i],A[i+C1]);
}
clrscrn();
printf(" Similar matrices have the same determinant.\n\n");
for(i=A0; i<ARRAY; i++)
{
printf(" A[%d] : \t\t\t det = %+.4f ",i,det_R(A[i]));
p_mR(A[i], S9,P2,C6);
}
for(i=A0; i<ARRAY; i++)
{
f_mR(A[i]);
f_mR(P[i]);
f_mR(InvP[i]);
}
f_mR(T);
}
/* ------------------------------------ */
int main(void)
{
time_t t;
srand(time(&t));
do{
fun();
}while(stop_w());
return 0;
}
/* ------------------------------------ */
/* ------------------------------------ */
Exemple de sortie écran :
Similar matrices have the same determinant.
A[0] : det = -79862024.0000
+85.00 -54.00 +86.00 -48.00
-54.00 +3.00 +47.00 -35.00
+86.00 +47.00 -9.00 +69.00
-48.00 -35.00 +69.00 +57.00
A[1] : det = -79862024.0000
+240.56 +44.37 -46.64 +188.50
+359.56 +113.00 -272.89 +320.00
+179.22 -16.69 -77.68 +121.75
-189.83 -53.91 +154.90 -139.87
A[2] : det = -79862024.0000
+181.29 +23.36 -312.68 +462.44
-51.05 +41.56 +64.40 -122.01
-124.44 -94.16 +373.60 -472.76
-165.11 -140.19 +294.87 -460.45
A[3] : det = -79862024.0000
-51.30 +105.44 +10.91 -85.72
-282.97 +373.18 +40.63 -179.11
+1123.47 -1961.77 -277.48 +910.44
+44.52 +90.00 +32.59 +91.60
A[4] : det = -79862024.0001
-4191.67 +10314.91 -4116.75 +211.26
-2481.48 +6061.31 -2327.34 +161.09
-1738.02 +4192.17 -1711.08 -8.87
+2135.43 -5140.07 +2041.85 -22.57
A[5] : det = -79862023.9959
-14261.93 -10194.16 +19401.56 -11357.64
+11024.67 +8802.95 -15929.35 +8829.18
+12584.22 +9632.03 -17725.96 +10039.98
+29201.64 +21783.09 -40683.17 +23320.94
A[6] : det = -79862024.0030
-2911.18 +2793.42 -310.30 -704.69
-3601.74 +2749.24 -677.99 -2829.99
-8828.40 -2372.44 -5497.38 -32413.94
+2881.63 -926.87 +1072.79 +5795.33
Press return to continue
Press X return to stop