% visualize sections previously saved
% Copyright Radu Dogaru, 1995, 2006, 
% e-mail: radu_d@ieee.org
%-------------------------------------------------------
% Appendix to Book:  
% R. Dogaru: "Universality and Emergent Computation in CNNs", World Scientific 2003
% http://www.worldscibooks.com/chaos/5053.html
%----------------------------------------------------------------------------------- 

function mvis_sec
load sec0 

[N N]=size(ed1);
amin=a(1); amax=a(2);
bmin=b(1); bmax=b(2);

clf
set(figure(1),'Position',[543 383 717 333]);
map=[0 0 1; 1 0 0; 0 0 0; 0 0 0; 0 .7 0; 1 1 1; 1 .7 0; 0 0 0; 0 1 1; 1 1 0; 0 0 0; 0 0 0; 0 1 0];
%blue   red  black  black  dark_green white orange1  black  yellow orange2 black  black  green 
%blue = passive 
%green (all 3 equilibria are unstable), dark green (one root and unstable)
%red, orange1, orange2, yellow : edge of chaos - at least one equilibria is stable, 
%white 3 equilibria: 1-passive, 1-stable & active and 1-unstable & active
%---------------------------------------------------------------------------


colormap(map);
subplot(1,2,1); image(ed1);axis('xy'); grid on; 
axis image; xlabel('a'); ylabel('b');
title(['eps=',num2str(e),' a=[',num2str(a),']  b=[',num2str(b),']  c=',num2str(c)]);
% prepare a matrix for color legend 
M=[0 1 2 3 4 5 6 7 8 9 10 11 12];
M=[M; M]+1;
subplot(1,2,2); image(M'); axis image; grid on;