在 scilab 中编写 runge kutta 4 数值方法来求解描述神经元膜活动的 cressman 模型中的方程组

Question

我有从 Hodgkin Huxley 模型到 Cressman 模型的二十个方程，我想分析通过数值方法获得的模型。

有人可以帮助我使用 runge kutta IV 方法通过 scilab 编写这些方程吗？

Answer 1

除了用户 ODE 函数之外，您无需编程即可使用，直接开箱即用，

y = ode(y0, t0, times, rtol, atol, odefunc)

其中

(t0,y0)

是初始点，

times

是在

中报告解决方案样本的次数列表，

function dydt = odefunc(t, y)

是实现导数的函数。

由于其背后有经过验证且稳健的算法，因此您可以假设实际的数值误差与公差以有意义的方式相关

rtol, atol

。