1#/usr/bin/python 2 3from __future__ import division; 4import sys 5 6rcorn = -2.0; 7icorn = -2.0; 8side = 4.0; 9iheight = 500; 10iwidth = 1000; 11maxits = 100; 12 13for ii in range(0, iheight+1): 14 for ir in range(0, iwidth+1): 15 cr = rcorn + (ir/iwidth) * side; 16 ci = icorn + (ii/iheight) * side; 17 18 zr = 0.0; 19 zi = 0.0; 20 21 # z := z^2 + c 22 iti = 0; 23 escaped = False; 24 for iti in range(0, maxits): 25 mag = zr*zr + zi+zi; 26 if mag > 4.0: 27 escaped = True; 28 break; 29 zt = zr*zr - zi*zi + cr; 30 zi = 2*zr*zi + ci; 31 zr = zt; 32 if (escaped): 33 sys.stdout.write("o"); 34 else: 35 sys.stdout.write("."); 36 sys.stdout.write("\n"); 37