bo198214 Wrote:Yes, but how do you compute the infinite case. You can not simply do

as this does not converge. So do you compute it as described here via the Jordan form?

Hmm, I observed that I'm not good in answering correctly in speedy discussions. But may be I've got it right in this case.

If the sequence of corresponding entries of the powers of (A-I)^n/n do not converge well with finite dimension, but at least alternate in sign, I try to approximate the sums of individual entries of the n'th powers of the matrix B=(A-I) by Euler-summation via partial sums, which is sometimes an option.

If the growthrate of the entries is below a certain constant limit constant, or is a decreasing sequence (which is for instance the case in the scalar series for the logarithm) and no further exotic behaviour has to be expected, then Euler-summation of order q can give the final approximation to arbitrary precision for each individual entry for of the intended result-matrix by finitely many partial sums (where "exotic" has to be precisely specified,well).

Is it that what you meant?

Gottfried

Gottfried Helms, Kassel