Barton Willis wrote:
>-----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
>
>
>
>
Ah. Groovy.
Joal Heagney