Простой ODE второго порядка в Scilab

Question

Я очень новичок в scilab и пытаюсь понять, как работает ODE для второго порядка ODE в SCILAB.

Мое уравнение: (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, вы бы лучше написали свою правую оду (т.е. dx / dt) следующим образом

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