def solve_T(self, P, V, quick=True):
a, b, Tc, m = self.a, self.b, self.Tc, self.m<EOL>if quick:<EOL><INDENT>x0 = R*Tc<EOL>x1 = V*b<EOL>x2 = x0*x1<EOL>x3 = V*V<EOL>x4 = x0*x3<EOL>x5 = m*m<EOL>x6 = a*x5<EOL>x7 = b*x6<EOL>x8 = V*x6<EOL>x9 = (x2 + x4 + x7 - x8)**<NUM_LIT:2><EOL>x10 = x3*x3<EOL>x11 = R*R*Tc*Tc<EOL>x12 = a*a<EOL>x13 = x5*x5<EOL>x14 = x12*x13<EOL>x15 = b*b<EOL>x16 = x3*V<EOL>x17 = a*x0<EOL>x18 = x17*x5<EOL>x19 = <NUM_LIT>*b*x16<EOL>x20 = -<NUM_LIT>*V*b*x14 + <NUM_LIT>*V*x15*x18 + x10*x11 + x11*x15*x3 + x11*x19 + x14*x15 + x14*x3 - <NUM_LIT:2>*x16*x18<EOL>x21 = V - b<EOL>x22 = <NUM_LIT:2>*m*x17<EOL>x23 = P*x4<EOL>x24 = P*x8<EOL>x25 = x1*x17<EOL>x26 = P*R*Tc<EOL>x27 = x17*x3<EOL>x28 = V*x12<EOL>x29 = <NUM_LIT>*m*m*m<EOL>x30 = b*x12<EOL>return -Tc*(<NUM_LIT>*a*m*x9*(V*x21*x21*x21*(V + b)*(P*x2 + P*x7 + x17 + x18 + x22 + x23 - x24))**<NUM_LIT:0.5>*(m + <NUM_LIT:1.>) - x20*x21*(-P*x16*x6 + x1*x22 + x10*x26 + x13*x28 - x13*x30 + x15*x23 + x15*x24 + x19*x26 + x22*x3 + x25*x5 + x25 + x27*x5 + x27 + x28*x29 + x28*x5 - x29*x30 - x30*x5))/(x20*x9)<EOL><DEDENT>else:<EOL><INDENT>return Tc*(-<NUM_LIT:2>*a*m*sqrt(V*(V - b)**<NUM_LIT:3>*(V + b)*(P*R*Tc*V**<NUM_LIT:2> + P*R*Tc*V*b - P*V*a*m**<NUM_LIT:2> + P*a*b*m**<NUM_LIT:2> + R*Tc*a*m**<NUM_LIT:2> + <NUM_LIT:2>*R*Tc*a*m + R*Tc*a))*(m + <NUM_LIT:1>)*(R*Tc*V**<NUM_LIT:2> + R*Tc*V*b - V*a*m**<NUM_LIT:2> + a*b*m**<NUM_LIT:2>)**<NUM_LIT:2> + (V - b)*(R**<NUM_LIT:2>*Tc**<NUM_LIT:2>*V**<NUM_LIT:4> + <NUM_LIT:2>*R**<NUM_LIT:2>*Tc**<NUM_LIT:2>*V**<NUM_LIT:3>*b + R**<NUM_LIT:2>*Tc**<NUM_LIT:2>*V**<NUM_LIT:2>*b**<NUM_LIT:2> - <NUM_LIT:2>*R*Tc*V**<NUM_LIT:3>*a*m**<NUM_LIT:2> + <NUM_LIT:2>*R*Tc*V*a*b**<NUM_LIT:2>*m**<NUM_LIT:2> + V**<NUM_LIT:2>*a**<NUM_LIT:2>*m**<NUM_LIT:4> - <NUM_LIT:2>*V*a**<NUM_LIT:2>*b*m**<NUM_LIT:4> + a**<NUM_LIT:2>*b**<NUM_LIT:2>*m**<NUM_LIT:4>)*(P*R*Tc*V**<NUM_LIT:4> + <NUM_LIT:2>*P*R*Tc*V**<NUM_LIT:3>*b + P*R*Tc*V**<NUM_LIT:2>*b**<NUM_LIT:2> - P*V**<NUM_LIT:3>*a*m**<NUM_LIT:2> + P*V*a*b**<NUM_LIT:2>*m**<NUM_LIT:2> + R*Tc*V**<NUM_LIT:2>*a*m**<NUM_LIT:2> + <NUM_LIT:2>*R*Tc*V**<NUM_LIT:2>*a*m + R*Tc*V**<NUM_LIT:2>*a + R*Tc*V*a*b*m**<NUM_LIT:2> + <NUM_LIT:2>*R*Tc*V*a*b*m + R*Tc*V*a*b + V*a**<NUM_LIT:2>*m**<NUM_LIT:4> + <NUM_LIT:2>*V*a**<NUM_LIT:2>*m**<NUM_LIT:3> + V*a**<NUM_LIT:2>*m**<NUM_LIT:2> - a**<NUM_LIT:2>*b*m**<NUM_LIT:4> - <NUM_LIT:2>*a**<NUM_LIT:2>*b*m**<NUM_LIT:3> - a**<NUM_LIT:2>*b*m**<NUM_LIT:2>))/((R*Tc*V**<NUM_LIT:2> + R*Tc*V*b - V*a*m**<NUM_LIT:2> + a*b*m**<NUM_LIT:2>)**<NUM_LIT:2>*(R**<NUM_LIT:2>*Tc**<NUM_LIT:2>*V**<NUM_LIT:4> + <NUM_LIT:2>*R**<NUM_LIT:2>*Tc**<NUM_LIT:2>*V**<NUM_LIT:3>*b + R**<NUM_LIT:2>*Tc**<NUM_LIT:2>*V**<NUM_LIT:2>*b**<NUM_LIT:2> - <NUM_LIT:2>*R*Tc*V**<NUM_LIT:3>*a*m**<NUM_LIT:2> + <NUM_LIT:2>*R*Tc*V*a*b**<NUM_LIT:2>*m**<NUM_LIT:2> + V**<NUM_LIT:2>*a**<NUM_LIT:2>*m**<NUM_LIT:4> - <NUM_LIT:2>*V*a**<NUM_LIT:2>*b*m**<NUM_LIT:4> + a**<NUM_LIT:2>*b**<NUM_LIT:2>*m**<NUM_LIT:4>))<EOL><DEDENT>
r'''Method to calculate `T` from a specified `P` and `V` for the SRK
        EOS. Uses `a`, `b`, and `Tc` obtained from the class's namespace.

        Parameters
        ----------
        P : float
            Pressure, [Pa]
        V : float
            Molar volume, [m^3/mol]
        quick : bool, optional
            Whether to use a SymPy cse-derived expression (3x faster) or 
            individual formulas

        Returns
        -------
        T : float
            Temperature, [K]

        Notes
        -----
        The exact solution can be derived as follows; it is excluded for 
        breviety.

        >>> from sympy import *
        >>> P, T, V, R, a, b, m = symbols('P, T, V, R, a, b, m')
        >>> Tc, Pc, omega = symbols('Tc, Pc, omega')
        >>> a_alpha = a*(1 + m*(1-sqrt(T/Tc)))**2
        >>> SRK = R*T/(V-b) - a_alpha/(V*(V+b)) - P
        >>> # solve(SRK, T)