Re: maxima bug



-----Joal Heagney  wrote: -----

>So is there a way to get an exponential of a
>matrix that uses the dot multiply?

Yes, Maxima has at least three packages for computing
the matrix exponential. The one I know the most about is
in the package 'linearalgebra.'  An example:
(using Maxima 5.9.3)

(%i1) display2d : false$
(%i2) load("linearalgebra")$
(%i3)  m : matrix([1,2],[2,1])$
(%i4) matrixexp(m,z);
(%o4)
matrix([%e^-z*(%e^(4*z)+1)/2,%e^-z*(%e^(4*z)-1)/2],[%e^-z*(%e^(4*z)-1)/2,%e^-z*(%e^(4*z)+1)/2])

(%i5) matrixfun(lambda([x],x*(x+1)),m);
(%o5) matrix([6,6],[6,6])
(%i6) % - (m.m + m);            <--- Note: It's m.m, not m * m.
(%o6) matrix([0,0],[0,0])
Barton