A=matrix([[1,2],[3,6]])
print "A is\n",A
A.columns()
m=A.nrows()
print "m = No. of rows =",m
n=A.ncols()
print "n = No. of columnss =",n
m==n
r=A.rank()
r
x=matrix([[-2],[1]])
x
y=matrix([[-3],[1]])
y
r==m
r==n
B=(A.transpose()*A)
k=B.det()
print "determinant of B =",k
C=(A*A.transpose())
l=C.det()
print "determinant of C =",l
print "So one sided inverses doesnot exist for given matrix \nA =\n",A

Strang-2.4.1 (last edited 2010-12-17 06:37:36 by trnkarthik)