x=linspace(0.99998,1.00002,50);
y=x.^3-3*x.^2+3*x-1;
% plot(x,y,':')
plot(x,y)

