clear all
clc

N = 5;


X = [-2 -1 1 3 4];
Y = [-90; -80; 72; 80; -180];

Wand = zeros(N,N);

for i = 1:N
   for z = 0:(N-1)
      Wand(i,z+1) = X(i)^z;
   end
end
disp(Wand)
A = (Wand'*Wand)^(-1)*Wand'*Y

syms x;
for i=1:N
y(i) = A(i)*x.^(i-1)

end
y=sum(y)
x=-3:0.1:5

y=eval(y)

plot(X,Y,'o',x,y)