@staticmethod<EOL><INDENT>def Melhem(self, T, full=True, quick=True):<DEDENT>
c1, c2 = self.alpha_function_coeffs<EOL>T, Tc, a = self.T, self.Tc, self.a<EOL>a_alpha = a*exp(c1*(-T/Tc + <NUM_LIT:1>) + c2*(-sqrt(T/Tc) + <NUM_LIT:1>)**<NUM_LIT:2>)<EOL>if not full:<EOL><INDENT>return a_alpha<EOL><DEDENT>else:<EOL><INDENT>da_alpha_dT = a*((-c1/Tc - c2*sqrt(T/Tc)*(-sqrt(T/Tc) + <NUM_LIT:1>)/T)*exp(c1*(-T/Tc + <NUM_LIT:1>) + c2*(-sqrt(T/Tc) + <NUM_LIT:1>)**<NUM_LIT:2>))<EOL>d2a_alpha_dT2 = a*(((c1/Tc - c2*sqrt(T/Tc)*(sqrt(T/Tc) - <NUM_LIT:1>)/T)**<NUM_LIT:2> + c2*(<NUM_LIT:1>/Tc - sqrt(T/Tc)*(sqrt(T/Tc) - <NUM_LIT:1>)/T)/(<NUM_LIT:2>*T))*exp(-c1*(T/Tc - <NUM_LIT:1>) + c2*(sqrt(T/Tc) - <NUM_LIT:1>)**<NUM_LIT:2>))<EOL>return a_alpha, da_alpha_dT, d2a_alpha_dT2<EOL><DEDENT>
r'''Method to calculate `a_alpha` and its first and second
        derivatives according to Melhem et al. (1989) [1]_. Returns `a_alpha`, 
        `da_alpha_dT`, and `d2a_alpha_dT2`. See `GCEOS.a_alpha_and_derivatives`  
        for more documentation. Two coefficients needed.

        .. math::
            \alpha = e^{c_{1} \left(- \frac{T}{Tc} + 1\right) + c_{2} 
            \left(- \sqrt{\frac{T}{Tc}} + 1\right)^{2}}

        References
        ----------
        .. [1] Melhem, Georges A., Riju Saini, and Bernard M. Goodwin. "A 
           Modified Peng-Robinson Equation of State." Fluid Phase Equilibria 
           47, no. 2 (August 1, 1989): 189-237. 
           doi:10.1016/0378-3812(89)80176-1.