num = [2 4]; den = [1 5 6 4]; sys = tf(num,den); step(sys); impulse(sys);