Mathematica 9 is now available
Services & Resources / Wolfram Forums
-----
 /
MathGroup Archive
2005
*January
*February
*March
*April
*May
*June
*July
*August
*September
*October
*November
*December
*Archive Index
*Ask about this page
*Print this page
*Give us feedback
*Sign up for the Wolfram Insider

MathGroup Archive 2005

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: Re: Solve or Reduce on a monstrosity of an expresssion (and a prize!)

  • To: mathgroup at smc.vnet.net
  • Subject: [mg57493] Re: [mg57487] Re: Solve or Reduce on a monstrosity of an expresssion (and a prize!)
  • From: Andrzej Kozlowski <akoz at mimuw.edu.pl>
  • Date: Sun, 29 May 2005 21:00:10 -0400 (EDT)
  • References: <200505230620.CAA04045@smc.vnet.net> <EEA0FBBD-9C31-419A-8D0B-7C73AE4DA32E@akikoz.net> <Pine.LNX.4.58.0505231817420.9452@boston.eecs.umich.edu> <458D701E-37FA-425F-89C4-52A5628E22CF@akikoz.net> <Pine.LNX.4.58.0505241736010.24188@boston.eecs.umich.edu> <1D91709F-1E4D-4B98-95F5-695F7BD65577@akikoz.net> <Pine.LNX.4.58.0505242049410.3527@boston.eecs.umich.edu> <F016112F-BF6B-4D4E-8A3C-1F38AB3DDC5F@akikoz.net> <Pine.LNX.4.58.0505251303210.3527@boston.eecs.umich.edu> <D20ACE19-9F36-432F-9C52-9D059659BEAA@akikoz.net> <d76ojk$7ru$1@smc.vnet.net> <200505290503.BAA25023@smc.vnet.net>
  • Sender: owner-wri-mathgroup at wolfram.com

Very impressive indeed. In fact, one of the most impressive things I  
have ever seen on the MathGroup.
This really deserves a prize.

Andrze


On 29 May 2005, at 14:03, Maxim wrote:

>
> On Fri, 27 May 2005 09:17:40 +0000 (UTC), Daniel Reeves
> <dreeves at umich.edu> wrote:
>
>
>> Well, it's for a theorem in my phd thesis about the efficacy of a  
>> game
>> theory approximation technique in a class of simple auction games.
>>
>> To open this up to non-Mathematica folks, here's the function (the
>> derivative of the original function I posted) in a platform-neutral
>> format:
>>   http://ai.eecs.umich.edu/people/dreeves/f.txt
>>
>> Again, the goal is to show that that expression is negative in
>> (0,(n-1)/n)
>> for all integers n>=2.  One way to show that is to show that the
>> expression has at most 2 roots and show that it is zero at x =  
>> (n-1)/n
>> and
>> at x = -(n+1)^2/(n*(n-1)).  (It is.)
>>
>>
>
> Sometimes we can make it easier for GroebnerBasis if we combine common
> subexpressions; also, with the proper choice of replacements, we  
> can make
> the coefficients integer:
>
> In[1]:=
> f[x_, n_] := 9/2/c[x, n]^2*(n + 1)b[x, n]^2 (x - d[x, n])(x - x*d 
> [x, n]
> + d[x, n] + d[x, n]^2 + n (d[x, n] - 1) (x + d[x, n]))
>
> a[x_, n_] := 9*(n + 1)^2 + Sqrt[3(n + 1)^3 (x^2 (n - 1) + 27(n + 1))];
>
> b[x_, n_] := (a[x, n](n - 1) x^2)^(1/3);
>
> c[x_, n_] := -3^(2/3) x^2 (n^2 - 1) + 3^(1/3)(x^2(n^2 - 1) (9 + 9n
> + Sqrt[3(n + 1) (x^2(n - 1) + 27(n + 1))]))^(2/3);
>
> d[x_, n_] := c[x, n]/(3 b[x, n] (n + 1));
>
> In[6]:=
> D[f[x, n], x] /.
>          27*(1 + n) + (-1 + n)*x^2 -> 3*a^2/(1 + n) /.
>        (-1 + n^2)*x^2*(9 + 9*n + 3*Sqrt[a^2]) |
>            (-1 + n)*x^2*(9*(1 + n)^2 + 3*Sqrt[(1 + n)^2*a^2]) ->
>          9*b^3 //
>      Refine[#, {n > 0, a > 0, b > 0}]& //
>    Together // Numerator // FullSimplify
>
> Out[6]=
> 4374*a*b^10 - 486*a*b^8*(-1 + n^2)*(a + 12*(1 - b + n))*x -  
> 4374*a*b^8*(-1
> + n^2)*x^2 - 81*b^5*(-1 + n^2)^2*(b^3 + 2*a^2*(b - 6*(1 + n))
> + 12*a*(3*b^2 + 2*b*(1 + n) - 3*(1 + n)^2))*x^3 + 1458*a*b^6*(-1
> + n^2)^2*x^4 + 27*b^3*(-1 + n^2)^3*(12*a^2*(1 + n) + b^2*(6 - b + 6*n)
> + 18*a*(b^2 + 2*(1 + n)^2))*x^5 - 162*a*b^4*(-1 + n^2)^3*x^6
> + 18*b^2*(a*(3 + a - 3*b) + 3*b + 3*(a + b)*n)*(-1 + n^2)^4*x^7 + 3* 
> (-1
> + n^2)^5*(b^2 + 2*a*(3 + a + 3*n))*x^9 + (-1 + n^2)^6*x^11
>
> In[7]:=
> GroebnerBasis[
>      {%,
>       (1 + n)*(27*(1 + n) + (-1 + n)*x^2) - 3*a^2,
>       (-1 + n^2)*x^2*(9 + 9*n + 3*a) - 9*b^3},
>      {n, x}, {a, b},
>      MonomialOrder -> EliminationOrder] //
>    Factor // Timing
>
> Out[7]=
> {114.312*Second, {(-1 + n)^7*(1 + n)^8*x^13*(1 - n + n*x)*(1 + 2*n  
> + n^2 -
> n*x + n^2*x)*(27 + 27*n - x^2 + n*x^2)}}
>
> In[8]:=
> Reduce[%[[2]] == 0 && n > 1, x, Reals]
>
> Out[8]=
> n > 1 && (x == (-1 - 2*n - n^2)/(-n + n^2) || x == 0 || x == (-1 +  
> n)/n)
>
> These are all the possible real roots (not necessarily the set of  
> roots of
> D[f[x, n], x] exactly, because of zeros in the denominator; also there
> could be parasite solutions due to the presence of radicals).
>
> Maxim Rytin
> m.r at inbox.ru
>
>


  • Prev by Date: Re: Re: FindRoot: Failed to converge to the
  • Next by Date: What is causing this error message?
  • Previous by thread: Re: Solve or Reduce on a monstrosity of an expresssion (and a prize!)
  • Next by thread: plot hyperbola