integrate p^2/(e^((p-mu)/T)-1)dp