Re: [sage-support] Formal power of one matrix

2013-01-24 Thread Stephen Montgomery-Smith
See attachment. Mathematica also has a MatrixExp command, and on 2 by 2 matrices the answers are really nice. For example, if the matrix has complex eigenvalues, the answers are written in terns of Sin and Cos rather than exponentials of complex numbers. On 01/24/13 11:20, Christophe BAL wr

Re: [sage-support] Formal power of one matrix

2013-01-24 Thread Christophe BAL
I'm curious to see this output. Can you send it ? 2013/1/24 Stephen Montgomery-Smith > On 01/24/13 08:57, Christophe BAL wrote: > >> Hello, >> I would like, if it is possible, to calculate the formal power of one >> matrix ? >> >> My attempt is after but it doesn't work... :-( >> >> Christophe >

Re: [sage-support] Formal power of one matrix

2013-01-24 Thread Stephen Montgomery-Smith
On 01/24/13 08:57, Christophe BAL wrote: Hello, I would like, if it is possible, to calculate the formal power of one matrix ? My attempt is after but it doesn't work... :-( Christophe == var('n') assume(n, 'integer') E = matrix([ [0 , 1 , 0 , 0 , 0 ],

Re: [sage-support] Formal power of one matrix

2013-01-24 Thread Ivan Andrus
I don't think such a thing is possible in the way you are hoping for. My suggestion would be to find the Jordan normal form (and save the conjugating matrix by passing `transformation=True`). Your matrix is diagonalizable so you can then take a formal exponential of that (though you may have to

[sage-support] Formal power of one matrix

2013-01-24 Thread Christophe BAL
Hello, I would like, if it is possible, to calculate the formal power of one matrix ? My attempt is after but it doesn't work... :-( Christophe == var('n') assume(n, 'integer') E = matrix([ [0 , 1 , 0 , 0 , 0 ], [1/4 , 0 , 3/4 , 0 , 0 ], [0 , 1/