This function generate an incidence matrix for given v,b,k and concurrence matrix
ibdgen(v,b,k,NNPo,ntrial,pbar)number of treatments
number of blocks
block size
specified concurrence matrix
number of trials, default is 5
logical, progress bar to be displayed. Default is FALSE
returns a v by b incidence matrix or a text that design was not found