如何在不限制SymPy整合？

Question

我给p(x) = 2ax^2 - b并要求两次整合，但我没有给限制。然而，当我只是用integrate，我没有给积分常数，这是下一步的关键。

由于我没有限制，使用scipy.integrate.dblquad没有工作。相反，我从SymPy进口integrate和做两个独立的积分。这是我有：

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来解决它，这确实增加了常数（见Solving Differential Equations）。

看来，他们做了文档中介绍这一点。人们需要添加常数自己或重构的问题，并要求不同的求解！