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;