def c_Duffy(z, m, h=h):
z, m = _check_inputs(z, m)<EOL>M_pivot = <NUM_LIT> / h  <EOL>A = <NUM_LIT><EOL>B = -<NUM_LIT><EOL>C = -<NUM_LIT><EOL>concentration = A * ((m / M_pivot)**B) * (<NUM_LIT:1> + z)**C<EOL>return concentration<EOL>
Concentration from c(M) relation published in Duffy et al. (2008).

    Parameters
    ----------
    z : float or array_like
        Redshift(s) of halos.
    m : float or array_like
        Mass(es) of halos (m200 definition), in units of solar masses.
    h : float, optional
        Hubble parameter. Default is from Planck13.

    Returns
    ----------
    ndarray
        Concentration values (c200) for halos.

    References
    ----------
    Results from N-body simulations using WMAP5 cosmology, presented in:

    A.R. Duffy, J. Schaye, S.T. Kay, and C. Dalla Vecchia, "Dark matter
    halo concentrations in the Wilkinson Microwave Anisotropy Probe year 5
    cosmology," Monthly Notices of the Royal Astronomical Society, Volume
    390, Issue 1, pp. L64-L68, 2008.

    This calculation uses the parameters corresponding to the NFW model,
    the '200' halo definition, and the 'full' sample of halos spanning
    z = 0-2. This means the values of fitted parameters (A,B,C) = (5.71,
    -0.084,-0.47) in Table 1 of Duffy et al. (2008).