@staticmethod<EOL><INDENT>def Mathias(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*(<NUM_LIT:1> + c1*(<NUM_LIT:1>-sqrt(Tr)) -c2*(<NUM_LIT:1>-Tr)*(<NUM_LIT>-Tr))**<NUM_LIT:2><EOL>if not full:<EOL><INDENT>return a_alpha<EOL><DEDENT>else:<EOL><INDENT>da_alpha_dT = a*(c1*(-sqrt(T/Tc) + <NUM_LIT:1>) - c2*(-T/Tc + <NUM_LIT>)*(-T/Tc + <NUM_LIT:1>) + <NUM_LIT:1>)*(<NUM_LIT:2>*c2*(-T/Tc + <NUM_LIT>)/Tc + <NUM_LIT:2>*c2*(-T/Tc + <NUM_LIT:1>)/Tc - c1*sqrt(T/Tc)/T)<EOL>d2a_alpha_dT2 = a*((<NUM_LIT:8>*c2/Tc**<NUM_LIT:2> - c1*sqrt(T/Tc)/T**<NUM_LIT:2>)*(c1*(sqrt(T/Tc) - <NUM_LIT:1>) + c2*(T/Tc - <NUM_LIT:1>)*(T/Tc - <NUM_LIT>) - <NUM_LIT:1>) + (<NUM_LIT:2>*c2*(T/Tc - <NUM_LIT:1>)/Tc + <NUM_LIT:2>*c2*(T/Tc - <NUM_LIT>)/Tc + c1*sqrt(T/Tc)/T)**<NUM_LIT:2>)/<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 Mathias (1983) [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 = \left(c_{1} \left(- \sqrt{\frac{T}{Tc}} + 1\right)
            - c_{2} \left(- \frac{T}{Tc} + 0.7\right) \left(- \frac{T}{Tc} 
            + 1\right) + 1\right)^{2}

        References
        ----------
        .. [1] Mathias, Paul M. "A Versatile Phase Equilibrium Equation of 
           State." Industrial & Engineering Chemistry Process Design and 
           Development 22, no. 3 (July 1, 1983): 385-91. 
           doi:10.1021/i200022a008.