Creating Mandelbrot in Python

Mandelbrot sets are fascinating mathematical objects that exhibit intricate and infinitely complex patterns. They are named after Benoit Mandelbrot, a mathematician who discovered and popularized them in the 1970s. In this article, we will explore the concept of Mandelbrot’s sets, their scientific significance, and how to generate them using Python. We will start with a basic implementation and gradually enhance it to create visually appealing and dynamically colored Mandelbrot’s images.

Creating Mandelbrot in Python

Understanding Mandelbrot Sets:

  • Mandelbrot sets are created by iterating a mathematical function on complex numbers.
  • Each point in the complex plane is tested to determine whether it remains bounded or escapes to infinity.
  • The boundary between bounded and unbounded points forms intricate fractal patterns.

Scientific Significance:

  • These sets have profound implications in various scientific fields, including mathematics, physics, and computer graphics.
  • They exhibit self-similarity, meaning they exhibit similar structures at different levels of magnification.
  • The fractal nature of Mandelbrot sets provides insights into chaos theory, dynamical systems, and complex analysis.

Purpose of Generating Mandelbrot Sets:

  • Exploring and visualizing the beauty of mathematical concepts.
  • Illustrating the power of computational tools in generating complex structures.
  • Uncovering the underlying patterns and symmetries in the natural world.

Mathematical Explanation:

Mandelbrot sets are generated by iterating a mathematical function on complex numbers. The formula used to calculate each point within the set is:

Zₙ = Zₙ₋₁² + C

Here, Zₙ is the complex number being iterated, Z₀ starts at 0, and C is the initial complex number for each point on the complex plane. The iteration is repeated for each point until either the absolute value of Zₙ exceeds 2 or a maximum number of iterations is reached.

To calculate the color or shade of each pixel in the Mandelbrot image, we track the number of iterations it takes for the absolute value of Zₙ to exceed 2. This iteration count provides a measure of how quickly or slowly each point escapes to infinity. Points that remain bounded within a certain threshold are considered part of the Mandelbrot’s set.

Steps to generate Pixels

Let’s break down the process step by step:

  1. Define the complex number C for each point on the complex plane. The position of each pixel in the image corresponds to a unique C value.
  2. Set Z to 0, as Z₀ = 0.
  3. Iterate the formula Zₙ = Zₙ₋₁² + C for a certain number of iterations or until the absolute value of Zₙ exceeds 2.
  4. Check the absolute value of Zₙ after each iteration. If it exceeds 2, the point is likely to escape to infinity and is not part of the Mandelbrot set. Exit the loop.
  5. If the maximum number of iterations is reached without the absolute value of Zₙ exceeding 2, the point is considered part of the Mandelbrot set. The number of iterations determines the color or shade of the corresponding pixel.

By repeating this process for every pixel on the complex plane, we generate a visual representation of the Mandelbrot set, where different colors or shades represent different escape times.

Explanation

It’s important to note that the complex plane corresponds to the x and y coordinates of the image, and the rectangular region of the plane displayed in the image can be adjusted by setting the minimum and maximum values for the x and y axes. This allows for zooming and exploring specific regions of the Mandelbrot set in greater detail.

The mathematical beauty of the Mandelbrot set lies in its intricate and infinitely complex patterns, which arise from the interaction between the real and imaginary components of complex numbers during the iteration process. Exploring and visualizing these patterns is not only aesthetically captivating but also provides insights into the principles of chaos theory, fractals, and complex analysis.

Generating Mandelbrot Sets with Python:

To start, let’s write a Python program that generates a basic black and white Mandelbrot set image. We will use the numpy library for numerical operations and the PIL (Python Imaging Library) library for image creation. Here’s the code:

   import numpy as np
   from PIL import Image

   def mandelbrot(c, max_iter):
       z = 0
       for i in range(max_iter):
           z = z * z + c
           if abs(z) > 2:
               return i
       return max_iter

   def create_mandelbrot(width, height, xmin, xmax, ymin, ymax, max_iter):
       image = Image.new('L', (width, height))
       pixels = image.load()
       for x in range(width):
           for y in range(height):
               zx = np.interp(x, [0, width], [xmin, xmax])
               zy = np.interp(y, [0, height], [ymin, ymax])
               c = zx + zy * 1j
               color = mandelbrot(c, max_iter)
               pixels[x, y] = color

       return image

   # Example usage
   width = 800
   height = 600
   xmin, xmax = -2.0, 1.0
   ymin, ymax = -1.5, 1.5
   max_iter = 256

   image = create_mandelbrot(width, height, xmin, xmax, ymin, ymax, max_iter)
   image.show()

