while loop conditional test problem

Peter Goodman schrieb: Hi, I am having a problem getting a while loop conditional test to behave as I expect. See the attached example. Any ideas what I'm doing wrong? Thanks, Pete.

You should expect function_error to always return 0.1, because you have used <= instead of <. So when i==x an additional 1/10 is added to i and you function should always return x+0.1.

And even this is not really true - only for values of x which are multiples of 0.1.

Concerning the rest and main part of your question the answer is: roundoff errors by numeric evaluations. You get the correct results by evaluating symbolical. When evaluating numerically set the number of decimal places to the highest possible numer (17) and have a look. Best guess is that numerical evaluation use the "normal" way to deal with numers and that is converting them to binary. 0.1 in binary has an infinite number of decimal point an has to be truncated.

By adding 0.1 ten times mathcad got a value slightly under 1 and therefore added an additional 0.1.

by adding 0.1 20 times it got a value slightly over 2 and the while loop stopped (correctly).

You can try working with integers as much as possible

But even then you get slight rounding errors if x is not integer.

As long as evaluations are done numerically you will run in such situations. You might try not to ask if a==b but if abs(a-b)<eps where eps is a small value like 10^-6. Or in your case you might ask for i<=x+eps.

WE