#####
##### CW(21,16) Dan Gordon 
##### P = 1, 2, 5, 7, 9, 10, 11, 12, 13, 18
##### N = 0, 3, 8, 14, 15, 17
#####

> P := [1, 2, 5, 7, 9, 10, 11, 12, 13, 18]:
> N := [0, 3, 8, 14, 15, 17]:
> W_21_16 := CirculantMatrixPN(21,16,P,N):

> writeto("W23.16.txt");