Mathc matrices/c11b6
Installer et compiler ces fichiers dans votre répertoire de travail.
![]() |
c0b6.c |
---|
/* ------------------------------------ */
/* Save as : c0b6.c */
/* ------------------------------------ */
#include "v_a.h"
/* ------------------------------------ */
void fun(void)
{
double ab[R6*C6]={
+754,+750,-983,-721,-737,+394,
+750,+886,+930,+794,-809,-245,
-983,+930,+990,-747,-37,-647,
-721,+794,-747,+886,-301,+682,
-737,-809,-37,-301,-491,+816,
+394,-245,-647,+682,+816,-33
};
int r = RC6;
double **A = ca_A_mR(ab,i_mR(r,r));
double **EignVector= i_mR(r,r);
double **T_EignVector= i_mR(r,r);
double **EigsValue = i_mR(r,C1);
double **T1= i_mR(r,r);
double **T2= i_mR(r,r);
clrscrn();
printf(" Copy/Past into the octave windows \n\n");
p_octave_mR(A,"a",S1,P0);
printf(" [V, E] = eigs (a,10) \n\n");
stop();
clrscrn();
eignVectoruv_mR(A,EignVector);
printf(" EignVector:");
p_mR(EignVector,S10,P4,C6);
printf(" EigsValue : T_EignVector * A * EignVector");
transpose_mR(EignVector,T_EignVector);
mul_mR(T_EignVector,A,T1);
mul_mR(T1,EignVector,T2);
p_mR(T2,S10,P4,C6);
eigsuv_mR(A,EigsValue);
printf(" EigsValue :");
p_mR(EigsValue,S10,P4,C1);
stop();
f_mR(A);
f_mR(EignVector);
f_mR(T_EignVector);
f_mR(EigsValue);
f_mR(T1);
f_mR(T2);
}
/* ------------------------------------ */
int main(void)
{
fun();
return 0;
}
/* ------------------------------------ */
/* ------------------------------------ */
Déclaration des fichiers h. Exemple de sortie écran :
/* ------------------------------------ */
Copy/Past into the octave windows
a=[
+754,+750,-983,-721,-737,+394;
+750,+886,+930,+794,-809,-245;
-983,+930,+990,-747,-37,-647;
-721,+794,-747,+886,-301,+682;
-737,-809,-37,-301,-491,+816;
+394,-245,-647,+682,+816,-33]
[V, E] = eigs (a,10)
Press return to continue.
/* ------------------------------------ */
EignVector:
-0.1631 +0.5430 +0.5025 -0.5953 -0.1793 +0.1986
+0.4504 +0.6251 -0.2601 +0.1579 +0.5032 +0.2462
+0.7343 -0.2242 +0.2854 +0.0874 -0.4360 +0.3625
-0.1972 +0.3543 +0.4424 +0.7597 -0.1989 -0.1522
-0.1891 -0.3656 +0.4749 +0.0991 +0.6016 +0.4829
-0.3959 +0.0699 -0.4209 +0.1613 -0.3508 +0.7157
EigsValue : T_EignVector * A * EignVector
+2337.8319 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000
-0.0000 +2099.7550 -0.0000 +0.0000 -0.0000 -0.0000
-0.0000 -0.0000 -1853.6437 +0.0000 +0.0000 +0.0000
-0.0000 +0.0000 +0.0000 +1635.5483 -0.0000 +0.0000
-0.0000 -0.0000 +0.0000 -0.0000 -1297.4337 +0.0000
-0.0000 -0.0000 +0.0000 -0.0000 +0.0000 +69.9422
EigsValue :
+2337.8319
+2099.7550
-1853.6437
+1635.5483
-1297.4337
+69.9422
Press return to continue.