R=matrix( [[1,-0.25,-0.25,0],[-0.25,1,0,-0.25],[-0.25,0,1,-0.25],[0,-0.25,-0.25,1] ])
print R
p=vector( [50,50,25,25])
x=R.solve_right(p)
print x