python解二阶微分方程组_用Python求解一阶和二阶微分方程组

我需要解决以下微分方程组：

$\frac{dx_1}{dt} = -k_1x_1+k_2x_2-(K_R)x_1y_1$

$\frac{dx_2}{dt} = k_1x_1-k_2x_2-k_3x_2-(K_R)x_2y_2$

$\frac{dx_3}{dt} = k_3x_3$

$\frac{dy_1}{dt} = -k_1y_1+k_2y_2-(K_R)x_1y_1$

$\frac{dy_2}{dt} = k_1y_1-k_2y_2-k_3y_2-(K_R)x_2y_2$

$\frac{dy_3}{dt} = k_3y_3$

$\frac{dz_1}{dt} = -k_1z_1+k_2z_2+(K_R)x_1y_1$

$\frac{dz_2}{dt} = k_1z_1-k_2z_2-k_3z_2+(K_R)x_2y_2$

$\frac{dz_3}{dt} = k_3z_3$

在t = 0时的初始条件为x2 =1。在时间t = 1时，将化合物y引入y2隔室中，y2 =10。KR的值为1e-3。

我已经解决了使用矩阵求幂的简单得多的系统，并且想知道是否可以使用类似的方法来解决上述系统。

我有一个隔间模型系统X，其简化版本如下所示：

然后，微分方程组为：