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