def Laliberte_heat_capacity(T, ws, CASRNs):
Cp_w = Laliberte_heat_capacity_w(T)<EOL>w_w = <NUM_LIT:1> - sum(ws)<EOL>Cp = w_w*Cp_w<EOL>for i in range(len(CASRNs)):<EOL><INDENT>d = _Laliberte_Heat_Capacity_ParametersDict[CASRNs[i]]<EOL>Cp_i = Laliberte_heat_capacity_i(T, w_w, d["<STR_LIT>"], d["<STR_LIT>"], d["<STR_LIT>"], d["<STR_LIT>"], d["<STR_LIT>"], d["<STR_LIT>"])<EOL>Cp = Cp + ws[i]*Cp_i<EOL><DEDENT>return Cp<EOL>
r'''Calculate the heat capacity of an aqueous electrolyte mixture using the
    form proposed by [1]_.
    Parameters are loaded by the function as needed.

    .. math::
        TODO

    Parameters
    ----------
    T : float
        Temperature of fluid [K]
    ws : array
        Weight fractions of fluid components other than water
    CASRNs : array
        CAS numbers of the fluid components other than water

    Returns
    -------
    Cp : float
        Solution heat capacity, [J/kg/K]

    Notes
    -----
    Temperature range check is not implemented.
    Units are Kelvin and J/kg/K.

    Examples
    --------
    >>> Laliberte_heat_capacity(273.15+1.5, [0.00398447], ['7647-14-5']) 
    4186.569908672113

    References
    ----------
    .. [1] Laliberte, Marc. "A Model for Calculating the Heat Capacity of
       Aqueous Solutions, with Updated Density and Viscosity Data." Journal of
       Chemical & Engineering Data 54, no. 6 (June 11, 2009): 1725-60.
       doi:10.1021/je8008123