-n - (n+1) *ProductLog[-1, -(n/(n+1)) e^(-n/(n+1))] from 0 to 4