int(1/(R+S_0 exp(-beta R/gamma)-N)) dR