Abstract: The spherical NFW profile provides a good approximation to the density distributions of dark matter halos. In dynamical equilibrium, this density profile should correspond to a stationary distribution function of dark matter in phase space. We derive this distribution function and show that it is a good match to those obtained from N-body simulations.