>>syms x n;
>>sol=limit(symsum(1/x,x,1,n)-log(n),n,inf)
sol =
eulergamma
>> vpa(sol,20)
ans =
0.57721566490153286061