All published worksheets from http://sagenb.org
var('x y') e1 = x + y == 3 e2 = x - y == 1 e3 = 3 * x - y == 5 solve([e1,e2,e3],x,y)