[T,Y] = ode15s('minus_1000_y',[0,1], [1]);
plot(T,Y)
hold
plot(T,zeros(size(T,1),1),'LineStyle','none','Marker','+',...
'MarkerEdgeColor','red','Markersize',4);
size(T)
