def _calc_gamma(self, V_a):
        a1, a2 = self._calc_strain_coefs()
        fstr_a = self._calc_fstrain(V_a)

        gamma_a = (2*fstr_a+1)*(a1+a2*fstr_a)/(6*(1+a1*fstr_a+0.5*a2*fstr_a**2))

        return gamma_a