## Friday, 6 April 2007

### Essay Backtrack - Julia Sets

I have been mentioning Julia Sets in my essay updates a lot, and yet I still have not really explained what they are in great detail, or given you any examples. This post aims to explain Julia Sets and some related results in a manner which is graspable by a first year mathematician or so. The prerequisite knowledge really for a basic understanding is Complex Numbers and Functions.

To find a Julia Set, we first need the function for which it is a Julia Set of. Although this concept can be extended to other cases (or reduced to in the case of $\mathbb{R}$), the most common case is that of a function $f:\mathbb{C}\rightarrow\mathbb{C}$. The function $f$ also needs to be quite smooth, specifically it needs to be holomorphic (you do not really need to understand this greatly for a basic understanding, but if you're interested wikipedia has the definition behind the link).

Now, the concept of a Julia Set is based upon the long term effect of repeated iterations of $f$ on points in the plane. That is to say, given a point $z\in\mathbb{C}$, we want to find out what happens as we take repeated iterations of $z$, that is $f(z),\,f^2(z)=f(f(z)),\ldots f^n(z)\ldots$ as $n\rightarrow\infty$. There are, in a simplistic case, 2 different things that can happen. The first is that the point tends off to infinity upon repeated iteration of $f$. The second is that it doesn't, and remains bounded. Now, this may not sound particularly exciting, but this concept can give some surprising results right from the offset.

Generally speaking, a good point to start at in the study of Julia Sets is looking at the Julia Set of a Polynomial Map - $f_c(z)=z^2+c$ (note that it is a simple exercise to notice that in fact any complex polynomial map of the form $ax^2+bx+c$ can be reduced to one of the form $z^2+c$). So, let us look at the most simple of these maps, that is $f_0(z)=z^2$. Now, let us try to find which points are in the set $\{z:\left\vert f^n(z)\right\vert. In this particular case, it is simple. As we know, if we take a point $re^{i\theta}\in\mathbb{C}$, then $(re^{i\theta})^n=r^ne^{in\theta}$, and so all points such that $r<1$ will move towards $0$, all points with $r>1$ will move toward infinity and all points with $r=0$ will stay on the unit circle. We can now say that the filled julia set $\mathcal{K}=\{re^{i\theta}:r\leq1\}$. Furthermore, the Julia set is $\mathcal{J}=\partial\mathcal{K}=\{re^{i\theta}:r=1\}$. Finally, in terms of definitions, the Fatou Set of this function is defined to be $\mathcal{F}=\mathbb{C}\setminus\mathcal{K}$, that is all the points which are not in the filled Julia Set.

So, having seen a very simple example, we may wonder what will the effect be when we change our value of $c$? Well, this is where we get our surprising result. Running a simple script through MATLAB (which I will write about in the future), we find that our Filled Julia set (with $c=0.52+0.25i$) now looks like this:

This is a set of fractal nature, that is to say that it has an infinitely fine structure (so if you zoom in on the edge of it, it remains detailed at arbitrary levels) and small parts of it bear a resemblance to the entire set.

Julia Sets were some of the first things I came across when writing my second year essay, and my interest in them has really pushed me towards learning about complex dynamics in general, thus my third year essay!

If you want to find out more about Julia Sets, and in fact fractals in general, I would suggest that the book Fractal Geometry: Mathematical Foundations and Applications by Kenneth Falconer is brilliant, and well worth a read for anyone with interest in this area.

Stephen said...

Hi

I liked this little article on Julia Sets, and am looking forward to your MATLAB code on creating such a Julia Set. I am neither a mathmetician, nor programmer (geophysicist actually) but find this subject extrememly interesting.

Thanks a lot

Steve

P.Bhattacharyya said...

Hello,
Good to see your notes on Julia set. Would be happy to get your MATLAB Codes on Julia set.Thank you.

It is interesting to point out that at the beginning, the definitions of Julia and Fatou sets were reversed. What is nowadays called the Julia set J(f) was called the Fatou set F(f) and what is nowadays called the Fatou set was called the Julia set.(see for example the papers of I.N.Baker,J.Myrberg,P.Bhattacharyya etc published in the 60s, and till late 70s.It is historical now.

The interest in the theory was revived after high speed computing facilities became easily available.

Evgenia said...

Hey,
thanks for your notes for julia sets!
they are really nice.
evgenia