ABSTRACTWe explore a function with two shape parameters for the dark-matter halo density profile subject to baryonic effects, which is a special case of the general Zhao family of models applied to simulated dark-matter haloes by Dekel et al. This profile has variable inner slope and concentration parameter, and analytic expressions for the gravitational potential, velocity dispersion, and lensing properties. Using the Numerical Investigation of a Hundred Astrophysical Objects cosmological simulations, we find that it provides better fits than the Einasto profile and the generalized NFW profile with variable inner slope, in particular towards the halo centres. We show that the profile parameters are correlated with the stellar-to-halo mass ratio Mstar/Mvir. This defines a mass-dependent density profile describing the average dark-matter profiles in all galaxies, which can be directly applied to observed rotation curves of galaxies, gravitational lenses, and semi-analytic models of galaxy formation or satellite–galaxy evolution. The effect of baryons manifests itself by a significant flattening of the inner density slope and a 20 per cent decrease of the concentration parameter for Mstar/Mvir = 10−3.5–10−2, corresponding to $M_{\rm star} \!∼ \! 10^{7-10}\, \mathrm{ M}_⊙$. The accuracy by which this profile fits simulated galaxies is similar to certain multiparameter mass-dependent profiles, but its fewer parameters and analytic nature make it most desirable for many purposes.