sys = tf([1 2],[1 2 7]);[u,t] = gensig('square',10,50);[ys,ts] = lsim(sys,u,t);hold onplot(t,u,'-.')plot(ts,ys,'r')