def calculate(self, T):
return Zabransky_quasi_polynomial(T, self.Tc, *self.coeffs)<EOL>
r'''Method to actually calculate heat capacity as a function of 
        temperature.

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

        Returns
        -------
        Cp : float
            Liquid heat capacity as T, [J/mol/K]