% Accmulation of errors
h=0.000001;
x = 1-20*h:h:1+20*h; % x is an array of length 41
plot(x,(x-1).^3,'.b',x,x.^3-3*x.^2+3*x-1,'.r')