A=matrix([[0,0,1,0,0],[0,0,0,1,0],[0,0,0,0,1],[0,0,0,0,0],[0,0,0,0,0]])
print A
A.jordan_form(transformation=True)
B=matrix([[1,-1,0,-1],[0,2,0,1],[-2,1,-1,1],[2,-1,2,0]])
print B
B.jordan_form(transformation=True)