Re: [Maxima-discuss] matrixexp

[Robert D. <rob...@gm...>]

On 2014-02-20, Raymond Rogers <ray...@gm...> wrote:

> load("simplify_sum")$load(linearalgebra)$
> H_Pascal : matrix([0,0,0,0,0],[1,0,0,0,0],[0,2,0,0,0],[0,0,3,0,0],[0,0,0,4,0]);
> H_Lag : matrix([0,0,0,0,0],[-1,1,0,0,0],[-1,0,1,0,0],[-1,0,0,1,0],[-1,0,0,0,1]);
> P_Pascal : matrixexp(H_Pascal);
> P_Lag : matrixexp(H_Lag);
> P_Lag_1 : sum((H_Lag^i)/(i!),i,0,4);

Raymond, don't you mean to write H_Lag^^i (matrix power) in the
summation? I find that when I write it that way, I get a result
very close to what matrixexp returns.

best,

Robert Dodier