mod coeffs.nb in StatThermalEq