SP function

Surplus production model with FMSY and MSY as leading parameters