b=conv([1 -4],[1 0.5]);
a=conv(conv([1 -2],[1 -1]),[1 -0.5]);
P=tf(b,a);