Rademacher function
y = RADMAC(k,x)
Alternatively, find such that ; where
Psi function Savitzky-Golay smoothing