a=matrix([[1,0,0,0],[1/4,1,0,0],[1/3,1/3,1,0],[1/2,1/2,1/2,1]])
print a
a.inverse()