The use of distributions leads naturally to the convolution algebras of linear time-invariant systems and the ones suitable for weakly nonlinear systems emerge as simple extensions to higher order distributions, without having to resort to ad hoc operators.
