x=linspace(0,2);
plot(x, -1+x.*(3+x.*(-3+x)))

xx=linspace(.99998,1.00002);
plot(xx, -1+xx.*(3+xx.*(-3+xx)),'.')