% The program creates the I-beam density distribution
%-----------------------------------------------------------------------------

X = zeros(64,64);
X(11:15,25:54) = ones(5,30);
X(51:55,25:54) = ones(5,30);
X(16:50,38:42) = ones(35,5);
X(28:32,11:15) = ones(5,5);

pcolor(X), colormap(1-gray), shading flat, axis('square'),axis('off')
colorbar
title('Original density distribution')

