monocromo.m

x=-2:0.002:1; %porción del eje real tomada
y=-1:0.002:1; %porción del eje imaginario tomada
itM=100; %"infinito"
[X,Y]=meshgrid(x,y); 
A=X+1j*Y; %la matriz A es la retícula de puntos que vamos a pintar
D=zeros(size(A)); %creamos D, la matriz que informa sobre A(i,j) está o no en el conjunto: D(i,j) vale 0 si no está y 1 si sí está
  
%% Averiguamos si los números están en el conjunto
for i=1:length(y)
    for j=1:length(x)
        z0=A(i,j); %tomamos el número que vamos a analizar en el bucle de abajo
        z=0; %inicializamos la variable
        c=0; %inicializamos el contador de iteraciones
        while abs(z)<2
            z=z^2+z0;
            c=c+1; 
            if c==itM %es decir, el número pertenece al conjunto
                D(i,j)=1;
                break
            end
        end
    end     
end
                  
%% Coloreamos la matriz D
 negro=[0,0,0];
 blanco=[1,1,1];
 map=[blanco %mapa de color, lo explicaremos más abajo
     negro];
imagesc(D)
colormap(map)
axis equal %ejes proporcionados
saveas(gcf,'efidelbrot.jpg') %para guardar la imagen