MathGroup Archive 2006

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

Search the Archive

Re: the ellipse and the circle

Jack Kennedy wrote:
> Greetings Math People,
> I'm trying to solve the following problem.
> Please see the figure at
> A circle of radius  r  sits on an ellipse with semi-major axis  a  and
> semi-minor axis  b  as shown in the figure. The center of the circle
> is  theta  radians around the ellipse. The radius  r  is small
> compared to  a  and  b.  In terms of  theta, a, b, and r, what are the
> coordinates of the two points of intersection (shown in red)?
> What follows is a transcript of my attempt. I did not include the
> results of the final Solve command because it is page upon pages.
> While I have to hand it to Mathematica for being able to produce such
> a fantastical answer, my gut tells me the "real" answer is nowhere
> near this complicated. I tried applying Simplify to the result, but I
> ran out of memory before it could finish (I have 1GB physical, 2GB
> virtual).
> Can anyone provide advice on how to coax Mathematica into a simpler
> answer? (Or confirm that the answer really is this complicated.)
> Thanks,
> J. Kennedy

Your equations are equivalent to finding the roots of a 4th order 
polynomial. By default, Solve will use radicals to express the root of a 
  quartic polynomial, and hence yields an enormous and useless result in 
this case. We can control this behavior by changing the options of Roots:

SetOptions[Roots, Cubics -> False, Quartics -> False]

Now using Solve will produce much smaller results at the cost of 
containing explicit Root objects.

Carl Woll
Wolfram Research

> In[1]:=
> $Version
> Out[1]=
> "5.1 for Microsoft Windows (October 25, 2004)"
> In[3]:=
> e = x^2/a^2 + y^2/b^2 == 1
> Out[3]=
> x^2/a^2 + y^2/b^2 == 1
> In[4]:=
> x0 = a*Cos[\[Theta]]
> Out[4]=
> a*Cos[\[Theta]]
> In[5]:=
> y0 = b*Sin[\[Theta]]
> Out[5]=
> b*Sin[\[Theta]]
> In[6]:=
> c = (x - x0)^2 + (y - y0)^2 == r^2
> Out[6]=
> (x - a*Cos[\[Theta]])^2 + (y - b*Sin[\[Theta]])^2 == r^2
> In[7]:=
> Solve[{e, c}, {x, y}]
> [thousand of lines deleted]

  • Prev by Date: Re: variable "K"? (Really strange behavior . . . )
  • Next by Date: Re: Re: variable "K"? (Really strange behavior . . . )
  • Previous by thread: the ellipse and the circle
  • Next by thread: Re: the ellipse and the circle