function Julijini_skupovi(c,k,Xr,Yr)
% 
% Crtamo Julijine skupove za dato c
%
% c - parametar jednacine z = z^2 + c
% k - broj iteracija
% Xr - funkcija se crta na intervalu [Xr(1,1) Xr(1,2)] 
% Yr - funkcija se crta na intervalu [Yr(1,1) Yr(1,2)]
%
% test:
% Julijini_skupovi(i,100,[-2 2],[-2 2])

clf;

n = 400;

x = linspace(Xr(1,1),Xr(1,2),n);
y = linspace(Yr(1,1),Yr(1,2),n);

% pravimo mrezu po kojoj ce se crtati skup
[X,Y] = meshgrid(x,y); 

% M je matrica sa podacima iz Julijinih skupova
M = zeros(length(X),length(Y));

for l = 1 : size(X,2)
    for j = 1 : size(Y,2)
       [w Iter] = Julia(X(l,j)+Y(l,j)*i , c, k);
       M(l,j) = M(l,j) + Iter;
    end
end
hold on;
%axis off;

colormap(jet);
pcolor(M);
shading interp;

hold off;

end
%**********************%

function [pri,it] = Julia(z,c,k)

R = max(abs(c),2);
i = 0;
while i < k
    if abs(z) > R
        pri = 1;
        it = i;
        return;
    end   
    z = z^2 + c;
    i = i + 1;
end
pri = 0;
it = i;
end