Неверный результат интеграции sympy

Question

Это выражение возвращает ноль, но не должно.
P = x^6-14x^4 + 49x^2-36
интегрировать (1/P, (x, 1/3, 1/2))
Я также использовал выражение для выражения без каких-либо результатов. Я делаю что-то неправильно или это ошибка?

Answer 1

Это работает:

from sympy import * 

x = symbols('x') 
P = x**6-14*x**4+49*x**2-36 
I = integrate(1/expand(P), (x, S.One/3, S.One/2)) 

Я получаю результат:

In [5]: I 
Out[5]: -3*log(3)/80 - log(7)/48 - log(2)/48 - log(8)/240 + log(10)/240 + log(4)/48 + 3*log(5)/80 

In [6]: I.n() 
Out[6]: -0.00601350282195297 

В альтернативной, вы можете запустить команду isympy -i, это будет работать запрос SymPy, который преобразует все Python целые числа к целям SymPy до того, как вход будет оцениваться синтаксическим анализатором SymPy.

Целочисленное деление Python отличается между Python 2 и Python 3, первым возвращает и целое число, второе возвращает число с плавающей запятой. Обе версии отличаются от целочисленного деления SymPy, которое возвращает фракции. Чтобы использовать разделение SymPy, вам нужно убедиться, что по крайней мере один из дивидендов и делителей - это объекты SymPy.