A=matrix([[1,-1],[3,6]])
print A
B=vector([0,18])
x=A.solve_right(B)
print x
var ('x,y')
p=implicit_plot(x-y==0,(x,-2,2),(y,-2,2))
show(p)