nbmod function

The Nicholson-Bailey model