Add quadratic term option

Add one more param: square.term.
If T, then return both interaction and quadratic term;
If F, then return only interaction term;
If 'only' (or anything except T or F), then return only quadratic term.