如何用数值方法求解该系统的任意数量的微分方程？

问题描述

我如何解决一个由k个微分方程组成的系统，并且每个方程中都包含导数？我正在尝试使用Scipy的resolve_ivp。

所有等式均具有以下形式：

equations

如何使用任何求解器对方程组进行数值求解？使用solve_ivp，看来您应该能够独立于其他方程式编写每个方程式，在这种情况下，当我们有两个以上的方程式时，这似乎是不可能的。

解决方法

如果设置C[i]=B[i,i]，则可以将方程式变换为线性系统B*z'=A。可以解决

zdot = numpy.linalg.solve(B,A)

使得导数是一个恒定线性系统的常数解，而z的结果解是线性z(t)=z(0)+zdot*t。