from math import log10, floor, isnan def dec_to_hex(num, byteWidth): result = [] while num > 0: pre = num % 16 num //= 16 cur = num % 16 num //= 16 result.append(16 * cur + pre) for n in range(byteWidth - len(result)): result.append(0) return result def hue_to_RGB(degree): degree %= 360 if 0 <= degree and degree < 60 or 300 <= degree and degree <= 360: r = 255 elif 60 <= degree and degree < 120: r = 4.25 * (120 - degree) elif 120 <= degree and degree < 240: r = 0 elif 240 <= degree and degree < 300: r = 4.25 * (degree - 240) if 0 <= degree and degree < 60: g = 4.25 * degree elif 60 <= degree and degree < 180: g = 255 elif 180 <= degree and degree < 240: g = 4.25 * (240 - degree) elif 240 <= degree and degree <= 360: g = 0 if 0 <= degree and degree < 120: b = 0 elif 120 <= degree and degree < 180: b = 4.25 * (degree - 120) elif 180 <= degree and degree < 300: b = 255 elif 300 <= degree and degree <= 360: b = 4.25 * (360 - degree) return [int(r), int(g), int(b)] name = input('Name: ') function = input('f(z, c) = ') initial = float(input('z_0 = ')) ict = int(input('Iteration count threshold = ')) mt = float(input('Magnitude threshold = ')) pixels_per_unit = int(input('Pixels per unit = ')) real_lower_bound = float(input('Real lower bound = ')) real_upper_bound = float(input('Real upper bound = ')) imaginary_lower_bound = float(input('Imaginary lower bound = ')) imaginary_upper_bound = float(input('Imaginary upper bound = ')) colored = True if input('Colored? (yes or no) ').lower() == 'yes' else False width = int((real_upper_bound - real_lower_bound) * pixels_per_unit) height = int((imaginary_upper_bound - imaginary_lower_bound) * pixels_per_unit) size = 3 * width * height + 26 print('Generating the image of the fractal...') f = open(name + '.bmp', 'wb') f.write(bytes([66, 77,] + dec_to_hex(size, 4) + [0, 0, 0, 0, 26, 0, 0, 0, 12, 0, 0, 0,] + dec_to_hex(width, 2) + dec_to_hex(height, 2) + [1, 0, 24, 0,])) for y in range(height): for x in range(width): c = complex(x / pixels_per_unit + real_lower_bound, y / pixels_per_unit - imaginary_upper_bound) z = initial for n in range(ict): try: z = eval(function) except: f.write(bytes([255, 255, 255,])) break else: if abs(z) > mt: rgb = hue_to_RGB(n * 6) if colored else [255, 255, 255,] f.write(bytes([rgb[2], rgb[1], rgb[0],])) break else: f.write(bytes([0, 0, 0,])) f.close() print('Done')

Generate Images of Fractals with Python (best in IDLE)