Update pnbd.R

LL2 to do the log exp sum trick