function y=harm(n)
    suma = 0;
    for k = 1:n do
        suma = suma + 1/k;
    end
    y = suma
endfunction

for k = 1:100 do
    x(k)= 100*k
    eul(k)=harm(100*k) - log(100*k)
end
blad = abs(eul-0.5772156649)
X = log(x)
Y = log(blad)
[a, b] = mnkw(X, Y)
disp(-a, "Rząd zbieżnosci metody obliczania przyblizenia stalej Eulera")