p = [1 -2];
for i = 2 : 9
  p = conv(p,[1,-2]);
endfor
p
x = 1.92 + (0:8000)/8000*0.16;

y = p(1);
for i = 2 : length(p)
  y = x.*y + p(i);
endfor

z = (x-2).^9;

plot(x,y,'b')
hold on
plot(x,z,'y','LineWidth',3)
hold off