The paper deals with the computation of delay Lyapunov matrices for linear exponentially stable time-invariant systems with several pointwise delays. It is shown that the matrices depend continuously on the system delays. Based on this theoretical result, an algorithm of approximation of the delay Lyapunov matrices with a prescribed accuracy is presented. The result is illustrated by an example.