function sum = hsum(n)
% Sum of the first n terms of the harmonic series

% Written by Paul Chew for CS100M, Feb 2006

sum = 0;
for i = 1:n
	sum = sum + 1/i;
end