Re: Solve or Reduce on a monstrosity of an expresssion (and a prize!)
- To: mathgroup at smc.vnet.net
- Subject: [mg57391] Re: Solve or Reduce on a monstrosity of an expresssion (and a prize!)
- From: Daniel Lichtblau <danl at wolfram.com>
- Date: Thu, 26 May 2005 04:31:44 -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>
- Sender: owner-wri-mathgroup at wolfram.com
Daniel Reeves wrote:
> yikes. ok, let's assume that all works. So it remains to show that both
> of the following expressions of n are zero for all n >= 2. Strangely,
> FullSimplify can do this for all specific n up to 113 but 114 seems to
> stump it.
>
> here are the expressions:
> [...]
I can point out some aspects to a probabilistic argument. Making it into
a full fledged proof would require figuring out some bounds that offhand
I don't know how to do, then verifying for appropriately many cases.
First you need to show that there are 2 roots for each x, at least when
n>=2. For this it suffices to show it holds for generic n. Since we
cannot seem to do that computationally, instead do it for "random" n.
Calling your expression fp, do, say,
fp7 = Simplify[fp /. n->7]
gb7 = GroebnerBasis[fp7,x]
Now check for a quadratic polynomial (in x) in gb7. If you find one,
that means there are at most 2 solutions when n=7. This can of course be
repeated for other values of n. The point is that, for integers larger
than 2, this computation appears to be feasible. Offhand I think this
idea (of verifying at sufficiently many values of n) suffices, though
there may be a lurking branch cut issue. Note that showing there are two
solutions for n=7 is not quite sufficient in and of itself because 7
might not be generic. In this setting, nongeneric points would lie on a
variety (roots of a polynomial), hence be a set of measure zero.
Verification that x->(n-1)/n and x->-(n+1)^2/(n*(n-1)) are the two
solutions can be done by making the substitutions into fp and verifying,
for sufficiently many n>=2, that these vanish. The reason is that
univariate algebraic functions vanishing at sufficiently many integer
points vanish everywhere that branches don't change (e.g. x-Sqrt[x^2]
vanishes on positive integers and indeed on all positive reals). I think
the branch cut issue goes away when you have the proviso that n>=2,
because for the substitutions above all radicands are positive in that
region. But maybe I'm missing something and there is more to it than this.
I should mention that if you did not know those were the appropriate
substitutions, they might be deduced by finding particular {x,n} pairs
and getting more values via homotopy continuation (basically, setting up
and solving an ODE via NDSolve). Values from this could then be used to
interpolate an implicit relation. This was the gist of some work I
described in a note last week:
http://forums.wolfram.com/mathgroup/archive/2005/May/msg00627.html
Offhand I do not know how well it would work for your particular
algebraic expression, but it strikes me as a plausible approach. (A
better one for this case might be to take an implicit plot of the zero
set of fp and guess the relation).
Daniel Lichtblau
Wolfram Research
- References:
- Solve or Reduce on a monstrosity of an expresssion (and a prize!)
- From: Daniel Reeves <dreeves@umich.edu>
- Solve or Reduce on a monstrosity of an expresssion (and a prize!)