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

Question

我有20个方程式，我从Hodgkin Huxley model到cressman model工作，我想分析我通过数字方法获得的模型。

我是数学家并且对编码知之甚少，所以如果有人能帮助我使用scilab使用runge kutta IV method对这些方程进行编码，我会非常gratefull。

如果您感兴趣，我可以做的就是分享我对神经元activiy的工作。

的问候。

Answer 1

除了用户ODE功能外，您可以直接开箱即用，无需编程即可使用

<div *ngFor="let joke of jokes">
<h3>{{joke.value}}</h3>
</div>

`,
providers: [RandomService]

})
export class PocetnaComponent  { 
jokes: Joke[];

constructor(private jokesService: RandomService){
this.jokesService.getRandomJokes().subscribe(jokes => {this.jokes = jokes});
}

}
interface Joke{
 id: number;
 value: string;
}

其中y = ode(y0, t0, times, rtol, atol, odefunc) 是初始点，(t0,y0)在times中报告解决方案样本的时间列表，而y是实现衍生工具的函数

由于背后有经过验证的强大算法，您可以假设实际数值误差与公差function dydt = odefunc(t, y)有意义地相关。

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