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

m=5;        % antall punkter

% w og p koordinatene til de gitte punktene:

w=[44; 61; 81; 113; 131];

p=[91; 98; 103; 110; 112]; 


figure(10)
plot(w, p,'or')
hold on;            %plotter de gitte punktene

X = [ones(m,1) log(w)]    % X er designmatrisen

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

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

t= linspace(40,140); 

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

plot(t,f,'b') 

 p_100= beta(1) + beta(2)*log(100)
 % approksimerer p n錼 w = 100



