Implicit form for a logarithmic spiral

408 Views Asked by At

I was wondering if there is an implicit form for a logarithmic spiral. For example, if

$$ x=e^{-t}\cos(t)\\y=e^{-t}\sin(t)$$

we can write $x^2+y^2=e^{-2t}$ and $y/x=\tan(t)$ which yields $$x^2+y^2=e^{-2\arctan(y/x)}$$ for $x\neq0$, but of course I would like to find one defined for any $(x,y)\in\mathbb{R}^2$. Any help?

2

There are 2 best solutions below

5
On

Rearranging the given equation gives the new equation $$y=x\tan{\bigg(\ln{\bigg(\frac{1}{\sqrt{x^2+y^2}}\bigg)}\bigg)}$$ which is defined for any $(x,y)\in\mathbb{R}^2$.

0
On

The implicit equation you give

$$x^2+y^2=e^{-2\arctan(y/x)}\tag{1}$$

is indeed not defined for $x\neq0$ but has the worst drawback that it generates a double spiral (the same issue met by the "solution" given by Peter Foreman) as shown here :

enter image description here

Fig. 1 : Formulation (1) provides a double spiral instead of a single one.

Nevertheless, there is a solution, suggested by @Blue (who so often has given very interesting remarks) : if we accept to consider "arctan2" as a function "like others", it suffices to transform (1) into

$$x^2+y^2=\exp(-2 \ \text{arctan2}(y,x))\tag{2}$$

(please note the arguments' order for atan2 : $y$, then $x$), or the equivalent form :

$$2 \ \text{arctan2}(y,x)+\log(x^2+y^2)=0\tag{3}$$

Another view of the issue : the spiral can be viewed as the projection of the intersection line of two surfaces as shown on Fig. 2 :

enter image description here

Fig. 2 : The logarithmic spiral as the vertical projection of the intersection line of surfaces with equations $z=$ln$(x^2+y^2)$ (a pseudo-cone) and $z=-2 \arctan \! 2(y,x)$

Remark : "arctan2", where the "2" means "with 2 arguments", finds its origin in numerical analysis (see for example https://en.wikipedia.org/wiki/Atan2 which traces it back to Fortran IV, i.e., around 1960). It is nowadays present in all scientific software, for example atan2 in Matlab, ArcTan[Y,X] with Mathematica... See its particular behavior along the negative $x$-axis "cut" in https://uk.mathworks.com/help/matlab/ref/atan2.html.