C
Almost all points of the complex plane are associated with one of the roots of a given polynomial in the following way: the point is used as starting value for Newton's iteration, yielding a sequence of points If the sequence converges to the root, then was an element of the region . However, for every polynomial of degree at least 2 there are points for which the Newton iteration does not converge to any root: examples are the boundaries of the basins of attraction of the various roots. There are even polynomials for which open sets of starting points fail to converge to any root: a simple example is, where some points are attracted by the cycle rather than by a root.
An open set for which the iterations converge towards a given root or cycle (that is not a fixed point), is a Fatou set for the iteration. The complementary set to the union of all these, is the Julia set. The Fatou sets have common boundary, namely the Julia set. Therefore, each point of the Julia set is a point of accumulation for each of the Fatou sets. It is this property that causes the fractal structure of the Julia set (when the degree of the polynomial is larger than 2).
To plot images of the fractal, one may first choose a specified number of complex points and compute the coefficients of the polynomial
d-1 | |
p(z)=z | |
1z |
+ … +pd-1z+pd:=(z-\zeta1)(z-\zeta2) … (z-\zetad)
zmn=z00+m\Deltax+in\Deltay; m=0,\ldots,M-1; n=0,\ldots,N-1
C
A generalization of Newton's iteration is
zn+1=zn-a
p(zn) | |
p'(zn) |
where is any complex number.[1] The special choice corresponds to the Newton fractal. The fixed points of this map are stable when lies inside the disk of radius 1 centered at 1. When is outside this disk, the fixed points are locally unstable, however the map still exhibits a fractal structure in the sense of Julia set. If is a polynomial of degree, then the sequence is bounded provided that is inside a disk of radius centered at .
More generally, Newton's fractal is a special case of a Julia set.
Serie : Other fractals where potential and trigonometric functions are multiplied.
The Nova fractal invented in the mid 1990s by Paul Derbyshire,[2] [3] is a generalization of the Newton fractal with the addition of a value at each step:[4]
zn+1=zn-a
p(zn) | |
p'(zn) |
+c=G(a,c,z)
The "Julia" variant of the Nova fractal keeps constant over the image and initializes to the pixel coordinates. The "Mandelbrot" variant of the Nova fractal initializes to the pixel coordinates and sets to a critical point, where[5]
\partial | |
\partialz |
G(a,c,z)=0.
Commonly-used polynomials like or lead to a critical point at .
In order to implement the Newton fractal, it is necessary to have a starting function as well as its derivative function:
\begin{align}f(z)&=z3-1\ f'(z)&=3z2\end{align}
The three roots of the function are
z=1, -\tfrac12+\tfrac{\sqrt3}{2}i, -\tfrac12-\tfrac{\sqrt3}{2}i
The above-defined functions can be translated in pseudocode as follows:
//3*z^2float2 Derivative (float2 z)It is now just a matter of implementing the Newton method using the given functions.
For each pixel (x, y) on the target, do: