# Set up BIBD  or Exist a BIBD?
> require(crossdes)
Loading required package: crossdes
Loading required package: AlgDesign
Loading required package: gtools

#problem 4.46  
> design=find.BIB(8,14,4)  #a=8, b=14, k=4. trt level: 1,2,3...8
> design
      [,1] [,2] [,3] [,4]
 [1,]    1    4    5    8
 [2,]    1    3    5    6
 [3,]    3    4    6    8
 [4,]    1    2    3    8
 [5,]    2    4    7    8
 [6,]    3    5    7    8
 [7,]    1    6    7    8
 [8,]    2    3    4    5
 [9,]    1    2    4    6
[10,]    1    2    3    7
[11,]    3    4    6    7
[12,]    2    5    6    7
[13,]    1    4    5    7
[14,]    2    5    6    8
> isGYD(find.BIB(8,14,4))

[1] The design is a balanced incomplete block design w.r.t. rows.

#problem 4.49
> isGYD(find.BIB(8,16,4))

[1] The design is neither balanced w.r.t. rows nor w.r.t. columns.