>>> system = Matrix(( (1, 4, 2), (-2, 1, 14))) >>> solve_linear_system(system, x, y) {x: -6, y: 2}
def express(a, b, name): sym = symbols(name) sol = solve(a-sym, b) assert len(sol) == 1 return (sym, sol[0])