Does there exist a causal nonlinear, time-invariant system that maps each input function $k \cdot\cos(fx)$ to $k \cdot\cos(2fx)$ for all choices of $k, f \in \Bbb R$?
If so, can one represent this system via some sort of differential equation, or in discrete terms, as a nonlinear difference equation?
This would be similar to the second Chebyshev polynomial, but that only works for $k=1$.