sympyから発生した例外を処理する方法

Question

sympyを学び始めています。私はsympyを使って1つの統合問題を試みましたが、私はこの例外を受け取ります：

File "/home/me/anaconda3/lib/python3.6/site- 
     packages/sympy/polys/domains/domain.py", 
     line 146, 
    in convert raise CoercionFailed("can't convert %s of type %s to %s" % 
    (element, type(element), self)) 

sympy.polys.polyerrors.CoercionFailed: can't convert (_x4*a + a)**0.5 of 
type 
<class 'sympy.core.power.Pow'> to RR(a,f,n,A,B) 
    [_A0,_A1,_A2,_A3,_A4,_A5,_A6,_A7,_A8,_A9,_A10,_A11, 
    _A12,_A13,_A14,_A15,_A16,_A17 
... 
raise ValueError("expected an expression convertible to a polynomial in %s, 
    got %s" % (self, expr)) 
    ValueError: expected an expression convertible to a polynomial in 
    Polynomial ring in 

完全なメッセージは非常に長いです。私は何か間違っている可能性が高いですか？私はsympyを学び始めたばかりです。期待される不定積分は

-((2*a^3*(2*B*(115 + 203*n + 104*n^2 + 16*n^3) + A* 
(301 + 478*n + 224*n^2 + 32*n^3))* 
Cos[e + f*x]*Hypergeometric2F1[1/2, -n, 3/2, 1 - Sin[e + f*x]]* 
(d*Sin[e + f*x])^n)/ 
    Sin[e + f*x]^n/(f*(3 + 2*n)*(5 + 2*n)*(7 + 2*n)* 
    Sqrt[a + a*Sin[e + f*x]])) - 
    (2*a^3*(2*B*(35 + 23*n + 4*n^2) + A*(77 + 50*n + 8*n^2))* 
    Cos[e + f*x]* 
    (d*Sin[e + f*x])^(1 + n))/(d*f*(3 + 2*n)*(5 + 2*n)*(7 + 2*n)* 
    Sqrt[a + a*Sin[e + f*x]]) - 
    (2*a^2*(2*B*(5 + n) + A*(7 + 2*n))*Cos[e + f*x]* 
    (d*Sin[e + f*x])^(1 + n)* 
    Sqrt[a + a*Sin[e + f*x]])/(d*f*(5 + 2*n)* 
    (7 + 2*n)) - (2*a*B*Cos[e + f*x]*(d*Sin[e + f*x])^(1 + n)* 
    (a + a*Sin[e + f*x])^(3/2))/(d*f*(7 + 2*n)) 

ですこれは私が

from sympy import * 
import traceback 
import logging 

logging.basicConfig(filename='sympy_err.log') 

x, A, B, f, e, c, d, a, b, n = symbols('x A B f e c d a b n', real=True) 
try: 
    integrate((d*sin(f*x+e))**n*(a+a*sin(f*x+e))**(5/2)*(A+B*sin(f*x+e)),x) 
except Exception as e: 
    logging.error(traceback.format_exc()) 

を使用したコードは、私の質問はされていますは、私はOKの上sympyからの例外をキャッチしています方法は何ですか？私は何の例外が襲撃されるのか分からないので。 sympyが積分を解くことができず、例外をスローすることができれば、これは何かが期待できるでしょうか？

私はLinuxを使用しています。ちょうどLinuxに最新のAnaconda 64bitをインストールしました。

code>conda list sympy 
# packages in environment at /home/me/anaconda3: 
# 
sympy      1.0      py36_0 

Answer 1

あなたが見ている例外は、SymPyのバグとみなされるべきです。

しかし、表現中の5/2は、Pythonによって浮動小数点数に評価されているため、SymPyにはより困難な時間があります。通常、浮動小数点の代わりに、有理数をSymPyで使用することをお勧めします。正しいバージョンはintegrate((d*sin(f*x+e))**n*(a+a*sin(f*x+e))**Rational(5, 2)*(A+B*sin(f*x+e)),x)です。しかし、現在のところ、この積分にはSymPyがハングアップしています（終了するのは面倒ではなく、答えを見つけることはできません）。

これは不足しています：あなたが見ている例外はバグであり、SymPyは現在のところ積分を計算できません。