Как решить ODE с произвольными константами

Question

Я пытаюсь решить следующий 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

Я хочу попытаться решить этот ODE сейчас для произвольных V0, a ​​и b. Есть ли способ, которым я могу сделать это?