def Dadgostar_Shaw_integral(T, similarity_variable):
a = similarity_variable<EOL>a2 = a*a<EOL>T2 = T*T<EOL>a11 = -<NUM_LIT><EOL>a12 = <NUM_LIT><EOL>a21 = <NUM_LIT><EOL>a22 = -<NUM_LIT><EOL>a31 = -<NUM_LIT><EOL>a32 = <NUM_LIT><EOL>constant = <NUM_LIT><EOL>H = T2*T/<NUM_LIT>*(a2*a32 + a*a31) + T2*<NUM_LIT:0.5>*(a2*a22 + a*a21) + T*constant*(a2*a12 + a*a11)<EOL>return H*<NUM_LIT><EOL>
r'''Calculate the integral of liquid constant-pressure heat capacitiy 
    with the similarity variable concept and method as shown in [1]_.

    Parameters
    ----------
    T : float
        Temperature of gas [K]
    similarity_variable : float
        similarity variable as defined in [1]_, [mol/g]

    Returns
    -------
    H : float
        Difference in enthalpy from 0 K, [J/kg]

    Notes
    -----
    Original model is in terms of J/g/K. Note that the model is for predicting
    mass heat capacity, not molar heat capacity like most other methods!
    Integral was computed with SymPy.

    See Also
    --------
    Dadgostar_Shaw
    Dadgostar_Shaw_integral_over_T

    Examples
    --------
    >>> Dadgostar_Shaw_integral(300.0, 0.1333)
    238908.15142664989

    References
    ----------
    .. [1] Dadgostar, Nafiseh, and John M. Shaw. "A Predictive Correlation for
       the Constant-Pressure Specific Heat Capacity of Pure and Ill-Defined
       Liquid Hydrocarbons." Fluid Phase Equilibria 313 (January 15, 2012):
       211-226. doi:10.1016/j.fluid.2011.09.015.