In this work, the problem of computation of the delay Lyapunov matrix is addressed for a scalar retarded-type equation with multiple arbitrary delays, commensurate or not. It is shown that this matrix, which is a function in the scalar case, constitutes a solution of the Fredholm integral equation of the second kind with a piecewise constant kernel. A numerical procedure for its computation is then developed based on the theory of Fredholm equations.