Floating-point accuracy

[Ed Hill]

On Wed, 2006-08-09 at 12:21 -0500, Joel Parker wrote:
> I'm having some trouble with floating-point accuracy in g77 and gfortran,
> and I'm starting to think I'm missing something obvious.

Please try:

      N1 = 2453957.09236D0
      N2 = 2415020.0D0
      N3 = 36525.0D0

and see if that helps.  :-)

Ed

-- 
Edward H. Hill III, PhD
office:  MIT Dept. of EAPS 54-1424;  77 Massachusetts Ave.
             Cambridge, MA 02139-4307
emails:  eh3@mit.edu                ed@eh3.com
URLs:    http://web.mit.edu/eh3/    http://eh3.com/
phone:   617-253-0098
fax:     617-253-4464