Try this —
M1 = readmatrix('https://www.mathworks.com/matlabcentral/answers/uploaded_files/833435/1to10_30s_15mm_diff_pivoted.txt')
% figure
% hold on
% for k = 1:size(M1,2)
% subplot(12,12,k)
% hf{k} = histfit(M1(:,k)); % Distribution Plot
% [pks{k},locs{k},fwhm{k},prm{k}] = findpeaks(hf{k}(2).YData, hf{k}(2).XData);
% pd{k} = fitdist(M1(:,k), 'Normal'); % Parameters
% grid
% title(sprintf('Column #%3d (\\mu = %5.2f, \\sigma = %5.2f FWHM = %7.2f)',k, pd{k}.mu, pd{k}.sigma,fwhm{k}))
% end
figure
hold on
for k = 1:size(M1,2)
subplot(122,2,k)
hf{k} = histfit(M1(:,k)); % Distribution Plot
[pks{k},locs{k},fwhm{k},prm{k}] = findpeaks(hf{k}(2).YData, hf{k}(2).XData);
pd{k} = fitdist(M1(:,k), 'Normal'); % Parameters
grid
title(sprintf('Column #%3d (\\mu = %5.2f, \\sigma = %5.2f FWHM = %7.2f)',k, pd{k}.mu, pd{k}.sigma,fwhm{k}))
end
hold off
figure
hold on
for k = 1:24:size(M1,2)
subplot(3,2,ceil(k/24))
hf{k} = histfit(M1(:,k)); % Distribution Plot
[pks{k},locs{k},fwhm{k},prm{k}] = findpeaks(hf{k}(2).YData, hf{k}(2).XData);
pd{k} = fitdist(M1(:,k), 'Normal'); % Parameters
grid
title(sprintf('Column #%3d (\\mu = %5.2f, \\sigma = %5.2f FWHM = %7.2f)',k, pd{k}.mu, pd{k}.sigma,fwhm{k}))
end
hold off
.
