% Eksempel -- Lineaer modell som i 6.6.8
clear all;

m=8;        % antall punkter

%x og y koordinatene til de gitte punktene:

x=[4; 6; 8; 10; 12; 14; 16; 18];

y=[1.58; 2.08; 2.5; 2.8; 3.1; 3.4; 3.8; 4.32]; 


figure(6)
plot(x, y,'or')
hold on;            %plotter de gitte punktene

X = [x x.^2 x.^3]    % X er designmatrisen

M = rref([X y])      % sjekker om systemet er konsistent og 
                     % om X har lin. uavh. kolonner

beta=X'*X \ X'*y     % gir minste kvadraters losning

t= linspace(4,18); 

f= beta(1)*t + beta(2)*t.^2 + beta(3)*t.^3;  %gir minste kvadraters tilpasning

plot(t,f,'b')   



