function s=harmonisk1(n)
s=sum((1:n).^-1);
