The Math Forum

Search All of the Math Forum:

Views expressed in these public forums are not endorsed by NCTM or The Math Forum.

Math Forum » Discussions » Inactive » comp.soft-sys.math.mathematica

Notice: We are no longer accepting new posts, but the forums will continue to be readable.

Topic: Low precision exponentiation
Replies: 5   Last Post: Feb 18, 2013 6:00 AM

Advanced Search

Back to Topic List Back to Topic List Jump to Tree View Jump to Tree View   Messages: [ Previous | Next ]
Albert Retey

Posts: 688
Registered: 7/15/08
Re: Low precision exponentiation
Posted: Feb 18, 2013 5:59 AM
  Click to see the message monospaced in plain text Plain Text   Click to reply to this topic Reply


> I am trying to evaluate 2.5^125 to high precision.
> R gives 5.527147875260445183346e+49 as the answer but Mathematica
> with N[2.5^125,30] gives 5.52715*10^49 and says that is to machine
> precision.

There are two problems here:

1) with the standard settings, Mathematica will print only 6 digits of
machine precision numbers. You can change that in the preferences tab or
with e.g.:

NumberForm[2.5^125, {16, 16}]

2) the input you are giving is interpreted to be machine precision in
the first place. You can get an exact result like this:


and a numeric result to the desired precision like this:


alternatively you can define the precision with this syntax:


> I am inexperienced at Mathematica. Am I doing something silly?

yes and no, Mathematica handles arbitrary precision numbers in a
somewhat unusual way. You might want to read the tutorials that are
linked to the documentation of N to learn about some details. On the
other hand, I think usually you can get away with the rule that one
should try to stay with exact results as long as possible when
requesting more than machine precision...



Point your RSS reader here for a feed of the latest messages in this topic.

[Privacy Policy] [Terms of Use]

© The Math Forum at NCTM 1994-2018. All Rights Reserved.