MathGroup Archive 2010

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

Search the Archive

might be useful

  • To: mathgroup at smc.vnet.net
  • Subject: [mg110293] might be useful
  • From: Peter Pein <petsie at dordos.net>
  • Date: Sat, 12 Jun 2010 05:31:33 -0400 (EDT)

As apology for my silly question, I'd like to provide a code snippet I
wrote a while ago.
Sometimes the problem arises to solve an equation/a set of equations
not for a variable, but for an expression. Examples:

In[1]:= Isolate[eqn:(_Equal|{_Equal ..}),expr_,elim_:{}]:=
Block[{tmp},
Solve[GroebnerBasis[Subtract@@@Flatten[{eqn,tmp==expr}],tmp,
Flatten[{First[Variables[expr]],elim}]]==0,tmp][[1]]/.{Rule->Equal,tmp->expr}]
In[2]:= Isolate[{x^2+3y^2==1,y-x==1/2},2x+y]
Out[2]= {2 x+y==1/8 (-5-3 Sqrt[13])}
In[3]:= Isolate[x^2+3y^2==1,2x+y]
Out[3]= {2 x+y==y-2 Sqrt[1-3 y^2]}
In[4]:= Isolate[{x^2+y^2+z^2==1,r Cos[t]+r Sin[t]==1},x+y,{r,t}]
Out[4]= {x+y==y-Sqrt[1-y^2-z^2]}
or
Isolate[{x^2+y^2+z^2==1,3z-(x+y)^2/2==0},x+y,z]
eliminates z but the answer is a bit lengthy.


If one or some of you want to improve it (and do this successfully),
please let me know.

Cheers,
Peter


  • Prev by Date: Find value of unknown const that causes integral to equal some value
  • Next by Date: Re: using rules for square roots
  • Previous by thread: Re: Find value of unknown const that causes integral to equal some value
  • Next by thread: Re: Handheld mathematica