Mathc matrices/a86
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),555.,1E-3);
double **cosh_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 **cosh_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(cosh,EigsValue,cosh_EigsValue);
//A == EigsVector * EigsValue * T_EigsVector
mul_mR(EigsVector,cosh_EigsValue,T1);
mul_mR(T1,T_EigsVector,cosh_A);
printf(" cosh(A):");
p_mR(cosh_A,S10,P4,C6);
/* ================================================================== */
eigs_V_mR(cosh_A,EigsVector);
transpose_mR(EigsVector,T_EigsVector);
/* EigsValue : T_EigsVector * A * EigsVector */
mul_mR(T_EigsVector,cosh_A,T1);
mul_mR(T1,EigsVector,EigsValue);
f_eigs_mR(acosh,EigsValue,cosh_EigsValue);
//A == EigsVector * EigsValue * T_EigsVector
mul_mR(EigsVector,cosh_EigsValue,T1);
mul_mR(T1,T_EigsVector,A);
printf(" acosh(cosh(A)) = A :");
p_mR(A,S10,P4,C6);
f_mR(A);
f_mR(cosh_A);
f_mR(EigsVector);
f_mR(T_EigsVector);
f_mR(EigsValue);
f_mR(cosh_EigsValue);
f_mR(T1);
}
/* ------------------------------------ */
int main(void)
{
time_t t;
srand(time(&t));
do
{
fun();
} while(stop_w());
return 0;
}
/* ------------------------------------ */
/* ------------------------------------ */
acosh(cosh(A)) = A Exemple de sortie écran :
------------------------------------
A :
+0.6608 +0.4029 -0.7211 +0.1469 -0.4697
+0.4029 +0.2884 -0.4317 +0.1532 -0.2475
-0.7211 -0.4317 +0.8667 -0.2079 +0.4934
+0.1469 +0.1532 -0.2079 +0.6138 -0.0042
-0.4697 -0.2475 +0.4934 -0.0042 +0.4031
cosh(A):
+1.9697 +0.5938 -1.0963 +0.2932 -0.6794
+0.5938 +1.3675 -0.6728 +0.2017 -0.4106
-1.0963 -0.6728 +2.2430 -0.3465 +0.7646
+0.2932 +0.2017 -0.3465 +1.2675 -0.1705
-0.6794 -0.4106 +0.7646 -0.1705 +1.4845
acosh(cosh(A)) = A :
+0.6608 +0.4029 -0.7211 +0.1469 -0.4697
+0.4029 +0.2884 -0.4317 +0.1532 -0.2475
-0.7211 -0.4317 +0.8667 -0.2079 +0.4934
+0.1469 +0.1532 -0.2079 +0.6138 -0.0042
-0.4697 -0.2475 +0.4934 -0.0042 +0.4031
Press return to continue
Press X return to stop