% Harjutustund 10, ylesanne 6
x=[1.2:0.1:2];
y=[9.08,10.43,11.9,13.48,15.19,17.03,19.01,21.13,23.39];


% võrdsed kaalud
plot(x,y,'o')
% kuupsplainiga
ys3=interp1(x,y,x,'spline');
hold on
plot(x,ys3)
% lähendame vähimruutude mõttes kuupfunktsiooniga

for i=1:4
  for j=1:4
      A(i,j)=sum(kaalud.*x.^(8-i-j));
  end
  B(i,1)=sum(kaalud.*y.*x.^(4-i));
end
A
B
C=A^(-1)*B
f=@(x) C(1)*x^3+C(2)*x^2+C(3)*x+C(4);
hold on
fplot(f,[0,4])

% leiame funktsiooni väärtused
% x=1.12
[interp1(x,y,1.12,'spline') f(1.12)]
% x=1.35
[interp1(x,y,1.35,'spline') f(1.35)]