Linear time-invariant delay-differential equation systems are approximately represented by a set of linear algebraic equations with the block pulse functions. A least squares estimate is then used to determine the unknown parameters. Examples with satisfactory results are given.