Python 如何求解任意常数的常微分方程

我正在尝试使用以下代码解决以下ODE：

from scipy.integrate import odeint
import matplotlib.pyplot as plt
import math


def model(V,t):
    dVdt=a*V-b*V*math.log(V)
    return dVdt
V0=1
a=1
b=1
t = np.linspace(0,20,50)

V=odeint(model,V0,t)

plt.plot(t,V)
plt.xlabel('time')
plt.ylabel('V(t)')
plt.show

plt.plot

---

我想现在尝试解决任意V0、a和b的ODE。我有什么办法可以做到这一点吗

你所说的任意V0，a&b是指符号积分而不是数字积分吗？如果是这样的话，看看symphy。你可以手工解决这个问题，很容易看到

d/dtlog（a-b*log（V））=-b

，所以

a-b*log（V）=（a-b*log（V0））*exp（-b*t）

等等。你心里有没有一个等式，这样的符号解在哪里（那么容易或根本）不可能实现？