option ps=60 ps=70; 
PROC IML;                   

egg1 = { 1  0  1  0  1,                
         0  1  1  1  0,
         1  0  1  0  1,
         1  0  1  1  0,
         0  1  1  0  0,
         1  0  1  0  1,
         1  0  0  1  1};
                                     

egg2 = { 1  0  1  0  1,
         1  1  1  1  0,
         0  1  1  1  1,
         1  0  1  1  1,
         0  1  1  0  0,
         1  1  1  0  1,
         1  1  1  1  0};

egg0 = { 1  1  1  0  1,            
         1  0  1  0  0,
         0  1  1  1  0,
         1  1  0  0  1,
         0  0  1  0  0,
         1  1  1  0  1,
         1  1  1  1  0};

trsegg1  = egg1`;                   
                                       

egg12    = egg1 + egg2;             
print      trsegg1, egg12;

rowi     = { 1 1 1 1 1 1 1 };        

stage0   = rowi*egg0;                 
stage1   = rowi*egg1;                 
stage2   = rowi*egg2;                 
eggn     = stage1 + stage2;          

rowij    = { 1 0 1 0 1 0 1,           
             0 1 0 1 0 1 0};          
                                              

eggij    = rowij*egg12;          

print      eggij, stage0, stage1,
                  stage2, eggn;

diffij   = egg0 - egg1;               
diff     = rowi * diffij* {1 1 1 1 1}`; 
print      diffij, diff;
run;