s=0; j=1; while (s<5) s=s+(1/j); j=j+1; end s1=s-(1/j); j=j-1; disp(s1) disp(j)