N =17;
data = [1:N+1;1:N+1];

subplot(1,3,1)
colormap(hsv(N))
pcolor(data)
set(gca,'XtickLAbel','')
title('Figura 1: Auto Limites')
caxis auto
colorbar

subplot(1,3,2)
pcolor(data)
axis off
title('Figura 2: Limites Extendidos')
caxis([-5,N+5])
colorbar

subplot(1,3,3)
pcolor(data)
axis off
title('Figura 3: Limites Restringidos')
caxis([5,N-5])
colorbar

pause
figure

map = jet(2000);
colormap(map);
subplot(1,2,1);
pcolor(R1,Z1,real(f1))
set(gca,'XtickLabel','')
title('Figura 1: Instante 0.01')
%caxis auto
caxis([0,1])
colorbar

subplot(1,2,2);
pcolor(R2,Z2,real(f2))
title('Figura 2: Instante 1.5');
cmin = min(min(real(f2)));
cmax = max(max(real(f2)));
caxis([0,1])
colorbar

%subplot(1,2,2);
%pcolor(R3,Z3,real(f3))
%title('Figura 2: Instante 5.0');
%cmin = min(min(real(f3)));
%cmax = max(max(real(f3)));
%caxis([cmin,cmax])
%colorbar

%subplot(1,2,2);
%pcolor(R4,Z4,real(f4))
%title('Figura 2: Instante 14.0');
%cmin = min(min(real(f4)));
%cmax = max(max(real(f4)));
%caxis([cmin,cmax])
%colorbar

%subplot(1,2,2);
%pcolor(R5,Z5,real(f5))
%title('Figura 2: Instante 50.0');
%cmin = min(min(real(f5)));
%cmax = max(max(real(f5)));
%caxis([0,1])
%colorbar

%subplot(1,2,2);
%pcolor(R6,Z6,real(f6))
%title('Figura 2: Instante 0.3');
%cmin = min(min(real(f6)));
%cmax = max(max(real(f6)));
%caxis([0,1])
%colorbar