function test
[t,y]=ode45(@ps3_1_3_hill_fun,[0,25],0);
plot(t,y)
function dydt=ps3_1_3_hill_fun(t,y)
H=0;
if 10<t<15
H=1;
end
dydt=1900*H+100-2*y;
--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 60.250.64.98
※ 文章網址: https://www.ptt.cc/bbs/MATLAB/M.1483808583.A.C8F.html