MathGroup Archive 2011

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

Search the Archive

Re: Eliminate

  • To: mathgroup at
  • Subject: [mg116178] Re: Eliminate
  • From: Daniel Lichtblau <danl at>
  • Date: Fri, 4 Feb 2011 01:39:19 -0500 (EST)

Francisco Javier Garc=EDa Capit=E1n wrote:
> Hello, as I reposted on the forum, what I really wanted is to eliminate
> v and w

(I knew that.)

> (in trying to do this, I put v=1-t, w=t, this is why I did the
> mistake).

You can do this as follows. Subtract one side from the other to get
polynomials. Then extract a Groebner basis with an elimination ordering.

polys == {(v + w)*(b^2*v^4 - a^2*v^3*w + 3*b^2*v^3*w +
     c^2*v^3*w - 3*a^2*v^2*w^2 +
     3*b^2*v^2*w^2 + 3*c^2*v^2*w^2 - a^2*v*w^3 + b^2*v*w^3 +
     3*c^2*v*w^3 + c^2*w^4) - x,
   -(w*(-(b^2*v^4) - 2*b^2*v^3*w - a^2*v*w^3 + b^2*v*w^3 +
      c^2*v*w^3 + c^2*w^4)) - y,
  -(v*(b^2*v^4 - a^2*v^3*w + b^2*v^3*w + c^2*v^3*w -
     2*c^2*v*w^3 - c^2*w^4)) - z};

Either of these works. The first is faster, in your example. Isn't
always that way, though.

Timing[elim2 == First[GroebnerBasis[polys, {x,y,z}, {v,w},

Timing[elim == First[GroebnerBasis[polys /. {a->4,b->-7,c->11},
   {x,y,z}, {v,w}, MonomialOrder->EliminationOrder]];]

> By the way, I contacted with Bernard Gibert and he told me that the
> equation that I am looking for is the quintic Q077
> and I have checked that the equation is that can be downloaded there:
> I am still interested in how Mathematica can arrive from parametric
> equations in terms of v,w  to a implicit equation in x,y,z (a,b,c here
> are constants)

As above. If you give explicit integer or rational values for those
constants then it becomes hugely faster, I might add.

Daniel Lichtblau
Wolfram Research

  • Prev by Date: Re: DesignerUnits 2011-01-08 for Mathematica 8, 7, 6
  • Next by Date: Re: Mixing optional parameters and options
  • Previous by thread: Re: Eliminate
  • Next by thread: Re: Eliminate