for k=1:50 do
    n=100*k
    x(k)=n
    y(k)=harm(k)
    z(k)=log(k)
    eul(k)=harm(k)-log(k)
end

blad=abs(eul-0.5772156649)

logEul=log(eul)
logBlad=log(blad)

[a1, a2] = mnkw(logEul, logBlad)

disp(-a1)