a = 0; n = 0; while (a<5) a = a + 1/(n+1); n = n + 1; end a = a - 1/n; n = n - 1; disp(a),disp(n)