function problem19(r)
% try, e.g., problem19(2e-5)
x1=1-r;
x2=1+r;
n=500;
x=linspace(x1,x2,n);
y=x.^3-3*x.^2+3*x-1;
plot(x,y,'.')