>> num=[1 2 1]; >> den=[1 1 1 1]; >> sys=tf(num, den) Transfer function: s^2 + 2 s + 1 ----------------- s^3 + s^2 + s + 1 >> bode(sys) >> step(sys) >> bode