I have tracked down an error in one of my SAGE scripts down to
following behavior of SAGE:

The following code does what you'd expect:

M= Matrix([[1,0,0],[0,1,0],[0,0,1]]);
M=M*91/100; #the line we modify later
print M;
N= M.inverse(); #or: N= M^(-1)
print N;
print N*M #result, should be neutral element


If you replace the line with the command M=M*91/100; by: M=M*0.91;
then following error message shows up:
ZeroDivisionError: self is not invertible

(I find it very interesting that the command: M=M*0.9; does not
produce this error)

Also interesting is that i can supress the error by inserting: M=
M*0.91; M= M.n(20); into the line. Unfortunately
this does not work in my original programm.


How can I invert a matrix with such elements?

--~--~---------~--~----~------------~-------~--~----~
To post to this group, send email to sage-support@googlegroups.com
To unsubscribe from this group, send email to 
sage-support-unsubscr...@googlegroups.com
For more options, visit this group at 
http://groups.google.com/group/sage-support
URLs: http://www.sagemath.org
-~----------~----~----~----~------~----~------~--~---

Reply via email to