%uloha 125

load skvrny

%prevod na vektor (po radcich)
s=data';    %z radku sloupce
s=s(:); %sloupce pod sebe do vektoru
s=s';   %radkovy vektor
s=s(1:end-2);   %vynechani poslednich dvou hodnot

%aproximace - MNÈ
d=length(s);
X=[sin(s') ones(d,1)];  %sloupcove
b=X\s';
s_aprox = X*b;  %aproximovane hodnoty

%graf
x=1:d;
plot(x,s,'kx',x,s_aprox,'k')

%extrapolace nasledujici hodnoty
s_11_2010=b(1)*sin(d+1) + b(2);
disp(sprintf('Listopad 2010: %f', s_11_2010))