%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%         10.gyakorlat        %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Nemlin. din. rsz-ek alapjai %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%        2014.12.03/04        %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%         Box-counting        %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

clear all;
figure
xlim([0 1000])              % x és y tengelyen meddig ábrázolunk
ylim([0 1000])              % de! x értéke -0.5 és 0.5 közötti!
hold on;

Accuracy = 0.01;

x = -0.5: 0.001: 0.5;
y = x';
[X,Y] = meshgrid(x,y);

z0 = X + 1i*Y;              % a komplex síkon vagyunk
original=1.*z0;
depth = 20;
for k = 1:depth
	z0 = z0 - (z0.*z0.*z0 + 1)./(3.*z0.*z0);
end
% más-más értéket/színt (20,40,60) adok a pontoknak attól függõen, hogy melyik megoldáshoz jut az adott pontból indítva a rendszer 
Red=(abs(z0 - (0.5000 - 0.8660i))<Accuracy)*20;
Green=(abs(z0 - (0.5000 + 0.8660i))<Accuracy)*40;
Blue=(abs(z0 - -1)<Accuracy)*60;
%ami nem konvergál adott iteráció szám alatt az marad 0:= sötétkék
image(Red+Green+Blue)   %a 3 db szín mátrixot összegzem és ábrázolom

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


%% összes sötétkék megszámolása:

% osszesitett = ;
% szamlalo_1 = 0;
% ;

%%

% % összes sötétkék megszámolása, kicsit másképp
% szamlalo = 0;
% for k = 1:1:1001
%     for j = 1:1:1001
%     
%         if abs(z0(k,j) - (0.5000 - 0.8660i))<Accuracy
%             
%         elseif abs(z0(k,j) - (0.5000 + 0.8660i))<Accuracy
%             
%         elseif abs(z0(k,j) - - 1)<Accuracy
%         
%         else
%             szamlalo = szamlalo+1;
%         end
%     end
% end

%%

%n = ;
%epsz = ;

%% egyenesek kirajzolása - mennyire osztjuk fel
% for i = 0:epsz:1000
%     i = i+0.5;
%     plot([i,i],[0,1001]);
%     plot([0,1001],[i,i]);   
% end

%%

% szamlalo_1 = 0;
% osszesitett = Red+Green+Blue;
% for a = 1:epsz:(1001-epsz)
%     for b = 1:epsz:(1001-epsz)
%         osszesitett_resz = osszesitett(a:(epsz+a-1), b:(epsz+b-1));
%         
%         reszeredmeny = sum(sum(osszesitett_resz == 0));
%         
%         if reszeredmeny > 0
%             szamlalo_1 = szamlalo_1 +1;
%                    
%         end
%         
%     end
% end

%% csak, hogy lássuk az eredményeket:
%% szamlalo (ha a másik számlálást nézzük)
szamlalo_1

% dimenzio meghatározása:

%dimenzio= 


