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