clear all; close all; clc;

maxIterations = 249;
divergenceLimit = 4;

resX = 2500;
resY = 2500;

startReal = -0.555;
startImag = -0.515;
endReal = -0.535;
endImag = -0.485;

stepsizeX = (abs(startReal) - abs(endReal))/resX;
stepsizeY = (abs(startImag) - abs(endImag))/resY;

canvas3d = zeros(resY,resX,3);

kolor = colormap(hot(250));

finished = false;
for x = 1 : resX
    for y = 1 : resY
        c = complex(startReal+stepsizeX*x, startImag+stepsizeY*y);
        z = 0.0;
        iterationCount = 0;
        finished = false;
        while(~finished)
            iterationCount = iterationCount + 1;
            z = z^2 + c;
            if(iterationCount > maxIterations) 
                finished = true;
                canvas3d(y,x,:) = 0;
            end    
            if(abs(z) > divergenceLimit)
                finished = true;
                canvas3d(y,x,:) = kolor(iterationCount,:);
            end        
        end
    end
end

h = image(canvas3d);%canvas./max(canvas(:))*350)

axis off;
axis image;
imwrite(canvas3d, './test.tiff', 'tif');