Computing the characteristic roots for delay differential equations