squarewave Cosine-like bipolar square wave.
y = squarewave(x, duty), where
x : angle in radians
y : duty cycle in [0, 1]
y = squarewave(x) to use a duty cycle of 1