def COSTALD(T, Tc, Vc, omega):
Tr = T/Tc<EOL>V_delta = (-<NUM_LIT> + <NUM_LIT>*Tr - <NUM_LIT>*Tr**<NUM_LIT:2><EOL>- <NUM_LIT>*Tr**<NUM_LIT:3>)/(Tr - <NUM_LIT>)<EOL>V_0 = <NUM_LIT:1> - <NUM_LIT>*(<NUM_LIT:1>-Tr)**(<NUM_LIT:1>/<NUM_LIT>) + <NUM_LIT>*(<NUM_LIT:1>-Tr)**(<NUM_LIT:2>/<NUM_LIT>)- <NUM_LIT>*(<NUM_LIT:1>-Tr) + <NUM_LIT>*(<NUM_LIT:1>-Tr)**(<NUM_LIT:4>/<NUM_LIT>)<EOL>return Vc*V_0*(<NUM_LIT:1>-omega*V_delta)<EOL>
r'''Calculate saturation liquid density using the COSTALD CSP method.

    A popular and accurate estimation method. If possible, fit parameters are
    used; alternatively critical properties work well.

    The density of a liquid is given by:

    .. math::
        V_s=V^*V^{(0)}[1-\omega_{SRK}V^{(\delta)}]

        V^{(0)}=1-1.52816(1-T_r)^{1/3}+1.43907(1-T_r)^{2/3}
        - 0.81446(1-T_r)+0.190454(1-T_r)^{4/3}

        V^{(\delta)}=\frac{-0.296123+0.386914T_r-0.0427258T_r^2-0.0480645T_r^3}
        {T_r-1.00001}

    Units are that of critical or fit constant volume.

    Parameters
    ----------
    T : float
        Temperature of fluid [K]
    Tc : float
        Critical temperature of fluid [K]
    Vc : float
        Critical volume of fluid [m^3/mol].
        This parameter is alternatively a fit parameter
    omega : float
        (ideally SRK) Acentric factor for fluid, [-]
        This parameter is alternatively a fit parameter.

    Returns
    -------
    Vs : float
        Saturation liquid volume

    Notes
    -----
    196 constants are fit to this function in [1]_.
    Range: 0.25 < Tr < 0.95, often said to be to 1.0

    This function has been checked with the API handbook example problem.

    Examples
    --------
    Propane, from an example in the API Handbook

    >>> Vm_to_rho(COSTALD(272.03889, 369.83333, 0.20008161E-3, 0.1532), 44.097)
    530.3009967969841


    References
    ----------
    .. [1] Hankinson, Risdon W., and George H. Thomson. "A New Correlation for
       Saturated Densities of Liquids and Their Mixtures." AIChE Journal
       25, no. 4 (1979): 653-663. doi:10.1002/aic.690250412