编写Python代码求解：捕食者-被捕食者方程组 $\left\{ \begin{matrix} {\frac{dx}{dt} = ax - bxy,} & {x(0) = 60} \\ {\frac{dy}{dt} = - cy + dxy,} & {y(0) = 30} \\ \end{matrix} \right.$ 其中 *x*(*t*) 表示第 *t* 个月时兔子的总体数量, *y*(*t*) 表示第 *t* 个月时狐狸的总体数量, 参数 *a*, *b*, *c*, *d*末知. 利用表 8.6 的 13 对观测值, 拟合上式 中的末知参数 *a*, *b*, *c*, *d*. 表 8.6种群数量观测值 $\begin{array} \hline t & 0 & 1 & 2 & 3 & 4 & 5 & 6 & 8 & 10 & 12 & 14 & 16 & 18 \\ \hline x(t) & 60 & 63 & 64 & 63 & 61 & 58 & 53 & 44 & 39 & 38 & 41 & 46 & 53 \\ y(t) & 30 & 34 & 38 & 44 & 50 & 55 & 58 & 56 & 47 & 38 & 30 & 27 & 26 \\ \hline \end{array}$

\\ {\frac{dy}{dt} = cxy - dy} \end{matrix} \right.$

抱歉，我可以回答您的问题。这是一个经典的捕食者-被捕食者方程组，可以使用Python的数值解法，如欧拉法或龙格-库塔法等来求解。您可以在互联网上搜索相关资料，了解更多关于此方程组的数值解法。