Re: NMinimize Bug in Mathematica 7.0?
- To: mathgroup at smc.vnet.net
- Subject: [mg96320] Re: NMinimize Bug in Mathematica 7.0?
- From: Jean-Marc Gulliet <jeanmarc.gulliet at gmail.com>
- Date: Wed, 11 Feb 2009 05:23:31 -0500 (EST)
- Organization: The Open University, Milton Keynes, UK
- References: <200902091033.FAA12301@smc.vnet.net> <gmrltf$9h0$1@smc.vnet.net>
In article <gmrltf$9h0$1 at smc.vnet.net>,
"Jaccard Florian" <Florian.Jaccard at he-arc.ch> wrote:
> I would avoid to use numerical methods without looking at a graphic.
> As you can see:
> Plot[x - 3*x^2 + x^4,{x, -2, 2}]
>
> Both answers are correct, they give approximations of local minimums.
> You can ask for a minimal and maximal bound for the local minimums like =
> this :
>
> In[8]:= NMinimize[x - 3*x^2 + x^4,{x, -2, -1}]
>
> Out[8]= {-3.513905038934789,
> {x -> -1.3008395739047898}}
>
> In[9]:= NMinimize[x - 3*x^2 + x^4,{x, 0.5, 1.5}]
>
> Out[9]= {-1.0702301817761541,
> {x -> 1.1309011226299863}}
>
> If you want the absolute minimum, use a non-numeric method :
>
> In[11]:= N[Minimize[x - 3*x^2 + x^4, x]]
>
> Out[11]= {-3.5139050389347894,{x -> -1.3008395659415772}}
NMinimize[] is designed to find *global* minima. The text cell just
above the input cell of the first example on the online documentation
for NMinimize[] reads:
"Find the global minimum of an unconstrained problem"
Then the following example is given,
In[1]:= NMinimize[x^4 - 3 x^2 + x, x]
Out[1]= {-1.07023, {x -> 1.1309}}
And the result is clearly wrong since it is not a *global* minimum
contrary to what is claimed in this very same documentation.
One should never forget that the Mathematica documentation centre is a
"live" environment, that is, one can re-evaluate the cells (great
feature), which also means that when the documentation for a given
release is generated anew, therefore the new bugs in Mathematica
functions are also carried on in the documentation itself.
Regards,
--Jean-Marc
- References:
- NMinimize Bug in Mathematica 7.0?
- From: appris@att.net
- NMinimize Bug in Mathematica 7.0?