print " given that aij = i+j and bij = (-1)^i+j "
a = matrix(2,2)
var('i')
var('j')
b = matrix(2,2)
for i in range(2):
    for j in range(2):
        a[i,j] = i+j
for i in range(1,3):
    for j in range(1,3):
        b[i-1,j-1] = (-1)**(i+j)
print a,
print " = A "
print ""
print b,
print " = B "
print "=============================="
C = a*b
print C,
print " is product AB "
print "=============================="
C = b*a
print C,
print " is product of BA "

Strang-1.4-6-U (last edited 2010-12-18 07:53:39 by ch.sandeep)