Output of Code:

In this code, we define two functions: mandelbrot(c, max_iter) and create_mandelbrot(width, height, xmin, xmax, ymin, ymax, max_iter). The mandelbrot() function calculates the iteration count for a given complex number c using the Mandelbrot algorithm. The create_mandelbrot() function generates the Mandelbrot image by iterating over each pixel, mapping it to the corresponding complex number, and assigning colors based on the iteration count.

You can modify the width, height, xmin, xmax, ymin, ymax, and max_iter variables to adjust the size, position, and level of detail of the generated Mandelbrot image.

This code generates a basic Mandelbrot image in black and white base based on the iteration count.

Enhancing the Mandelbrot Image:

While the black and white Mandelbrot image is fascinating, let’s enhance it by adding random colors to each point. This will create visually appealing Mandelbrot sets with vibrant hues. Here’s the updated code snippet:

   import numpy as np
   from PIL import Image
   import random

   def mandelbrot(c, max_iter):
       z = 0
       for i in range(max_iter):
           z = z * z + c
           if abs(z) > 2:
               return i
       return max_iter

   def create_mandelbrot(width, height, xmin, xmax, ymin, ymax, max_iter):
       image = Image.new('RGB', (width, height))
       pixels = image.load()

       base_color = (random.randint(0, 255), random.randint(0, 255), random.randint(0, 255))

       for x in range(width):
           for y in range(height):
               zx = np.interp(x, [0, width], [xmin, xmax])
               zy = np.interp(y, [0, height], [ymin, ymax])
               c = zx + zy * 1j
               color = mandelbrot(c, max_iter)
               shade = int((color / max_iter) * 255)
               r = int(base_color[0] * shade / 255)
               g = int(base_color[1] * shade / 255)
               b = int(base_color[2] * shade / 255)
               pixels[x, y] = (r, g, b)

       return image

   # Example usage
   width = 800
   height = 600
   xmin, xmax = -2.0, 1.0
   ymin, ymax = -1.5, 1.5
   max_iter = 256

   image = create_mandelbrot(width, height, xmin, xmax, ymin, ymax, max_iter)
   image.show()

In this updated code, we introduce a base_color variable that contains a random RGB color tuple. The shade of the color is calculated based on the iteration count, similar to the previous version. By adjusting the RGB values with the shade, we create different shades of the base color for each pixel, resulting in visually distinct Mandelbrot sets.

output:

Formula of Mandelbrot Sets:

The formula used to calculate the Mandelbrot set is z = z^2 + c, where z is the complex number being iterated and c is the initial complex number for each point on the complex plane. This iteration is repeated until either the absolute value of z exceeds 2 or the maximum number of iterations is reached. The iteration count determines the color or shade of the corresponding pixel in the Mandelbrot image.

Conclusion:

Creating Mandelbrots sets in Python offers a fascinating journey into the realm of fractals, complexity, and computational art. Through this article, we have explored the scientific significance of Mandelbrot sets and provided a step-by-step guide to generating and enhancing Mandelbrot images using Python. By manipulating color, zooming, and angles, we can unlock a vast range of visual possibilities, allowing us to appreciate the beauty and complexity hidden within these mathematical wonders. Have fun exploring the infinite depths of the Mandelbrot set and let your creativity soar!

Your Task:

Remember, the true beauty of Mandelbrot sets lies in their infinite intricacy and self-similarity, so feel free to experiment with different parameters, zoom levels, and color schemes to discover new and captivating images. Enjoy the journey! Share your Results with us in comment section or on our Facebook Group

Leave a Comment