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

g=[10 0];
f=[1 -1];
R=tf(g,f,1);