def _calc_strain_coefs(self):
        V0, gamma0, gammap0 = self.eos_mod.get_param_values(
            param_names=['V0','gamma0','gammap0'])

        a1 = 6*gamma0
        a2 = -12*gamma0 +36*gamma0**2 -18*gammap0
        return a1, a2