How could I speed up this loop made to plot using the pcolor function?

Hi, i am trying to plot by pcolor a matrix called H, of size 7x7x20001, where 20001 represents a temporary dimension. My idea is to be able to see the evolution of the graph in each step, that is why I have used the drawnow. The problem is that the complete realization of the loop is slow, therefore, the speed of evolution of the plotted graph is also slow. By the way, I am not sure that the use of the pcolor function in itself is fast or not, therefore I do not know if the origin of my problem is really in how I did the loop or in how the pcolor function works.

set(groot,'defaultAxesFontSize',12)
for k=1:length(r)
 pcolor(H(:,:,k)); 
 colorbar;
 xlabel('Dirección x');
 ylabel('Dirección y');
 colormap(jet); 
 caxis([23 100]);
 shading interp;
 axis equal tight; 
 drawnow
end

Answer 1

[KSSV]

A = rand(7,7,20001) ;
[m,n,p] = size(A) ;
[X,Y] = meshgrid(1:n,1:m) ;

set(groot,'defaultAxesFontSize',12)
h = pcolor(A(:,:,1)) ;
colorbar;
xlabel('Dirección x');
ylabel('Dirección y');
colormap(jet);
% caxis([23 100]);
shading interp;
axis equal tight;
for i = 1:p
    set(h,'CData',A(:,:,i)) ;
    title(num2str(i))
    drawnow
end

[Luca Ribba Esteva]

Wow bro, with your code its a lot more faster, you are the mvp, thx you very much :)