To compute |φ(z)| at a point z we must simply follow the
orbit zk = fcok(z) until at a step
k we reach some point which is very far from 0 where
fc is close to fo (and zk
to wk). Then we make k itterations
fo-1 (in the "reverse" direction) i.e.
|φ(z)| = limk→∞ |f ok(z)|1/2k.
G(z) = limn→∞ ln|f on(z)|
This implies also that
G(f(z)) = 2 G(z).
It is evident that far from charged body equipotential curves are circles. Let R is a large circle then
G(f -ok(R)) = G(R)/2k = const
therefore preimages of the circle are equipotential curves for the Julia set (see the picture to the left for R = 10).
Electric field lines (orthogonal to equipotential curves) are called
external rays. They are preimages of the field lines of charged circle
r exp(2πit), r>1. Each such ray Rt
can be specified by its angle at infinity t.
Due to isomorphism φ(z) external rays dynamics under f
is the same as dynamics of the straight rays under z2. E.g.
f(Rt ) = R2t .
To the left every point is painted in the red, green and blue colors if for r > 10 its orbit falls into the sectors (0o, 120o), (120o, 240o), (240o, 360o) correspondigly. With enough patience you can see equipotential curves and orthogonal field lines. Exteral rays with rational angles t = 0, 1/4, 1/3, 1/2... are marked here.
At last to the left external rays p/12 for p = 0,2,..11
are plotted by means of the
Wolf Jung's algorithm.
As since I used only simplest Jung's magic this applet is a bit
buggy for large q.
If you need flexible and accurate program to explore external rays and much
more use his Mandel.
I am grateful to Adam Majewski for he convinced me to translate Wolf Jung's sources into Java.
|In particular, if we choose a constant Go > G(0), then G(z) = Go is a simple closed curve, canonically parametrized by the angle of the corresponding dynamic ray. In particular, the critical value z = fc(0) = c has a well defined external angle, which we denote by t(c). However, for Go = G(0) this points G(z) = G(0) make a figure eight curve (see to the left). The open set G(z) < G(0) splits as a disjoint union Uo and U1, where the Ub are the regions enclosed by the two lobes of this figure eight.|
According to Douady and Hubbard ψc is a conformal
isomorphism from the complement of M onto the complement of the closed
unit disk. Therefore it determines equipotential curves and external rays for
the Mandelbrot set.
It is known that an external ray whose angle is rational actually lands on
M. That is
 J. Milnor Dynamics in One Complex Variable: Introductory Lectures. Preprint ims90-5