sympyのsolve関数で求めた複数の不等式の解の共通範囲を求めたい．

前提・実現したいこと

python3でプログラミングをしています．sympyモジュールをインポートしてsolve関数を使って複数の不等式の解を出力しました．出力した複数の解の共通範囲を求めて出力したいと考えています．具体的にはソースコードのSolved1変数とSolved2変数の共通範囲を求めたいです．

発生している問題・エラーメッセージ

不等式の解はOr型の変数です．しかしOr型の変数の値を取り出したり，書き換えたりするためのメソッドや関数が調べても全く分かりませんでした．

該当のソースコード

python3
1import sympy as sp
2x = sp.Symbol('x')
3Solved1 = sp.solve([x ** 2 + x - 6 > 0], [x])
4Solved2 = sp.solve([x ** 2 + x - 12 < 0], [x])

試したこと

調べ方が悪かったのか，インターネットで検索してもいい情報が見つかりませんでした．

補足情報（FW/ツールのバージョンなど）

取り出したり，書き換える方法がなかったらsolve関数を使わずに一から不等式を解くプログラムを作らないといけないと考えていて少し億劫です．助言していただけると嬉しいです．よろしくお願いします．

Jupyter Notebook 6.0.1

行動規範の内容に同意します

回答1件

ベストアンサー

「実数区間　積集合」で検索したところ、sympyで実数区間の積集合を計算するには、以下の二通りがありそうです。（他にも在るかも

Intervalで区間の型にして、&で積集合をとる
座標にしてintersectionを使う（試してません）

1で書いてみたのが下記です。なんか無理やり感があるので、参考程度ということで…

python
1import sympy as sp
2x = sp.Symbol('x')
3sol1 = sp.solve([x ** 2 + x - 6 ], [x])
4sol2 = sp.solve([x ** 2 + x - 12 ], [x])
5# 不等式をやめて、以下の処理で区間に変換しています
6
7interval1_1 = sp.Interval(-sp.oo,sol1[0][0])
8interval1_2 = sp.Interval(sol1[1][0],sp.oo)
9print(interval1_1,interval1_2)
10# Interval(-oo, -3) Interval(2, oo)
11
12interval2 = sp.Interval(sol2[0][0], sol2[1][0])
13print(interval2)
14# Interval(-4, 3)
15
16print(interval1_1&interval2, interval1_2&interval2)
17# Interval(-4, -3) Interval(2, 3)