% EXS 4.3 av Erik Dalsryd 2017-09-28
clear all, close all, clc

t = [0  1  3  4]';
y = [1  8  11  20]';

plot (t, y, 'o')
hold on

disp ('Linjär')
A = [t.^0  t]   % Skapa matisen
c = A\y     % Beräkna c i y=c1 + c2*t
r = y - A*c   % Beräkna residualen
plot (t, A*c )

disp ('Vågrät')
A = ones (size(t))
c = A\y     % Beräkna c i y=c
r = y - A*c
plot (t, A*c, '--')

disp ('Genom origo')
A = t
c = A\y     % Beräkna c i y=c*t
r = y - A*c
plot (t, A*c, '-.')

legend ('Mätpunkter','Linjär','Vågrät','Genom origo')
xlabel ('tid'), ylabel ('y')