clear all, close all

x = [91 121 152 182 213 244]';
y = [13.2 15.8 18.0 18.4 16.6 14.1]';

plot(x,y,'-*'), xlabel('dag nr'), ylabel('dagens längd')

% det overbestammda systemet
A = [ones(size(x)) x x.^2]

% mkm
Am = A'*A;
ym = A'*y;
c = inv(Am)*ym; % inte bra!

% modellfunktionen
f = @(t) c(1) + c(2)*t + c(3)*t.^2;

% plotta
figure()
xx = linspace(x(1),x(end),100);
subplot(2,1,1)
plot(x,y,'-*',xx,f(xx),'r'); xlabel('dag nr'), ylabel('dagens längd')
subplot(2,1,2)
plot(x,y-f(x),'-*'); title('residualen');