Hi, Looking again at the help text ('help binary matrix'), there was no bug in the first place -- I thought that for the non-uniform matrix the columns were numbered y1 to yN, so with N columns -- but actually they are numbered y0 to yN -- so N+1 columns. so the number N+1 is correct.