Mathc matrices/a85
Installer et compiler ces fichiers dans votre répertoire de travail.
|
c00a.c |
|---|
/* ------------------------------------ */
/* Save as : c00a.c */
/* ------------------------------------ */
#include "v_a.h"
/* ------------------------------------ */
#define RCA RC5
/* ------------------------------------ */
/* ------------------------------------ */
void fun(void)
{
double **A = rEsymmetric_mR(i_mR(RCA,RCA),455.,1E-3);
double **sin_A = i_mR(RCA,RCA);
double **EigsVector = i_mR(RCA,RCA);
double **T_EigsVector = i_mR(RCA,RCA);
double **EigsValue = i_mR(RCA,RCA);
double **sin_EigsValue = i_mR(RCA,RCA);
double **T1 = i_mR(RCA,RCA);
clrscrn();
printf(" A :");
p_mR(A,S10,P4,C6);
eigs_V_mR(A,EigsVector);
transpose_mR(EigsVector,T_EigsVector);
/* EigsValue : T_EigsVector * A * EigsVector */
mul_mR(T_EigsVector,A,T1);
mul_mR(T1,EigsVector,EigsValue);
f_eigs_mR(sin,EigsValue,sin_EigsValue);
//A == EigsVector * EigsValue * T_EigsVector
mul_mR(EigsVector,sin_EigsValue,T1);
mul_mR(T1,T_EigsVector,sin_A);
printf(" sin(A) :");
p_mR(sin_A,S10,P4,C6);
/* ================================================================== */
eigs_V_mR(sin_A,EigsVector);
transpose_mR(EigsVector,T_EigsVector);
/* EigsValue : T_EigsVector * A * EigsVector */
mul_mR(T_EigsVector,sin_A,T1);
mul_mR(T1,EigsVector,EigsValue);
f_eigs_mR(asin,EigsValue,sin_EigsValue);
//A == EigsVector * EigsValue * T_EigsVector
mul_mR(EigsVector,sin_EigsValue,T1);
mul_mR(T1,T_EigsVector,A);
printf(" asin(sin(A)) = A :");
p_mR(A,S10,P4,C6);
f_mR(A);
f_mR(sin_A);
f_mR(EigsVector);
f_mR(T_EigsVector);
f_mR(EigsValue);
f_mR(sin_EigsValue);
f_mR(T1);
}
/* ------------------------------------ */
int main(void)
{
time_t t;
srand(time(&t));
do
{
fun();
} while(stop_w());
return 0;
}
/* ------------------------------------ */
/* ------------------------------------ */
asin(sin(A)) = A : Exemple de sortie écran :
------------------------------------
A :
+0.3672 -0.1919 -0.1067 -0.0639 -0.0243
-0.1919 +0.3287 +0.1053 +0.2621 +0.0517
-0.1067 +0.1053 +0.2286 +0.2096 +0.0909
-0.0639 +0.2621 +0.2096 +0.4114 +0.1167
-0.0243 +0.0517 +0.0909 +0.1167 +0.4063
sin(A) :
+0.3478 -0.1702 -0.0922 -0.0427 -0.0155
-0.1702 +0.2987 +0.0845 +0.2292 +0.0374
-0.0922 +0.0845 +0.2126 +0.1849 +0.0783
-0.0427 +0.2292 +0.1849 +0.3720 +0.0976
-0.0155 +0.0374 +0.0783 +0.0976 +0.3891
asin(sin(A)) = A :
+0.3672 -0.1919 -0.1067 -0.0639 -0.0243
-0.1919 +0.3287 +0.1053 +0.2621 +0.0517
-0.1067 +0.1053 +0.2286 +0.2096 +0.0909
-0.0639 +0.2621 +0.2096 +0.4114 +0.1167
-0.0243 +0.0517 +0.0909 +0.1167 +0.4063
Press return to continue
Press X return to stop