% Harmonic Series. Computer Problem (18) Section 1.3 Burden

% Euler Constant

clear;
format short e;
s=0;
N=1.;
g=0;
while N<=500000 & abs(0.5772156649-g)>=10^(-6)
    s=s+1./N;
    g=s-log(N);
    N=N+1.;
end
disp([N-1,g]);