sympy 통합에서 발생한 예외 처리 방법

Question

sympy를 배우기 시작했습니다. sympy를 사용하여 하나의 통합 문제를 시도했지만이 예외가 발생합니다.

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를 배우기 시작했다. 내가 확인을 위 sympy에서 예외를 잡기하고있는 방법이 있나요 :이

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()) 

예상 역도는

-((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)) 

내 질문입니다 내가

를 사용하는 코드는? 예외가 무엇이 될지 모르겠다. 그리고 sympy가 적분을 해결할 수 없다면 예외가 던져지기를 기대하는 것입니까?

저는 Linux에 최신 Anaconda 64bit를 설치했습니다.

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

Answer 1

예외는 SymPy에서 버그로 간주됩니다.

표현식에있는 5/2이 파이썬에 의해 부동 소수점 숫자로 평가된다는 점에 유의하십시오. 그러면 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는 현재 어쨌든 정수를 계산할 수 없습니다.