% Interpolacja
clear all;
X = input('wektor a');
Y = input('wektor b');
for i = 1 : length(X)
	for j = 1 : length(Y)
		V(i,j) = X(i)^(j-1);
	end
end
A = V^(-1) * Y';

IX = min(X) : 0.1 : max(X);
IY = zeros(1, length(IX));
for j = 1 : length(X)
    IY = IY + A(j) * IX.^(j - 1);
end

plot(X,Y,'or');
hold on;
plot(IX,IY);
title('W(x) = a_0 + a_1 * x + _2 * x^2 + ... + a_n * x^(n-1)')
hold off;