a1=matrix([[0,0,0,1],[0,0,2,0],[0,3,0,0],[4,0,0,0]])
print a1
a1.inverse()
a2=matrix([[1,0,0,0],[-1/2,1,0,0],[0,-2/3,1,0],[0,0,-3/4,1]])
a2.inverse()
var('a,b,c,d')
a3=matrix([[a,b,0,0],[c,d,0,0],[0,0,a,b],[0,0,c,d]])
a3.inverse()

Strang-1.6_4-U (last edited 2010-12-17 12:33:52 by sarada)