n=input('enter how many terms you want for your sum of sequence')
sumo=0
for ii = [1:n]
    sumo=sumo+(1/ii);
    if mod(ii,1000)==0
        fprintf('%f\n',sumo)
    end
end
sumo