Scilab中的简单二阶ODE

Question

我是 scilab 的新手，我正在尝试了解 ODE 如何在 SCILAB 中用于二阶 ODE。

我的方程式是： (d^2y)/(dt^2 )+2 dy/dt + y=0

我想绘制结果。 到目前为止，我已经明白我需要将方程分解为两个一阶方程。所以程序必须是这样的：

function dx=f(x, y)
    dx(1)=2*y+1;
    dx(2)=y;
endfunction

你能检查一下这是否正确，我不确定dx(1)和dx(2)是否正确。

谢谢。

Answer 1

如果 x_1 是 y 而 x_2 是 dy/dt，你宁愿像这样写你的 ode 右手边（即 dx/dt）

function dxdt=f(t, x)
    dxdt(1) =  x(2);
    dxdt(2) = -x(1)-2*x(2);
endfunction