corrected (q=e^{-s/3} in (0,1]): q**3/18 - q**2/12 + t**3 - 3*t**2 + t*(q**2/4 + 2) - 2/9
printed   (q=e^{+s/3} in [1,oo)): q**3/18 - q**2/12 + t**3 - 3*t**2 + t*(q**2/4 + 2) - 2/9
