R=matrix([[1,0,2,3],[0,1,4,5],[0,0,0,0]])
R
b=vector([0,0,0])
b
x=R.solve_right(b)
x
y=R.transpose()
y
c=vector([0,0,0,0])
z=y.solve_right(c)
z