Search All of the Math Forum:

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

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

Topic: Problem finding maximum
Replies: 5   Last Post: Aug 15, 2012 5:28 AM

 Messages: [ Previous | Next ]
 Peter Pein Posts: 1,147 Registered: 5/4/05
Re: Problem finding maximum
Posted: Aug 15, 2012 5:28 AM

Am 06.08.2012 10:40, schrieb Cisco Lane:
> Hmm - My Mathematica (Mathematica 7 Home Edition, Ver 7.0.1.0, Mac OS X X86 (32-bit)) gives a different answer:
>
> f[x_, a_] = (a^3 - 6 x - a^2 (4 + x) + a (2 + 12 x - 4 x^2))/(8 a);
> aa1 = .7481;
> NMaximize[{Abs[f[x, aa1]], 0 <= x <= aa1}, x]
>
> {0.0274936, {x -> 0.403948}}
>
> Any idea why? I must use an automatically selected start point, because the values of aa1 vary widely.
>
> Using the rational aa1 and Maximize works:
>
> f[x_, a_] = (a^3 - 6 x - a^2 (4 + x) + a (2 + 12 x - 4 x^2))/(8 a);
> aa1 = 7481/10000;
> N[Maximize[{Abs[f[x, aa1]], 0 <= x <= aa1}, x]]
>
> {0.0540933, {x -> 0.}}
>
> but I am not sure why. If this method is reliable, I could use it, but it seems klugy.
>

Seemingly there is a bug in version 7:
NMaximize[{Abs[f[x, aa1]], 0 <= x <= aa1}, x,
Method -> {Automatic, "InitialPoints" -> List /@ Range[0, 3] aa1/3}]

{0.0274936, {x -> 0.403948}}

Even if the global maximum location is in the set of initial points, the
default method fails in this case.

This does not happen any more in version 8.
If an upgrade is no option, try Bob's tip using specific Methods.

You can learn a lot about Mathematica 7 if you try every example in the
tutorial [1] written for version 8, I guess ;)

Peter

[1]
http://www.wolfram.com/learningcenter/tutorialcollection/ConstrainedOptimization/

Date Subject Author
8/5/12 Bob Hanlon
8/5/12 Frank K
8/5/12 Dana DeLouis
8/6/12 Cisco Lane
8/15/12 Peter Pein