def c_Prada(z, m, h=h, Om_M=Om_M, Om_L=Om_L):
z, m = _check_inputs(z, m)<EOL>x = (<NUM_LIT:1.> / (<NUM_LIT:1.> + z)) * (Om_L / Om_M)**(<NUM_LIT:1.> / <NUM_LIT>)<EOL>intEQ12 = np.zeros(len(x))  <EOL>for i in range(len(x)):<EOL><INDENT>temp = integrate.quad(lambda v: (v / (<NUM_LIT:1> + v**<NUM_LIT>))**(<NUM_LIT>), <NUM_LIT:0>, x[i])<EOL>intEQ12[i] = temp[<NUM_LIT:0>]<EOL><DEDENT>Da = <NUM_LIT> * ((Om_M / Om_L)**(<NUM_LIT:1.> / <NUM_LIT>)) * (np.sqrt(<NUM_LIT:1.> + x**<NUM_LIT>) /<EOL>(x**(<NUM_LIT>))) * intEQ12<EOL>y = (<NUM_LIT>) / (h * m)<EOL>sigma = Da * (<NUM_LIT> * y**<NUM_LIT>) / (<NUM_LIT:1.> + (<NUM_LIT> * y**<NUM_LIT>) + (<NUM_LIT> * y**<NUM_LIT>))<EOL>c0 = <NUM_LIT><EOL>c1 = <NUM_LIT><EOL>alpha = <NUM_LIT><EOL>x0 = <NUM_LIT><EOL>s0 = <NUM_LIT>  <EOL>s1 = <NUM_LIT>  <EOL>beta = <NUM_LIT><EOL>x1 = <NUM_LIT><EOL>cmin = c0 + (c1 - c0) * ((<NUM_LIT:1.> / np.pi) * np.arctan(alpha * (x - x0)) + <NUM_LIT:0.5>)<EOL>smin = s0 + (s1 - s0) * ((<NUM_LIT:1.> / np.pi) * np.arctan(beta * (x - x1)) + <NUM_LIT:0.5>)<EOL>cmin1393 = c0 + (c1 - c0) * ((<NUM_LIT:1.> / np.pi) * np.arctan(alpha *<EOL>(<NUM_LIT> - x0)) + <NUM_LIT:0.5>)<EOL>smin1393 = s0 + (s1 - s0) * ((<NUM_LIT:1.> / np.pi) * np.arctan(beta *<EOL>(<NUM_LIT> - x1)) + <NUM_LIT:0.5>)<EOL>B0 = cmin / cmin1393<EOL>B1 = smin / smin1393<EOL>sigma_prime = B1 * sigma<EOL>A = <NUM_LIT><EOL>b = <NUM_LIT><EOL>c = <NUM_LIT><EOL>d = <NUM_LIT><EOL>Cs = A * ((sigma_prime / b)**c + <NUM_LIT:1.>) * np.exp(d / (sigma_prime**<NUM_LIT>))<EOL>concentration = B0 * Cs<EOL>return concentration<EOL>
Concentration from c(M) relation published in Prada et al. (2012).

    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.
    Om_M : float, optional
        Matter density parameter. Default is from Planck13.
    Om_L : float, optional
        Cosmological constant density parameter. Default is from Planck13.

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

    Notes
    ----------
    This c(M) relation is somewhat controversial, due to its upturn in
    concentration for high masses (normally we expect concentration to
    decrease with increasing mass). See the reference below for discussion.

    References
    ----------
    Calculation based on results of N-body simulations presented in:

    F. Prada, A.A. Klypin, A.J. Cuesta, J.E. Betancort-Rijo, and J.
    Primack, "Halo concentrations in the standard Lambda cold dark matter
    cosmology," Monthly Notices of the Royal Astronomical Society, Volume
    423, Issue 4, pp. 3018-3030, 2012.