function [p_shapes, d, r_shape, p_num] = LPmake_Figure9a()

% LPMAKE_FIGURE9A  defines some arrays needed to generate Figure 9a
% using 4 copies each of 8 octominoes to tile a 16 x 24 region with 
% an 8 x 16 notch

% Create binary matrices for the octominoes
for page = 1:8;
    p_shapes(:,:,page) = zeros(16,24);
end
p_shapes(1:4,1:4,1) = [ 1 1 1 0; 1 1 1 0; 1 0 0 0; 1 0 0 0];
p_shapes(1:2,1:5,2) = [1 1 1 1 0; 1 1 0 1 1];
p_shapes(1:4,1:4,3) = [1 0 0 0; 1 0 0 0; 1 0 1 0; 1 1 1 1];
p_shapes(1:3,1:4,4) = [1 1 1 1; 1 0 1 1; 0 0 0 1];
p_shapes(1:4,1:5,5) = [1 0 0 0 0; 1 0 0 0 0; 1 0 0 0 0; 1 1 1 1 1];
p_shapes(1:4,1:3,6) = [1 0 0; 1 1 0; 1 1 0; 1 1 1];
p_shapes(1:3,1:5,7) = [1 1 1 1 1; 0 1 0 0 0; 0 1 1 0 0];
p_shapes(1:4,1:3,8) = [1 0 1; 1 1 1; 1 0 0; 1 1 0];
% Create d and p_num
d = 4*ones(1,8);
p_num = 8;
% Create the binary matrix for the region
r_shape = ones(16,24);
r_shape(9:16,1:16) = 0;
% Create the linear system
[ a, b, var_p ]  = polyomino_multihedral_matrix ( r_shape, p_num,...
    p_shapes, d );
[m,n] = size(a);
% Create and save the LP file
polyomino_lp_write ( 'Figure9a.lp', '\ Figure9a', m, n, a, b );
end