function s=harm(n)
k=1:n;
k=1./k;
s=sum(k);
