Ask Your Question

iris's profile - activity

2017-08-01 05:04:36 -0600 received badge  Popular Question (source)
2016-10-16 09:29:00 -0600 received badge  Editor (source)
2016-10-16 09:27:01 -0600 commented answer n-th power of matrices

Thank you for answering. However, I think you misinterpreted my question. To be clear, in my question n is a variable (as stated), while the matrix I want to calculate the n-th power of is known. I reformulated my question.

2016-10-16 07:28:01 -0600 asked a question n-th power of matrices

Is there any way to calculate the n-th power of a (upper unitriangular) matrix in Sage? Here n is an integer variable.

For example, if y=matrix([[1,0,0],[0,1,1],[0,0,1]]), then I want to obtain a formula in n for y^n. In this case this would be y^n=matrix([[1,0,0],[0,1,n],[0,0,1]]).

I tried the following: sage: y=matrix([[1,0,0],[0,1,1],[0,0,1]]); var('n'); sage: y^n

This resulted in the error: "NotImplementedError: non-integral exponents not supported"

Adding sage: assume(n, 'integer') has no effect at all.