[?] Invert a matrix and then get an approximation using Neumann Series?

[gardete (Pedro Gardete)]

Hi,

Below is a simple snippet of code performing the inverse of matrix I-M, and then an approximation using Neumann Series.

Despite the fact that all eigenvalues of M are in the unit circle, I can’t get it to match past the second decimal, independently of the number of terms I add.

My guess is the calculation precision is insufficient? Any ideas how to solve this?

Thanks!

[sghuisman (Sander Huisman)]

I think you need MatrixPower rather than the regular Power… I would try it out if you would have posted your code as text…

[gardete (Pedro Gardete)]

Genius! Thanks,
Pedro