maxima/trunk/matrixexp.patch

   1 --- a/share/linearalgebra/matrixexp.lisp
   2 +++ b/share/linearalgebra/matrixexp.lisp
   3 @@ -138,8 +138,8 @@
   4            (print `(ratvars = ,$ratvars gcd = '$gcd algebraic = ,$algebraic))
   5            (print `(ratfac = ,$ratfac))
   6            (merror "Unable to find the spectrum")))
   7 -
   8 -    (setq res ($fullratsimp (ncpower (sub (mult z ($ident n)) mat) -1) z))
   9 +
  10 +    (setq res ($fullratsimp ($invert_by_lu (sub (mult z ($ident n)) mat) '$crering) z))
  11      (setq m (length sp))
  12      (dotimes (i m)
  13        (setq zi (nth i sp))