A=[1 1;1 3;2 -1;3 1];
b=[3 7 0.2 5]';

%zadatak 4.21

% mnozenjem sa A' obe strane od Ax=b, dobijamo normalni sistem A'*A*x=A'*b
% dimenzije nxn, za cije resavanje dalje koristimo operator \

coeffs = ((A'*A)\(A'*b))