x=[-2 -1 2]';
y=[-2 1 -1]';

A=[x.^2 x ones(size(x))];
A=vander(x);

c=A\y

plot(x,y,'o');


hold on; grid on;
X=linspace(-3,3,200);
Y=polyval(c,X);
plot(X,Y);