w=tf([1 0 2],[1 0 4 1])
w2=feedback(w,1)
step(w)