n=100;
start=-3;
bixg=5;
x=start+(bixg/n)*[1:n];
y=x.*x.*x+1.5*x.*x-3*x+.5;
plot([x(1) x(n)],[0 0],x,y)