map = [1 1 1 
0 0 1
1 0 0
0 0 0];
image_size = [1000 1000];   % Image size
limits = [-1.5 1.5 -1 1];   % Real and imaginary limits
c = -0.4 + 0.6i;            % Fixed complex number
total_iterations = 300;
hold on
colour = julia(c, total_iterations, image_size, limits);
plot(-0.629372,0.735756,'g*')
I = imagesc(colour);
colormap(map)
myaa('publish')
hold off