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

% Euler Constant

format short;
%format short e;
%format long;
%format long e;
s=0.;
k=1;
for N=1:5000
    s=s+1./N;
    if mod(N,100)==0
        g=s-log(N + 1./2.);
        disp([N,g]);
    end
end