This code computes the approximation of the convolution integral (*) between two functions, f(t) and g(t), sampled at the times t1,t2,...,tm (unit spacing), via the trapezoidal quadrature rule. To compute the integral with other than unit spacing, mu