Mathc complexes/a197
Apparence
Installer et compiler ces fichiers dans votre répertoire de travail.
c01a.c |
---|
/* ------------------------------------ */
/* Save as : c01a.c */
/* ------------------------------------ */
#include "w_a.h"
/* ------------------------------------ */
#define RA R3
#define CA C3
#define Cb C1
/* ------------------------------------ */
int main(void)
{
double a[RA*(CA*C2)] ={ 1,2, 3,4, 5,6,
5,4, 1,3, 6,8,
7,2, 5,1, 1,1};
double b[RA*(Cb*C2)] ={ 1,2,
3,4,
5,6 };
double **A = ca_A_mZ(a,i_mZ(RA,CA));
double **B = ca_A_mZ(b,i_mZ(RA,Cb));
double **Inv = i_mZ(RA,CA);
double **X = i_mZ(RA,Cb);
double **T = i_mZ(RA,Cb);
clrscrn();
printf(" We want to find X such as, \n\n");
printf(" AX = B \n\n");
printf(" If A is a square matrix and, \n\n");
printf(" If A has an inverse matrix, \n\n");
printf(" you can find X by this method\n\n");
printf(" X = inv(A) B \n\n\n");
printf(" To verify the result you can \n\n");
printf(" multiply the matrix A by X. \n\n");
printf(" You must refind B. \n\n");
getchar();
clrscrn();
printf(" A :\n");
p_mZ(A, S5,P0, S4,P0, C6);
printf(" B :\n");
p_mZ(B, S5,P0, S4,P0, C6);
getchar();
clrscrn();
printf(" inv(A) :\n");
pE_mZ(inv_mZ(A,Inv), S12,P4, S12,P4, C4);
printf(" X = inv(A) * B :\n");
p_mZ(mul_mZ(Inv,B,X), S10,P4, S8,P4, C4);
getchar();
clrscrn();
printf(" B :\n");
p_mZ(B, S5,P0, S4,P0, C6);
printf(" AX :\n");
p_mZ(mul_mZ(A,X,T), S5,P0, S4,P0, C6);
f_mZ(T);
f_mZ(X);
f_mZ(B);
f_mZ(Inv);
f_mZ(A);
stop();
return 0;
}
/* ------------------------------------ */
/* ------------------------------------ */
Exemple de sortie écran :
----------------------------------
We want to find X such as,
AX = B
If A is a square matrix and,
If A has an inverse matrix,
you can find X by this method
X = inv(A) B
To verify the result you can
multiply the matrix A by X.
You must refind B.
----------------------------------
A :
+1 +2i +3 +4i +5 +6i
+5 +4i +1 +3i +6 +8i
+7 +2i +5 +1i +1 +1i
B :
+1 +2i
+3 +4i
+5 +6i
----------------------------------
inv(A) :
-1.2917e-01 +8.1007e-02i +8.5597e-02 -6.6602e-02i +8.5160e-02 -3.3778e-04i
+1.8258e-01 -8.6352e-02i -1.5122e-01 +8.1911e-02i +9.3058e-02 -3.8517e-02i
+3.7156e-03 -6.3244e-02i +7.1142e-02 -3.4463e-02i -8.4514e-02 +1.6084e-02i
X = inv(A) * B :
+0.6598 +0.4745i
+0.2704 +0.2854i
-0.0376 -0.3013i
----------------------------------
B :
+1 +2i
+3 +4i
+5 +6i
AX :
+1 +2i
+3 +4i
+5 +6i
Press return to continue.