Множество Жюлиа

iter = 30;  % число итераций
npix = 800; % число точек по каждой из осей координат

% область
dl = 1.25;
x = linspace(-dl,+dl,npix);
y = linspace(-dl,+dl,npix);

% массив с точками исследуемой области
[X,Y] = meshgrid(x,y);
Z = X + i*Y;

% начальнаЯ точка
C = 0.27334-0.00742*i;

for l = 1:iter
    Z = Z.^2 + C;       
end;

W = exp(-abs(Z));

imagesc(W);

colormap(prism(1024));
axis ('equal','off')