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!
