% LKNM
x = [1 4 16 64]';
y = [4 38 60 22]';
% f(x) = a * x + b * sqrt(x) + c * log2(x)
A(:,1) = x;
A(:,2) = sqrt(x);
A(:,3) = log2(x);
A
% papirra  2 sor:
A'*A
A'*y
mo = (A'*A)\(A'*y)
% papirra:
a = mo(1)
b = mo(2)
c = mo(3)
% papirra:
% f(x) = -2 * x + 6 * sqrt(x) + 17 * log2(x)
function y = f(x)
  y = -2 * x + 6 * sqrt(x) + 17 * log2(x);
endfunction
t = 0:0.1:65;
plot(t,f(t),x,y,'*')