proc factex; factors x1 x2 / nlev = 3; output out=can x1 nvals = (-1 0 1) x2 nvals = (-1 0 1); proc print data=can; run; proc optex data=can seed=57922; model x1 x2 x1*x2 x1*x1 x2*x2; generate n=100 method=fedorov; output out=des; proc print data=des; run;