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

Question

我有20个方程式，从Hodgkin Huxley model到cressman model，我想分析一下我通过数值方法得到的模型。

我是一个数学家，对编码了解不多，所以如果有人能用runge kutta IV method帮助我用scilab编写这些方程式，我将非常感激。

如果你感兴趣，我能做的就是分享我在神经元activiy上的工作。

致以问候。

Answer 1

除了用户ODE函数之外，您可以在不编写任何程序的情况下直接开箱即用，

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

其中(t0,y0)是初始点，times是在y中报告解决方案样本的次数列表，function dydt = odefunc(t, y)是实现导数的函数。

因为它背后有一个经过验证的健壮算法，所以您可以假设实际的数值误差与公差rtol, atol有意义地相关。

有关示例，请参阅内置文档和在线示例，如Solving ODE in Scilab、Using scilab to solve and plot differential equations