Вот простой ODE, который нужно решить: dq / dt = L q с L=[-1/400 0; 1 -3/400]
, а q - вектор 2 * 1.
Программа SCILAB для решения этой проблемы начальных значений:
% equation rhs
function [f] = State(t,q,L);
f = L*q
endfunction;
T = 200.0
q0 = [ 0.1; 0.2 ];
[qT] = ode(q0, 0, T, list(State, L));
Как видите, логика кода проста, но я не могу понять, как преобразовать последнюю строкук коду Matlab, потому что я не использовал SCILAB раньше.В противном случае преобразование его в код Matlab должно быть простым, например, изменив endfunction
в определении функции на end
и добавив или удалив некоторое ;
.Не могли бы вы помочь мне преобразовать язык SCILAB, особенно list(...)
, в язык Matlab?Заранее спасибо.