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

clear;
x=1.;
s=1.;
while x<=5000.
    x=x+1.;
    s=s+1./x;
end
disp([x,s]);

% Euler Constant

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