sum(n*ω*e^(-β*ω*n), n from 0 to +inf)