  We present a procedure for computing the convolution of exponential signals
without the need of solving integrals or summations. The procedure requires the
resolution of a system of linear equations involving Vandermonde matrices. We
apply the method to solve ordinary differential/difference equations with
constant coefficients.
