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