We extend for the first time the applicability of the optimal homotopy asymptotic method (OHAM) to find the algorithm of approximate analytic solution of delay differential equations (DDEs). The analytical solutions for various examples of linear and nonlinear and system of initial value problems of DDEs are obtained successfully by this method. However, this approach does not depend on small or large parameters in comparison to other perturbation methods. This method provides us with a convenient way to control the convergence of approximation series. The results which are obtained revealed that the proposed method is explicit, effective, and easy to use.