< Mathc matrices
Installer et compiler ces fichiers dans votre répertoire de travail.
eignvect.c ' |
---|
/* ------------------------------------ */
/* Save as : eignvect.c */
/* ------------------------------------ */
#include "v_a.h"
/* ------------------------------------ */
void fun(int r)
{
double **A = rsymmetric_mR(i_mR(r,r),999.);
double **EignVector= i_mR(r,r);
double **EigsValue = i_mR(r,C1);
clrscrn();
printf(" Copy/Past into the octave windows \n\n");
p_octave_mR(A,"a",1,0);
printf(" [V, E] = eigs (a,10) \n\n");
stop();
clrscrn();
eignVectoruv_mR(A,EignVector);
printf(" EignVector:");
p_mR(EignVector,S5,P6,C6);
eigsuv_mR(A,EigsValue);
printf(" EigsValue :");
p_mR(EigsValue,S13,P6,C1);
f_mR(A);
f_mR(EignVector);
f_mR(EigsValue);
}
/* ------------------------------------ */
int main(void)
{
time_t t;
srand(time(&t));
do
{
fun(rp_I(4)+R2);
} while(stop_w());
return 0;
}
/* ------------------------------------ */
/* ------------------------------------ */
Déclaration des fichiers h.
Exemple de sortie écran :
/* ------------------------------------ */
Copy/Past into the octave windows
a=[
-445,-853,+194;
-853,-911,-921;
+194,-921,-977]
[V, E] = eigs (a,10)
Press return to continue.
/* ------------------------------------ */
EignVector:
+0.326731 +0.713242 -0.620106
+0.738677 +0.216582 +0.638317
+0.589579 -0.666616 -0.456091
EigsValue :
-2023.400085
-885.338473
+575.738557
Press return to continue
Press X to stop
Cet article est issu de Wikibooks. Le texte est sous licence Creative Commons - Attribution - Partage dans les Mêmes. Des conditions supplémentaires peuvent s'appliquer aux fichiers multimédias.