使用odeint求解耦合ODE组异常：咖啡豆烘焙模型问题

我正在为课程作业搭建咖啡豆烘焙的常微分方程组（ODE）模型，方程组里只有T_b和X_b是变量，其余都是常数参数。但在用scipy.integrate.odeint求解这个系统时，得到的T_b和X_b曲线居然是恒定值，这完全不符合实际烘焙过程中温度和水分变化的逻辑。

以下是我目前写的代码片段：

from scipy.integrate import odeint
import numpy as np
import matplotlib.pyplot as plt
# 编写豆温T_b的微分方程函数
def deriv(X,t):
    T_b, X_b = X...

内容的提问来源于stack exchange，提问作者L Buck