Bonjour,
je cherche le polynôme caractéristique d'une matrice (6,6). Cette
matrice est pleine mais résulte d'un problème aux valeurs propres
généralisé tridiagonal du type (K-omega^2M)(Phi)= (0)
M définie postive, K semi-définie positive
ramené sous la forme standard
(M^{-1}K-omega^2I)(Phi)=(0)
J'exécute avec Maxima :
Ca5:matrix([153*Z/209,-194*Z/209,52*Z/209,-14*Z/209,4*Z/209,-X/209]);
Ca5:addrow(Ca5,[-97*Z/209,179*Z/209,-104*Z/209,28*Z/209,-8*Z/209,2*X/209]);
Ca5:addrow(Ca5,[26*Z/209,-104*Z/209,155*Z/209,-98*Z/209,28*Z/209,-7*Z/209]);
Ca5:addrow(Ca5,[-7*Z/209,28*Z/209,-98*Z/209,155*Z/209,-104*Z/209,26*Z/209]);
Ca5:addrow(Ca5,[2*Z/209,-8*Z/209,28*Z/209,-104*Z/209,179*Z/209,-97*Z/209]);
Ca5:addrow(Ca5,[-Z/209,4*Z/209,-14*Z/209,52*Z/209,-194*Z/209,153*Z/209]);
charpoly(Ca5,l);
factor(%);
expand(%);
la commande charpoly entraîne :
Universal error handler called recursively (:ERROR NIL MAXIMA::DIMNARY
"" "Couldn't protect")
Universal error handler called recursively (:ERROR NIL MAXIMA::DIMNARY
"" "Couldn't protect")
Maxima encountered a Lisp error:
Error in MAXIMA::DIMNARY [or a callee]: Caught fatal error [memory may
be damaged]
Automatically continuing.
To reenable the Lisp debugger set *debugger-hook* to nil.
Merci de m'aider à améliorer ma pratique de MAXIMA,
cordialement
Dude007
------------------------------------------------------------------------------
This SF.Net email is sponsored by the Verizon Developer Community
Take advantage of Verizon's best-in-class app development support
A streamlined, 14 day to market process makes app distribution fast and easy
Join now and get one step closer to millions of Verizon customers
http://p.sf.net/sfu/verizon-dev2dev
_______________________________________________
Maxima-lang-fr mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/maxima-lang-fr
