clear all
close all
clc
M=input('Podaj stopien wielomianu m' );
X=input('Podaj wektor X: ');
Y=input('Podaj wektor Y: ');

N=size(X,2);
N2=size(Y,2);
if( N ~= N2 )
    error ( ' Liczba el. X != liczbie el. Y ');
end
if( M > N-1 )
    error ( 'Stopien wielomianu jest wiekszy od liczby wezlow ');
end
M=M+1;
V=ones(M);
TMP=zeros(M);

for j = 1:N

    for k = 1:M  %Macierz KSI
        for i = 1:M
            V(k,i) = X(j)^(k+i-2);
        end
    end

    for k = 1:M  %Dodawanie do macierzy tymczasowej
        for i = 1:M
            TMP(k,i) = TMP(k,i) + V(k,i);
        end
    end

end

F=zeros(M,1);
Y=Y';

for j = 1:N;
    for k=1:M
        F(k,1)= F(k,1) + X(j)^(k-1)*Y(j);
    end
end

A=(TMP^(-1))*F;
w=0;
x=min(X):0.1:max(X);

for j = 1:M
   w=w+A(j)*x.^(j-1);
end

disp(A)
plot(x,w,'b-')
title('Aproksymacja')
hold on
plot(X,Y,'or')