Si=inline('a*x.^3 + b*x.^2 + c*x + d','x','a','b','c','d');

a=-1;
b=-3;
c=-1;
d=2;
e=1;
f=-3;
g=-1;
h=2;

x1=linspace(-1,0,100);
x2=linspace(0,1,100);

plot(x1,Si(x1,a,b,c,d),x2,Si(x2,e,f,g,h),'LineWidth',2);

x=[-1 0 1];
y=[1 2 -1];
hold on;
plot(x,y,'r.','MarkerSize',20);