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 <denis.thenail@free.fr>
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 <denis.thenail@free.fr <mailto:denis.thenail@free.fr>>


   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
   Maxima-lang-fr@lists.sourceforge.net
   <mailto:Maxima-lang-fr@lists.sourceforge.net>