/*given A=[[1,2],[3,6]], m=n=2, r=1*/

A=matrix([[1,2],[3,6]])
A.columns()
m=A.nrows()
n=A.ncols()
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 "B:",k
C=(A*A.transpose())
l=C.det()
print "C:",l
print "So for given matrix one sided inverses doesn't exist \n A =\n",A

Strang-2.4-1 (last edited 2010-12-17 10:40:32 by sanne_abhi)