Bonjour, quel type de machine et quel OS travaillez-vous ? Pour ce qui me concerne, c'est UBUNTU 9.10 sur un PC avec pour proc un athlon 64 X2 4400+.
Merci d'avance Stavros Macrakis a écrit : > Salut, > > Je n'arrive toujours pas a reproduire le probleme, neanmoins je > remarque (ce que j'aurais du remarquer avant...) que l'erreur a ete > signalee par dimnary, donc c'est un probleme de formattage et non pas > de calcul. Sans pouvoir resoudre ce probleme, je vous conseillerais > simplement de ne pas imprimer l'expression %o6 (utiliser "$" au lieu > de ";" a la fin de la commande %i6) ou bien de l'imprimer en mode > display2d:false. > > -s > > 2009/12/16 denis.thenail <[email protected] > <mailto:[email protected]>> > > Bonjour, > > la version de MAXIMA est 5.17.1 > > j'ai construit la matrice en plusieurs étapes, c'est à dire > plusieurs scripts > > voici celui qui à partir de la matrice tridiagonale construit la > matrice à inverser : > Ka5:matrix([X,-X,0,0,0,0]); > Ka5:addrow(Ka5,[-X,2*X,-X,0,0,0]); > Ka5:addrow(Ka5,[0,-X,2*X,-X,0,0]); > Ka5:addrow(Ka5,[0,0,-X,2*X,-X,0]); > Ka5:addrow(Ka5,[0,0,0,-X,2*X,-X]); > Ka5:addrow(Ka5,[0,0,0,0,-X,X]); > Ma5:matrix([2*Y,Y,0,0,0,0]); > Ma5:addrow(Ma5,[Y,4*Y,Y,0,0,0]); > Ma5:addrow(Ma5,[0,Y,4*Y,Y,0,0]); > Ma5:addrow(Ma5,[0,0,Y,4*Y,Y, 0]); > Ma5:addrow(Ma5,[0,0,0,Y,4*Y,Y]); > Ma5:addrow(Ma5,[0,0,0,0,Y,2*Y]); > Ia5:invert(Ma5); > Aa5:Ia5 . Ka5; > > où Ka5 est la matrice de raideur et Ma5 une matrice de masse. > A part ça, oui, je sais que je peux écrire la marice en une seule > fois. > > voici le script complet > > Ka5:matrix([X,-X,0,0,0,0],[-X,2*X,-X,0,0,0],[0,-X,2*X,-X,0,0],[0,0,-X,2*X,-X,0], > [0,0,0,-X,2*X,-X],[0,0,0,0,-X,X]); > Ma5:matrix([2*Y,Y,0,0,0,0],[Y,4*Y,Y,0,0,0],[0,Y,4*Y,Y,0,0],[0,0,Y,4*Y,Y, > 0], > [0,0,0,Y,4*Y,Y],[0,0,0,0,Y,2*Y]); > Ia5:invert(Ma5); > Aa5:Ia5 . Ka5; > charpoly(Aa5,l); > factor(%); > expand(%); > > et la réponse de MAXIMA : > (%i1) batch("test.mac"); > > batching #p/home/denis/MAXIMA/test.mac > (%i2) Ka5 : matrix([X, - X, 0, 0, 0, 0], [- X, 2 X, - X, 0, 0, 0], > [0, - X, 2 X, - X, 0, 0], [0, 0, - X, 2 X, - X, 0], [0, 0, 0, - X, > 2 X, - X], > [0, 0, 0, 0, - X, X]) > [ X - X 0 0 0 0 ] > [ ] > [ - X 2 X - X 0 0 0 ] > [ ] > [ 0 - X 2 X - X 0 0 ] > (%o2) [ ] > [ 0 0 - X 2 X - X 0 ] > [ ] > [ 0 0 0 - X 2 X - X ] > [ ] > [ 0 0 0 0 - X X ] > (%i3) Ma5 : matrix([2 Y, Y, 0, 0, 0, 0], [Y, 4 Y, Y, 0, 0, 0], > [0, Y, 4 Y, Y, 0, 0], [0, 0, Y, 4 Y, Y, 0], [0, 0, 0, Y, 4 Y, Y], > [0, 0, 0, 0, Y, 2 Y]) > [ 2 Y Y 0 0 0 0 ] > [ ] > [ Y 4 Y Y 0 0 0 ] > [ ] > [ 0 Y 4 Y Y 0 0 ] > (%o3) [ ] > [ 0 0 Y 4 Y Y 0 ] > [ ] > [ 0 0 0 Y 4 Y Y ] > [ ] > [ 0 0 0 0 Y 2 Y ] > (%i4) Ia5 : invert(Ma5) > [ 362 97 26 7 2 1 ] > [ ----- - ----- ----- - ----- ----- - ----- ] > [ 627 Y 627 Y 627 Y 627 Y 627 Y 627 Y ] > [ ] > [ 97 194 52 14 4 2 ] > [ - ----- ----- - ----- ----- - ----- ----- ] > [ 627 Y 627 Y 627 Y 627 Y 627 Y 627 Y ] > [ ] > [ 26 52 182 49 14 7 ] > [ ----- - ----- ----- - ----- ----- - ----- ] > [ 627 Y 627 Y 627 Y 627 Y 627 Y 627 Y ] > (%o4) [ ] > [ 7 14 49 182 52 26 ] > [ - ----- ----- - ----- ----- - ----- ----- ] > [ 627 Y 627 Y 627 Y 627 Y 627 Y 627 Y ] > [ ] > [ 2 4 14 52 194 97 ] > [ ----- - ----- ----- - ----- ----- - ----- ] > [ 627 Y 627 Y 627 Y 627 Y 627 Y 627 Y ] > [ ] > [ 1 2 7 26 97 362 ] > [ - ----- ----- - ----- ----- - ----- ----- ] > [ 627 Y 627 Y 627 Y 627 Y 627 Y 627 Y ] > (%i5) Aa5 : Ia5 . Ka5 > [ 153 X 194 X 52 X 14 X 4 X X ] > [ ----- - ----- ----- - ----- ----- - ----- ] > [ 209 Y 209 Y 209 Y 209 Y 209 Y 209 Y ] > [ ] > [ 97 X 179 X 104 X 28 X 8 X 2 X ] > [ - ----- ----- - ----- ----- - ----- ----- ] > [ 209 Y 209 Y 209 Y 209 Y 209 Y 209 Y ] > [ ] > [ 26 X 104 X 155 X 98 X 28 X 7 X ] > [ ----- - ----- ----- - ----- ----- - ----- ] > [ 209 Y 209 Y 209 Y 209 Y 209 Y 209 Y ] > (%o5) [ ] > [ 7 X 28 X 98 X 155 X 104 X 26 X ] > [ - ----- ----- - ----- ----- - ----- ----- ] > [ 209 Y 209 Y 209 Y 209 Y 209 Y 209 Y ] > [ ] > [ 2 X 8 X 28 X 104 X 179 X 97 X ] > [ ----- - ----- ----- - ----- ----- - ----- ] > [ 209 Y 209 Y 209 Y 209 Y 209 Y 209 Y ] > [ ] > [ X 4 X 14 X 52 X 194 X 153 X ] > [ - ----- ----- - ----- ----- - ----- ----- ] > [ 209 Y 209 Y 209 Y 209 Y 209 Y 209 Y ] > (%i6) charpoly(Aa5, l) > 153 X 179 X 155 X > (%o6) (----- - l) ((----- - l) ((----- - l) > 209 Y 209 Y 209 Y > 2 > 155 X 153 X 179 X 18818 X > ((----- - l) ((----- - l) (----- - l) - --------) > 209 Y 209 Y 209 Y 2 > 43681 Y > 179 X 153 X > 2 52 X (----- - l) 2 104 X > (----- - l) > 20176 X 209 Y 5044 X 209 Y > 26 X (-------- - ----------------) 104 X (-------- - > -----------------) > 2 209 Y 2 209 Y > 43681 Y 43681 Y > + ---------------------------------- + > ------------------------------------) > 209 Y 209 Y > 179 X > 14 X (----- - l) 2 > 155 X 209 Y 5432 X > + (7 X (- (----- - l) (---------------- - --------) > 209 Y 209 Y 2 > 43681 Y > 179 X > 2 52 X (----- - l) > 20176 X 209 Y > 98 X (-------- - ----------------) > 2 209 Y > 43681 Y > - ----------------------------------))/(209 Y) > 209 Y > 153 X > 28 X (----- - l) 2 > 155 X 209 Y 1358 X > + (28 X (- (----- - l) (---------------- - --------) > 209 Y 209 Y 2 > 43681 Y > 153 X > 2 104 X (----- - l) > 5044 X 209 Y > 98 X (-------- - -----------------) > 2 209 Y > 43681 Y > - -----------------------------------))/(209 Y) > 209 Y > 2 > 153 X 179 X 18818 X > 98 X ((----- - l) (----- - l) - --------) > 209 Y 209 Y 2 > 43681 Y > + (98 X (- ----------------------------------------- > 209 Y > 179 X 153 X > 14 X (----- - l) 2 28 X (----- - l) > 2 > 209 Y 5432 X 209 Y 1358 X > 26 X (---------------- - --------) 104 X (---------------- - > --------) > 209 Y 2 209 Y > 2 > 43681 Y > 43681 Y > + ---------------------------------- + > -----------------------------------)) > 209 Y 209 Y > 179 X > 2 52 X (----- - l) > 20176 X 209 Y > 28 X (-------- - ----------------) > 2 209 Y > 155 X 43681 Y > /(209 Y)) + (2 X (- (----- - l) (---------------------------------- > 209 Y 209 Y > 179 X > 2 4 X (----- - l) > 155 X 1552 X 209 Y > - (----- - l) (-------- - ---------------)) > 209 Y 2 209 Y > 43681 Y > 179 X > 14 X (----- - l) 2 > 155 X 209 Y 5432 X > - (104 X (- (----- - l) (---------------- - --------) > 209 Y 209 Y 2 > 43681 Y > 179 X > 2 52 X (----- - l) > 20176 X 209 Y > 98 X (-------- - ----------------) > 2 209 Y > 43681 Y > - ----------------------------------))/(209 Y) > 209 Y > 179 X > 14 X (----- - l) 2 > 209 Y 5432 X > 28 X (---------------- - --------) > 209 Y 2 > 43681 Y > - (98 X (---------------------------------- > 209 Y > > 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. > (%i7) factor(%) > Erreur de segmentation > > merci pour votre attention > > > Stavros Macrakis a écrit : > > Je n'arrive pas a reproduire votre probleme. Pourriez-vous > nous envoyer une copie complete de votre interaction avec > Maxima, y compris la version? > > Par ailleurs, savez-vous que vous pouvez ecrire la matrice en > une commande, a savoir: > > Ca5:matrix([153*Z/209,-194*Z/209,52*Z/209,-14*Z/209,4*Z/209,-X/209], > > [-97*Z/209,179*Z/209,-104*Z/209,28*Z/209,-8*Z/209,2*X/209], > ... ) > > Merci, > > -s > > 2009/12/15 denis.thenail <[email protected] > <mailto:[email protected]> <mailto:[email protected] > <mailto:[email protected]>>> > > > 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] > <mailto:[email protected]> > <mailto:[email protected] > <mailto:[email protected]>> > > https://lists.sourceforge.net/lists/listinfo/maxima-lang-fr > > > > ------------------------------------------------------------------------------ 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
