Numerical Solution of Delay DEs