Как интегрировать без ограничений в SymPy?

Question

Мне дают p(x) = 2ax^2 - b и просят интегрировать его дважды, но мне не дают ограничений.Однако, когда я просто использую integrate, мне не дают константу интеграции, что крайне важно для следующего шага.

Поскольку у меня нет ограничений, использование scipy.integrate.dblquad не работает,Вместо этого я импортирую integrate из SymPy и делаю два отдельных интеграла.Вот что у меня есть:

from sympy import integrate

x = Symbol('x')
a = 240
b = 160

f = 2*a*x**2 - b

g = integrate(f) 

h = integrate(g)

Что я хочу получить:

g = 160*x**3 - 80*x**2 + C
h = 40*x**4 - 80*x**2 + Cx + c

, где C и c - это константы интеграции.Мой код в настоящее время создает уравнения без констант.

Answer 1

Из документов:

Обратите внимание, что SymPy не включает константу интеграции.Если вы хотите, вы можете добавить его самостоятельно или перефразировать вашу проблему как дифференциальное уравнение и использовать dsolve для ее решения, которое добавляет константу (см. Решение дифференциальных уравнений ).

Похоже, что они описали это в документации.Нужно сами добавить константу или переформулировать вопрос и задать другой решатель!