From: saxena@ramanujan.cs.albany.edu (Tushar Saxena) Newsgroups: sci.math.symbolic Subject: Re: Q: Mma (or other) Eliminate Date: 24 Jan 1997 16:33:32 GMT The following variable elimination problem is easily solved by using multivariate resultants. One of the method, which we have found to be quite efficient, is based on the Dixon formulation. See the following papers (available from http://www.cs.albany.edu/~saxena) to read about it: Kapur D., Saxena T. and Yang L., Algebraic and Geometric Reasoning using Dixon Resultants, {\em Proc. ACM International Symposium on Symbolic and Algebraic Computation}, Oxford, England, July 1994.\\ Kapur D and Saxena T., Comparison of Various Multivariate Resultant Formulations, {\em Proc. ACM International Symposium on Symbolic and Algebraic Computation}, Montreal, July 1995.\\ Kapur D and Saxena T., Sparsity Considerations in the Dixon Resultant Formulation, {\em Proc. ACM Symposium on Theory of Computing}, Philadelphia, May 1996.\\ This method is quite easy to implement in your favourite computer algebra system, and even a naive implementation should work surprisingly well! You said: >Mathematica gives no answer after 24 hours in a Pentium or Sun: I used my MAPLE implementation of this method to solve your problem in less than 50 seconds. And that was after substituting a symbolic parameter for 3^(1/2). The result was a degree 8 polynomial in $x,y$ and had 81 terms. To save bandwidth, I am sending you the polynomial by email. Tushar Saxena Research Consultant General Electric, Corporate Research & Development -------------------------------------------------- In article <32E372D7.158D@uvigo.es>, Francisco Botana wrote: >Trying to obtain in Mathematica 2.2.3 the envelope of a family of lines >I've got some problems. > >If the lines are ellipses of expression > >with focus in the ellipse > >e1=-500 - 80*a + 4*a^2 + 9*b^2 > >e22=-80 + 21*a^2 + 32*b + 31*b^2 - 32*3^(1/2)*x + > 10*3^(1/2)*b*x + 21*x^2 - 32*y - 62*b*y - > 10*3^(1/2)*x*y + 31*y^2 + > a*(32*3^(1/2) - 10*3^(1/2)*b - 42*x + 10*3^(1/2)*y) > >Mathematica gives no answer after 24 hours in a Pentium or Sun: > >Eliminate[{e1==0,e22==0, > Dt[e1,a,Constants->{x,y,b}]* > Dt[e22,b,Constants->{x,y,a}]== > Dt[e22,a,Constants->{x,y,b}]* > Dt[e1,b,Constants->{x,y,a}]}, > {a,b}] > >... > >My questions are: > >How can I speed up the answer (if Mma can obtain it)? >Is there any other symbolic package capable of coping with this problem? >How? > >Thanks in advance. >-- > Francisco Botana > Depto. de Matematica Aplicada EUIT Forestal-Universidad de Vigo > Campus de A Xunqueira, s/n 36005 PONTEVEDRA (Spain) > tel +34-86-801931 fax +34-86-801907