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, multiply the result by the spacing increment.