TA Algorithm to find a lean design
Runs n=9, Factors s=4.
Initial Design:
 1 3 1 2
 1 2 3 3
 3 2 1 1
 2 2 2 2
 2 1 1 3
 2 3 3 1
 1 1 2 1
 3 1 3 2
 3 3 2 3
Criterion: maximize Det(X(5)'X(5))*...*Det(X(9)'X(9))
Search Type:(a) - f1(x1) f2(x1), ..., f1(xs) f2(xs)  
Result Design:
 1 3 1 2
 1 2 3 3
 3 2 1 1
 3 1 3 2
 2 1 1 3
 2 3 3 1
 1 1 2 1
 3 3 2 3
 2 2 2 2

Runs Permutation: 1  2  3  8  5  6  7  9  4
Factors Permutation: 1  2  3  4

Max Product Det: 0.13925
============================================
Expanded Matrix:
  1 -1 -1  1 -1 -1 -1  0  2
  1 -1 -1  0  2  1 -1  1 -1
  1  1 -1  0  2 -1 -1 -1 -1
  1  1 -1 -1 -1  1 -1  0  2
  1  0  2 -1 -1 -1 -1  1 -1
  1  0  2  1 -1  1 -1 -1 -1
  1 -1 -1 -1 -1  0  2 -1 -1
  1  1 -1  1 -1  0  2  1 -1
  1  0  2  0  2  0  2  0  2
Det(X'X)^(1/k) -         Std    EachStd   Original
                        ---------------------------
Top  5 Left  5 Corner: 0.415244 0.755082 4.192963
Top  6 Left  6 Corner: 0.648054 0.908560 6.000000
Top  7 Left  7 Corner: 0.624481 0.813921 6.357804
Top  8 Left  8 Corner: 0.828607 0.923999 7.896444
Top  9 Left  9 Corner: 1.000000 1.000000 10.227532
   Det Product: 0.139245745      0.51594291      12917.5989
============================================
Expanded Matrix for initial design:
  1 -1 -1  1 -1 -1 -1  0  2
  1 -1 -1  0  2  1 -1  1 -1
  1  1 -1  0  2 -1 -1 -1 -1
  1  0  2  0  2  0  2  0  2
  1  0  2 -1 -1 -1 -1  1 -1
  1  0  2  1 -1  1 -1 -1 -1
  1 -1 -1 -1 -1  0  2 -1 -1
  1  1 -1 -1 -1  1 -1  0  2
  1  1 -1  1 -1  0  2  1 -1
For inital design
Det(X'X)^(1/k) -         Std    EachStd   Original
                        ---------------------------
Top  5 Left  5 Corner: 0.415244 0.775519 4.192963
Top  6 Left  6 Corner: 0.514361 0.759539 4.762203
Top  7 Left  7 Corner: 0.624481 0.808239 6.357804
Top  8 Left  8 Corner: 0.696772 0.792550 6.640092
Top  9 Left  9 Corner: 1.000000 1.000000 10.227532
   Det Product: 0.092935385     0.37731975      8621.4629