TOPICS
Search

Devil's Curve


DevilsCurve

The devil's curve was studied by G. Cramer in 1750 and Lacroix in 1810 (MacTutor Archive). It appeared in Nouvelles Annales in 1858. The Cartesian equation is

 y^4-a^2y^2=x^4-b^2x^2,
(1)

equivalent to

 y^2(y^2-a^2)=x^2(x^2-b^2),
(2)

the polar equation is

 r^2(sin^2theta-cos^2theta)=a^2sin^2theta-b^2cos^2theta,
(3)

and the parametric equations are

x=costsqrt((a^2sin^2t-b^2cos^2t)/(sin^2t-cos^2t))
(4)
y=sintsqrt((a^2sin^2t-b^2cos^2t)/(sin^2t-cos^2t)).
(5)

The curve illustrated above corresponds to parameters a^2=1 and b^2=2.

It has a crunode at the origin.

Devil's curve animation

For 0<|a|<|b|, the central hourglass is vertical, while for |a|>|b|>0, it is horizontal. When |a|=|b|>0, the equation factors as

 (y^2-x^2)(x^2+y^2-a^2)=0,
(6)

so the degenerate curve is a circle together with the two diagonal lines y=+/-x.

Let M=max(|a|,|b|) and m=min(|a|,|b|), with 0<m<M. The two bounded lobes have total area

 A=(M^2+m^2)tan^(-1)(m/M)-(M^2-m^2)tanh^(-1)(m/M).
(7)

Each lobe has area A/2. The other real branches are unbounded and are not included. As m->M, this expression has limit (piM^2)/2, which is the limit of the two selected lobes rather than the area of the entire disk. (E. Weisstein, Sep. 23, 2026).

The total perimeter L of the same two bounded lobes can be written using a six-variable Lauricella function. Set q=(m^2)/(M^2) and

 P_q(v)=q+(1-4q-q^2)v+3q(1+q^2)v^2+(-q^2-4q^3+q^4)v^3+q^3v^4,
(8)

and define z_1,...,z_4, without choosing an ordering, by

 (P_q(v))/q=product_(j=1)^4(1-z_jv).
(9)

Then

 L=2pimsqrt(q)F_D^((6))(1/2;-1/2,-1/2,-1/2,-1/2,3/2,1;1;z_1,z_2,z_3,z_4,q,-q).
(10)

(E. Weisstein, Sep. 23, 2026). Here the Lauricella function is taken on the explicit Euler integral branch

 F_D^((n))(alpha,beta_1,...,beta_n;gamma;w_1,...,w_n)=(Gamma(gamma))/(Gamma(alpha)Gamma(gamma-alpha))int_0^1v^(alpha-1)(1-v)^(gamma-alpha-1)product_(j=1)^n(1-w_jv)^(-beta_j)dv.
(11)

For the expression for L, the product of the first four square root factors is chosen to be the positive real value sqrt((P_q(v))/q) on 0<=v<=1. The remaining factors 1-qv and 1+qv are positive there, and P_q(v)>0. Thus this branch gives a finite real length for the bounded lobes only, not for the unbounded branches.

ElectricMotor

A special case of the Devil's curve is the so-called "electric motor curve":

 y^2(y^2-96)=x^2(x^2-100)
(12)

(Cundy and Rollett 1989).


See also

Area, Barbell Graph, Butterfly Curve, Dumbbell Curve, Eight Curve, Lauricella Functions, Lemniscate, Perimeter, Piriform Curve, Pitchfork Bifurcation, Teardrop Curve

Explore with Wolfram|Alpha

References

--. Nouvelle Annales, p. 317, 1858.Cramer, G. Introduction a l'analyse des lignes courbes algébriques. Geneva, Switzerland, p. 19, 1750.Cundy, H. and Rollett, A. Mathematical Models, 3rd ed. Stradbroke, England: Tarquin Pub., p. 71, 1989.Gray, A. Modern Differential Geometry of Curves and Surfaces with Mathematica, 2nd ed. Boca Raton, FL: CRC Press, pp. 92-93, 1997.Lacroix, S. F. Traité du calcul différentiel et intégral, Vol. 1. Paris, France, p. 391, 1810.Lawrence, J. D. A Catalog of Special Plane Curves. New York: Dover, pp. 151-152, 1972.MacTutor History of Mathematics Archive. "Devil's Curve." https://mathshistory.st-andrews.ac.uk/Curves/Devils/.Smith, D. E. History of Mathematics, Vol. 2: Special Topics of Elementary Mathematics. New York: Dover, p. 328, 1958.

Referenced on Wolfram|Alpha

Devil's Curve

Cite this as:

Weisstein, Eric W. "Devil's Curve." From MathWorld--A Wolfram Resource. https://mathworld.wolfram.com/DevilsCurve.html

Subject classifications