%interpolace polynomem

%zadana (zmerena) x a y
% x=[0.1 1 1.1 1.2 1.3 2    3    4    5   ]';
% y=[1   1 0.9 0.2 0.1 0.01 0.01 0.01 0.01]';
%v techto bodech chceme spocitat interpolaci
x = [0.5 1 2.5 3.8 4.5]';
y = [4.305 3.36 -7.875 27.72 125.55]';

xToPrint=0.1:0.01:5;


%sestrojime vandermondovu matici, matlab si ji nadefinoval jinak, tak ji
%preklopime
matice=fliplr(vander(x));

%resime soustavu matice * koeficienty = y
koeficienty=matice\y;

%spoctememe si hodnoty v xToPrint, opet treba preklopit
yToPrint=polyval(flipud(koeficienty),xToPrint);

%zobrazime originalni data jako cervene krizky a navic interpolaci jako caru
plot(x,y,'rx',xToPrint,yToPrint,'g-');
title('Polynom stupne 8');