func6p3try.m

来自「数值分析最常用的四十种算法」· M 代码 · 共 6 行

function dy = func6p3try(t,y)
dy = zeros(3,1);    % a column vector
dy(1) = y(2) * y(3);
dy(2) = -y(1) * y(3);
dy(3) = -0.51 * y(1) * y(2);