На одном листе построить 3 графика

clear all
file_d = dir('d*.txt');
file_h = dir('H*.txt');
file_c = dir('Z*.txt');
 
figure
set(gcf,'Position',[1 1 1600 900])
 
for i = 1:2
%  открываем файлы и считываем информацию     
    filename_d = {file_d.name}; 
    data_d = importdata(char(filename_d(i)));
    filename_h = {file_h.name};
    data_h = importdata(char(filename_h(i)));
    filename_c = {file_c.name};
    data_c = importdata(char(filename_c(i)));   
%  создаем линейную дву-мерную картезианскую сетку  
    x.YLim = [0,13];
    x = gca; 
%  выбираем цвет
    colormap(jet(100))
%  строим 1 график     
    contourf(linspace(0,150,length(data_d(1,:))), linspace(9,12,length(data_d(:,1))), data_d, 200, 'EdgeColor', 'none')
    line([0 150],[12 12], 'LineWidth', 3, 'color', 'black')
    line([0 150],[9 9],   'LineWidth', 3, 'color', 'black')
    colorbar
%  строим 2 график     
    contourf(linspace(0,150,length(data_h(1,:))), linspace(5,8 ,length(data_h(:,1))), data_h, 200, 'EdgeColor', 'none')
    line([0 150],[8 8],  'LineWidth',3,'color','black')
    line([0 150],[5 5],  'LineWidth',3,'color','black')
    colorbar
%  строим 3 график     
    contourf(linspace(0,150,length(data_c(1,:))), linspace(1,4 ,length(data_c(:,1))), data_c, 200, 'EdgeColor', 'none')
    line([0 150],[4 4],  'LineWidth',3,'color','black')
    line([0 150],[1 1],  'LineWidth',3,'color','black')
    colorbar
%    saveas(gcf,[char(filename_h),'.png'])
end
 
clear i x file

clear, clc
%
file_d = dir('d*.txt');
file_h = dir('H*.txt');
file_c = dir('Z*.txt');
%
figure
set(gcf,'Position',[1 1 1600 900])
%
for i = 1:20000
%  открываем файлы и считываем информацию     
    filename_d = {file_d.name}; 
    data_d = importdata(char(filename_d(i)));  
    filename_h = {file_h.name};
    data_h = importdata(char(filename_h(i)));
    filename_c = {file_c.name};
    data_c = importdata(char(filename_c(i)));   
%  строим 1 график 
    subplot(311);
    ax1 = gca;              %  создаем линейную дву-мерную картезианскую сетку    
    ax1.YLim = [0,3];
    ax1.XLim = [0,150];
    contourf(linspace(0, 150,length(data_d(1,:))), linspace(0, 3,length(data_d(:,1))), data_d, 128, 'EdgeColor', 'none')
    line([0 150],[0 0], 'LineWidth', 3, 'color', 'black')
    line([0 150],[3 3], 'LineWidth', 3, 'color', 'black')
    colormap(ax1,hot(128))   %  выбираем цвет
    colorbar
    title('Density evolution') 
%  строим 2 график     
    subplot(312);
    ax2 = gca;              %  создаем линейную дву-мерную картезианскую сетку    
    ax2.YLim = [0,3];
    ax2.XLim = [0,150];
    contourf(linspace(0, 150,length(data_c(1,:))), linspace(0, 3 ,length(data_c(:,1))), data_c, 128, 'EdgeColor', 'none')
    line([0 150],[0 0],  'LineWidth',3,'color','black')
    line([0 150],[3 3],  'LineWidth',3,'color','black')
    colormap(ax2,hot(128))   %  выбираем цвет
    colorbar
    title('Current evolution') 
%  строим 3 график 
    subplot(313); 
    ax3 = gca;              %  создаем линейную дву-мерную картезианскую сетку    
    ax3.YLim = [0,3];
    ax3.XLim = [0,150];
    contourf(linspace(0, 150,length(data_h(1,:))), linspace(0, 3 ,length(data_h(:,1))), data_h, 128, 'EdgeColor', 'none')
    line([0 150],[0 0],  'LineWidth',3,'color','black')
    line([0 150],[3 3],  'LineWidth',3,'color','black')
    colormap(ax3,hot(128))   %  выбираем цвет
    colorbar
    title('M-field evolution') 
    saveas(gcf,[char(filename_d(i)),'.bmp'])
end
 
clear i x file

for i = 1:20000
   tic % засекаем время
     % считаем/рисуем 
   T(i) = toc; % конец отсчета времени, сохраняем значение
end