A Mandelbrot Fractal in Python
I coded up this Mandelbrot fractal in python while watching TV. Not sure if is helpful for anybody but you may want to take a look at it and enjoy the bauty of chaotic dynamical systems. The code is posted here and below too. Here is a fascinating high quality version of it.
# Mandelbrot set # By Mark Alen # linux_jvm@yahoo.com # April 2012 import ImageDraw from PIL import Image, ImageFilter from math import log white = (255, 255, 255) width = 5000 height = width image1 = Image.new("RGB", (width, height), white) draw = ImageDraw.Draw(image1) # http://en.wikipedia.org/wiki/Mandelbrot_set for xpix in range(1,width+1): for ypix in range(1,height+1): x0 = (xpix*1.0/width*3.5) -2.5 y0 = (ypix*1.0/height*2)-1 x = 0 y = 0 iteration = 0 max_iteration = 1000 while ( (x*x + y*y) < 4) & (iteration < max_iteration ): xtemp = x*x - y*y + x0 y = 2*x*y + y0 x = xtemp iteration = iteration + 1 mycol =int(255.0*(1-log(iteration/1000.0*255+1)/log(256))) color = (255-mycol,255-mycol,255-mycol) if iteration == max_iteration: color = white draw.point((xpix,ypix), color) if ((xpix*(100/1)%(width))==0): print (xpix*ypix*1.0)/((width*height))*100.0, "%" filename = "mandel.png" image1.save(filename) print "Done!"