% Harmonic Series. Computer Problem (7) Section 2.2 Cheney-Kincaid

% Euler Constant

s=0.;
k=1;
for N=1:5000
    s=s+1./N;
    g=s-log(N);
end
disp([N,g]);