NB function

The Nicholson-Bailey model