def c_DuttonMaccio(z, m, h=h):
z, m = _check_inputs(z, m)<EOL>a = <NUM_LIT> + <NUM_LIT> * np.exp(-<NUM_LIT> * (z**<NUM_LIT>))  <EOL>b = -<NUM_LIT> + <NUM_LIT> * z                         <EOL>logc200 = a + b * np.log10(m * h / (<NUM_LIT>**<NUM_LIT:12>))  <EOL>concentration = <NUM_LIT>**logc200<EOL>return concentration<EOL>
Concentration from c(M) relation in Dutton & Maccio (2014).

    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
    ----------
    Calculation from Planck-based results of simulations presented in:

    A.A. Dutton & A.V. Maccio, "Cold dark matter haloes in the Planck era:
    evolution of structural parameters for Einasto and NFW profiles,"
    Monthly Notices of the Royal Astronomical Society, Volume 441, Issue 4,
    p.3359-3374, 2014.