% Eksempel -- Lineaer modell med polynomer i P_2

m=4;        % antall punkter

x=[0; 1; 2; 3]; y=[1; 2; 1; 5]; %x og y koordinatene til de gitte punktene

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

X=[ones(m,1) x x.^2]     % X er designmatrisen

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

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

t=linspace(-1, 4); 

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

plot(t,p,'b')  % plotter m.kv. polynomet 



