私はPythonでこの小さなテストプログラムを実行して、どのようsolve
に機能するかを確認しnsolve
ました。
from sympy import *
theta = Symbol('theta')
phi = Symbol('phi')
def F(theta,phi):
return sin(theta)*cos(phi)+cos(phi)**2
def G(phi):
return ((1 + sqrt(3))*sin(phi) - 4*pi*sin(2*phi)*cos(2*phi))
solution1 = solve(F(pi/2,phi),phi)
solution2 = solve(G(phi),phi)
solution3 = nsolve(G(phi),0)
solution4 = nsolve(G(phi),1)
solution5 = nsolve(G(phi),2)
solution6 = nsolve(G(phi),3)
print solution1, solution2, solution3, solution4, solution5, solution6
そして、私はこの出力を取得します:
[pi/2, pi] [] 0.0 -0.713274788952698 2.27148961717279 3.14159265358979
solve の最初の呼び出しで、対応する関数の 2 つの解が得られました。しかし、2番目のものではありません。なぜだろう?nsolve
初期テスト値で機能するようですが、その値に応じて、異なる数値解が得られます。nsolve
たった1行で、または別の関数を使用して、または別の関数を使用して、すべての数値解のリストを取得する方法はありますか?