plot(tout,y), grid on, xlabel('time'), ylabel('output')