Date: Jan 19, 2013 4:11 PM
Author: S.K.Mody
Subject: Re: How to avoid wrong answers in simple arithmetic expressions.

dpb <> wrote in message <kdeaav$3e3$>...
> On 1/19/2013 7:53 AM, dpb wrote:
> ...

> >
> > Scaling...[and other stuff elided for time being]
> >

> What happens if you were to multiply v by 100 and do the same w/ the
> input variable internal to the function before you start your loop?
> --

Thanks for the suggestion.

Yes, I did try this, and it gives the correct result although I'm not sure I understand why. My first thought was that operations on numbers which differed only in the base 10 exponent should give results of equal accuracy.

However, looking up the IEEE 754 floating point floating point representation. I see that the exponent is a power-of-2 exponent rather than a power-of-10. So, for example, 0.09 and 9 have very different mantissa's. Nevertheless both 0.09 and 9 do not have exact IEEE 754 representations. More generally, most integers do not have an exact representation in this format.

So, does Matlab treat small integers in a special way?
Otherwise why would multiplying the numbers by 100 in the original problem guarantee a correct solution?

Best Regards.