subplot(311), plot(time, Y_delay), grid on, ylabel('output delay')
subplot(312), plot(time, Y_lagging1), grid on, ylabel('output lagging1')
subplot(313), plot(time, Y_lagging2), grid on, ylabel('output lagging2'), xlabel('time [sec]')