Using pseudospectral methods, a nonlinear delay equation can be approximated with a finite dimensional system of ODEs. I will discuss in which sense the ODE system approximates the dynamical properties of the delay equation, and how the technique can be applied to differential and integro-differential equations, finite and infinite delays, and state-dependent delays. I will present some examples mainly coming from population dynamics.