%**********************************************************
%**               pokusy s fraktály                      **
%**         konkrétne s juliovími mnozinami              **
%** Leos Marsalek (c) 2002-2003 Leos.Marsalek@tiscali.cz **      
%**********************************************************
tic;                        %spusteni vnitrnich hodin
obr=0;                      
c=-0.6699 + 0.3626i;        
%vyber nahodneho komlexniho cisla - na tomto cisla zalezi vzhled
                            %fraktalu
for x=1:320                 %krokovani a vypocet barvy pixelu
   
x_v=(x-160)/1000;        %krok pixelu
    
for y=1:240             
       y_v=(y-120)/1000;
       q=(x_v)+(y_v)*i;     
       
for b=1:512          % maximalni pocet interakci
           
q=(q)^b+c;       %vlastni vypocet fraktalu
           
if abs(q)>2
           obr(y,x)=b;          
           
break
           end
       end
   end
end
imageview(obr);             %vykresleni obrazku
toc;                        %zastaveni a vypsani delky vypoctu. Na Duronu 850 trva vypocet kolem 130s