%Problem 4  
sysa = tf([1],[1 -1 3])
step(sysa)

sysb = tf([1 -3],[1 5 4])
step(sysb)

sysc = tf([3 2],[1 8 7])
step(sysc)

sysd = tf([4.5],[1 6 25])
step(sysd)