% Evaluating p = (x-2)^9 in x = 1.920, 1.921, ... , 2.080

p = [1 -18 144 -672 2016 -4032 5376 -4608 2304 -512];

x = linspace(1.92,2.08,161);
y = polyval(p,x);
z = (x-2).^9;

plot(x,y,'b',x,z,'r')