% Show successive Julia Sets: 
% The parameters can be selected w>=14; p>=500.

disp(' ')
disp('The simulation shows several Julia Sets in the Figure above')
disp(' ')
disp('To Stop Demo, Type "CTRL+C" in the Command Window')
disp(' ')

c = 0.2+0.65i; 
w = 14; 
p = 500; 

JuliaSet(c,w,p);
title('c = 0.2+0.65i, Julia Set.')
pause(1)
clear all

c =-0.17+0.78i; 
w = 16; 
p = 500;

JuliaSet(c,w,p);
title('c = -0.17+0.78i, Julia Attractor.')
pause(1)

clear all

c = -0.12375+0.56508i; 
w = 16; 
p =500;

JuliaSet(c,w,p);
title('c= -0.12375+0.56508i, basin with an attractor.')
pause(1)
clear all

% Three attractors

c =-0.12+0.74i; 
w = 16; 
p = 500;

JuliaSet(c,w,p);
title('c = -0.12+0.74i.')
pause(1)
clear all

% Three attractor with thin spiral

c = -0.11+0.6557i; 
w = 16; 
p = 500;

JuliaSet(c,w,p);
title('c = -0.11+0.6557i)')
pause(1)
clear all

% Three attractors

c = -0.194+0.6557i; 
w = 16; 
p = 500;

JuliaSet(c,w,p);
title('c = -0.194+0.6557i, Cantor set.') 
pause(1)
clear all

% Dendrites

c = 0+i; 
w = 16; 
p = 500;

JuliaSet(c,w,p);
title('c = 0+i, filament dendrites.')
pause(1)

clear all

c = -0.125+0.i; 
w = 16; 
p = 500;

JuliaSet(c,w,p);
title('c = -0.125+0.i, parabolic case.')
pause(1)

clear all

c = 0.11031-0.67037i; 
w = 20; 
p = 700;

JuliaSet(c,w,p);
title('c = 0.11031-0.67037i. Dust of Fatou.')
pause(1)

clear all

% Other parabolic case
c = 0.27334+0.00742i; 
w = 16; 
p = 500;

JuliaSet(c,w,p);
title('0.27334+0.00742i.')


% "Complex and Chaotic Nonlinear Dynamics. 
% Advances in Economics and Finance, 
% Mathematics and Statistics" 
% T.Vialar, Springer 2009.
% Copyright(c).