Esercizi e sorgenti dei programmi
Visualizzazione dell'insieme di Mandelbrot in linguaggio Python
Di seguito riportiamo la codifica di un programmino Python che, utilizzando la libreria Graphics visualizza l'insieme di Mandelbrot. Il programmino è stato sviluppato come esercizio in laboratorio.
Codifica in linguaggio Python 3.x
from graphics import * import numpy as np def punto(re, im, zre0, zim0, maxCicli): zre = zre0 zim = zim0 k = 0 while zre**2 + zim**2 <=4.0 and k < maxCicli: temp = zre**2 - zim**2 +re zim = 2.0*(zim*zre) + im zre = temp k = k+1 return k win = GraphWin("Mandelbrot", 800, 600, autoflush=False) win.setCoords(-2, -1.0, 0.7, 1.0) max = 100 deltaRe = float(2.7/800.0) deltaIm = float(2.0/600.0) cre = -2 while cre < 0.7: cim = 2 while cim > -1: c = punto(cre, cim, 0, 0, max) if c >= max: win.plot(cre,cim,color_rgb(0,0,0)) else: win.plot(cre,cim,color_rgb(int(255*c/max), 0, 255-int(255*c/max))) cim = cim-deltaIm cre = cre+deltaRe update() win.getMouse() win.close()
Il programma produce il seguente output: