我正在尝试解决两个一阶ODE的系统，其中除了CA和T之外的所有值（即q，CAi，V， k，w，p，H，U）是常量，CA和T是时间的函数。方程式如下：
dCA/dt = q*(CAi - CA)/V - k*CA
dT/dt = w*(Ti - T)/(V*p) + H*k*CA/(V*p*C) + U*(Tc - T)/(V*p*C)
python中解决方案的任何想法？我看着使用scipy，但我发现的所有例子都非常深入和复杂，并且适用于二阶ODE系统。