sub6p3.m

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

function dy = sub6p3(t,y)
dy = zeros(3,1);    % a column vector
dy(1) = y(2);
dy(2) = y(3);
dy(3) = -1 * y(1) -2* y(2)-2*y(3)+80;