SymPy function with 2 variables

Specific Domains Optimization (Mathematical)
2

Please, read PSA: make it easier to help you, especially, points 1 and 3.

From the documentation you should use solveset, no?