BdzHmgnS=[4 9 19 21 24 27 28 29 36 37 38 49 53 55 57 69 77 78 79 81 ...
  84 93 105 106];
BdzHmgnM=[3 15 18 22 35 54 67 68 76 80 92 103 104];
BdzHmgnL=[5 6 20 25 32 33 34 52 61 66 71 73 82 87 72 100 101 102 110 111 112];
BdzHmgn=sort([BdzHmgnS BdzHmgnM BdzHmgnL]);
BdzSemiHmgn=[1 7 8 12 13 16 17 23 26 30 31 51 56 60 64 65 74 75 ...
    83 85 94 95 98 109];
BdzNonHmgn=[2 10 11 39 40 41 42 43 44 45 46 47 48 50 58 59 62 63 ...
    70 86 88 89 90 91 96 97 99 107 108];
BdzAll=sort([BdzHmgn BdzSemiHmgn BdzNonHmgn]);

for i=1:length(BdzHmgnS)
  BdzHmgnSIdx(i)=find(BdzHmgnS(i)==BdzAll);
end
for i=1:length(BdzHmgnM)
  BdzHmgnMIdx(i)=find(BdzHmgnM(i)==BdzAll);
end
for i=1:length(BdzHmgnL)
  BdzHmgnLIdx(i)=find(BdzHmgnL(i)==BdzAll);
end
for i=1:length(BdzHmgn)
  BdzHmgnIdx(i)=find(BdzHmgn(i)==BdzAll);
end
for i=1:length(BdzSemiHmgn)
  BdzSemiHmgnIdx(i)=find(BdzSemiHmgn(i)==BdzAll);
end
for i=1:length(BdzNonHmgn)
  BdzNonHmgnIdx(i)=find(BdzNonHmgn(i)==BdzAll);
end
