Re: Eliminate
- To: mathgroup at smc.vnet.net
- Subject: [mg116178] Re: Eliminate
- From: Daniel Lichtblau <danl at wolfram.com>
- 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}, MonomialOrder->EliminationOrder, CoefficientDomain->RationalFunctions]];] 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 > > http://mail.google.com/mail/?shva==1#inbox/12deb2a54c93072f > > and I have checked that the equation is that can be downloaded there: > > http://bernard.gibert.pagesperso-orange.fr/curves/Resources/Q077.rtf > > 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