We present a methodology for analyzing robust independent-of-delay and delay-dependent stability of equilibria of systems described by nonlinear Delay Differential Equations by algorithmically constructing appropriate Lyapunov-Krasovskii functionals using the sum of squares decomposition of multivariate polynomials and semidefinite programming. We illustrate the methodology using an example from population dynamics.