Martin Rubey wrote:
> I tried to demonstrate Cayley Hamilton in Sage, but failed.  Here is what I
> tries:
> 

Alternatively, you could do everything in polynomial rings:

sage: R.<x1,x2,x3,x4>=QQ[]
sage: m=matrix(R, 2, [x1,x2,x3,x4])
sage: m

[x1 x2]
[x3 x4]
sage: m.charpoly()
x^2 + (-x1 - x4)*x - x2*x3 + x1*x4
sage: f=m.charpoly()
sage: f(m)

[0 0]
[0 0]


Jason


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

Reply via email to