sub6p3.m

数值分析最常用的四十种算法

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;