plot 1 f^3/(e^(1 f) -1) from 0 to 20