🚀 go-pugleaf

*** REPLACE THIS LINE WITH YOUR MESSAGE ***

Before starting to argue as to which machine has the exact result, we should
first discuss the algorithm.
 
Calculating the harmonic by:
 
    for  ( i=1; i <= 10000; i++ )
        j += ( 1.0 / i );
 
is clumsy. When i reaches the big numbers, a lot of small significant digits
will be lost when we perform the addition.    For example, lets suppose that
the machine has 8 decimal digits of accuracy, and that the sum for
i = 1 to 9499 is  9.7325632.   Adding 1.0/9500 (1.0526316E-4) to our sum
will be done by first adjusting the small number to 0.0001053, and then
adding it to 9.7325632.  We have lost quite a few significant digits in
the process.
 
A better way is to write:
 
    for  ( i=10000; i > 0; i-- )
        j += ( 1.0 / i );
 
In this case, the Vax/780 comes with the answer 9.787604 (and not 9.787613)!
I would be curious to see what are the results on other systems.

			Jean-Luc Brouillet