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