sympy で solveset の連立方程式 の interval の 使い方を教えて下さい。
SymPy の solveset は基本的に1変数の方程式を解く関数 であり、定義域として domain=S.Reals や Interval を指定できます。
連立方程式を解く場合は nonlinsolve または solve を使用しますが、これらは solveset のように直接 domain=Interval(...) を引数として受け取りません。そのため、連立方程式の解を指定した Interval（範囲）内に絞り込みたい場合は、「一度解を求めてから Interval に含まれるかフィルタリングする」 または ConditionSet / 論理積を使う のが一般的です。
(以下、見出しのみ)
方法1：解を求めてから Interval でフィルタリングする（最も実用的）...
方法2：1変数ずつ solveset と Interval を組み合わせる...
方法3：ConditionSet と ImageSet（周期関数や無限解を扱う場合）...