clf;
plot (t, y, t, u*0.1);
hold on;
plot([0, 1, 1, 5], [1.1, 1.1,1.05, 1.05], 'r--');
plot([1, 1, 5], [0, 0.95, 0.95], 'r--');
