式 :k**2-40>0
条件 :kは整数,k>0
ほしい解:k>=7
プログラム①:symbol定義は、だめでした。
プログラム②:Unionの要素の取り出し方を教えて下さい。
プログラム③:solveの結果の要素の取り出し方を教えて下さい。
よろしくお願いします
プログラム①
from sympy import *
k = Symbol('k',integer=True,positive=True)
ineq = 'k**2-40>0'
ans1=solve_univariate_inequality(sympify(ineq), k, relational=False)
print('#',ineq,':',ans1)
#
# (以下エラー)抜粋
# line 633, in solve_univariate_inequality
# raise ValueError(filldedent('''
# ValueError:
# k**2 - 40 > 0 contains imaginary parts which cannot be made 0 for any
# value of k satisfying the inequality, leading to relations like I < 0.
プログラム②
from sympy import *
k = Symbol('k')
ineq = 'k**2-40>0'
ans1=solve_univariate_inequality(sympify(ineq), k, relational=False)
print('#',ineq,':',ans1)
# k**2-40>0 : Union(Interval.open(-oo, -2*sqrt(10)), Interval.open(2*sqrt(10), oo))
プログラム③
from sympy import *
k = Symbol('k')
ineq = 'k**2-40>0'
print('#',ineq,':',solve(ineq))
# k**2-40>0 : ((-oo < k) & (k < -2*sqrt(10))) | ((k < oo) & (2*sqrt(10) < k))