[Date Index]
[Thread Index]
[Author Index]
Re: Root finding needs higher accuracy
*To*: mathgroup at smc.vnet.net
*Subject*: [mg123210] Re: Root finding needs higher accuracy
*From*: DrMajorBob <btreat1 at austin.rr.com>
*Date*: Sat, 26 Nov 2011 05:09:41 -0500 (EST)
*Delivered-to*: l-mathgroup@mail-archive0.wolfram.com
*References*: <201111250957.EAA11210@smc.vnet.net>
*Reply-to*: drmajorbob at yahoo.com
The poster wanted 25 digit accuracy, and WorkingPrecision->20 yields only
6 correct digits:
f[x_] := Log[x] + x^2/(2 Exp[1]) - 2 x/Sqrt[Exp[1]] + 1
root = x /. FindRoot[f12[x], {x, 2}, WorkingPrecision -> 20]
Replace[RealDigits@N[root, 20] - RealDigits@N[Sqrt@E, 20] //
First, {x : Longest[0 ..], ___} :> Length@{x}]
1.6487220277297822898
6
Bobby
On Fri, 25 Nov 2011 03:57:35 -0600, Bill Rowe <readnews at sbcglobal.net>
wrote:
> On 11/24/11 at 6:57 AM, hwborchers at googlemail.com (HwB) wrote:
>
>> I would like to numerically find the root of the following function
>> with up to 20 digits.
>
>> f12[x_] := Log[x] + x^2 / (2 Exp[1]) - 2 x / Sqrt[Exp[1]] + 1
>
>> This problem is known to be difficult for solvers in double
>> precision arithmetics. I thought it should be easy with Mathematica,
>> but the following attempts were not successful.
>
>> SetPrecision[ x /. FindRoot[f12[x], {x, 1.0, 3.4}, Method ->
>> "Brent", AccuracyGoal -> Infinity, PrecisionGoal -> 20], 16]
>> #1.648732212532746
>> SetPrecision[ x /. FindRoot[f12[x], {x, 1.0, 3.4}, Method ->
>> "Secant", AccuracyGoal -> Infinity, PrecisionGoal -> 20], 16]
>> #1.648710202030051
>
>> The true root obviously is Sqrt[Exp[1]]//N = 1.648721270700128...
>
>> The symbolic solver explicitely says it cannot solve this
>> expression. What do I need to do to get a much more exact result out
>> of Mathematica?
>
> On my machine the simple naive approach seems to work fine. That is
>
> In[25]:= FindRoot[f12[x], {x, 2}, WorkingPrecision -> 20]
>
> Out[25]= {x->1.6487220277297822898}
>
> In[26]:= $Version
>
> Out[26]= 8.0 for Mac OS X x86 (64-bit) (October 5, 2011)
>
>
--
DrMajorBob at yahoo.com
Prev by Date:
**Re: How to integrate a function over a polygon**
Next by Date:
**Re: Using Equal with Real Numbers**
Previous by thread:
**Re: Root finding needs higher accuracy**
Next by thread:
**Re: Root finding needs higher accuracy**
| |