function Mandelbrotov_skup(brTc, brIt)

% brIt - broj iteracija u Mandelbrotovom skupu
% brTc - broj tacaka unutar slike
% test:
%           Mandelbrotov_skup(1000,51)
% Generisemo z = 0 sa njegovim realnim i imaginarnim delom
% z(n+1) = z(n)^2 + c

clf;
zRe = zeros(brTc,brTc);
zIm = zeros(brTc,brTc);
% Generisemo konstantu c sa njenim realnim i imaginarnim delom
cRe = repmat(linspace(-1.5,0.5,brTc),brTc,1);
cIm = repmat(linspace(-1,1,brTc)',1,brTc);

% Iterativna formula 
for j = 1 : brIt
    % Racunamo q = z*z + c u kompleksnom prosotoru
    % q je konanca promenljiva u koju pakujemo rezultate
    qRe = zRe .* zRe - zIm .* zIm + cRe;
    qIm = 2 .* zRe .* zIm + cIm;
    % Da bismo izbegli numericku divergenciju ogranicicemo  q na [-5 5]    
    zRe = qRe;
    q_veci_od_pet =  find(qRe > 5.);
    zRe(q_veci_od_pet) = 5.;
    q_manji_od_pet = find(qRe < -5.);
    zRe(q_manji_od_pet) = -5.;
   
    zIm = qIm;
    z_veci_od_pet = find(qIm > 5.);
    zIm(z_veci_od_pet) = 5.;
    z_manji_od_pet = find(qIm < -5.);
    zIm(z_manji_od_pet) = -5.;
end

% Graficki prikazujemo
slika = log( sqrt(zRe.*zRe+zIm.*zIm) + 1);
imagesc(slika); % crta sliku sa koloritom