from sympy import symbols, Eq, solve
x, y = symbols('x y')
eq1 = Eq(2*x  + y - 1)  # i.e.  2x+y−1=0  is eq 1
eq2 = Eq(x + y - 5)    # i.e.  x+y−5=0   is eq 2
sol = solve((eq1, eq2),(x, y))
print(sol)
''' 
Output:
 SymPyDeprecationWarning(
{x: -4, y: 9}
'''