% Sauer 2.3.2
% Denna kod skrevs på övningen 2017-09-20 i SF1514/18/19

clear

% (a) Linjär interpolation
% p(x) = c1 + c2*x

% Lös ut c1 och c2
A = [1 -10; 1 -5];
b = [1.95; 3.01];
c = A\b

% Evaluera polynomet
x = -8;
p = c(1) + c(2)*x

% Följande skrevs inte på själva övningen!
% (b) Kvadratisk interpolation
% p(x) = c1 + c2*x + c3*x^2

% Lös ut c1, c2 och c3
A = [1 -10 (-10)^2;
     1 -5 (-5)^2;
     1 0 0];
b = [1.95; 3.01; 4.58];
c = A\b

% Evaluera polynomet
x = -8;
p = c(1) + c(2)*x + c(3)*x^2

% I praktiken är det lättare att använda Matlabs inbyggda
% funktioner polyfit och polyval. Prova t.ex.
    coeff = polyfit([-10;-5;0], b, 2)
% Vad blir skillnaden på coeff och c?
