jueves, 23 de enero de 2014

Minimo fractal de Mandelbrot

Mandelbrot minimo
Mandelbrot minimo

Hace unos meses que he tenido que retomar el uso de octave y recordé que hace unos años hice un pequeño script que calcula en blanco y negro el fractal de mandelbrot en 17 lineas de octave, incluyendo borrado y mostrado en pantalla y pausa pertinente al final del script. Se puede reducir mas el tamaño del código quitando por ejemplo las variable y poniendo directamente los valores pero el núcleo del script crea el fractal en unas 10-11 lineas. ¿Seria posible hacerlo en menos?:

#!/usr/bin/octave
# laantiguafragua.blogspot.com
clf;
puntos=400;
puntosx=linspace(-2.1,0.9,puntos);
puntosy=linspace(-1.5,1.5,puntos);
[X, Y]=meshgrid(puntosx, puntosy);
c=X+Y*i;
Z=zeros(puntos);
itr=100; # Numero de iteraciones
for k=1:itr
Z=(Z.^2)+c;
end;
W=exp(-abs(Z));
[J, map]=gray2ind(W);
imshow(J, map);
pause;