Added Pella-Tomlinson form. Still need to fix Binf and shape of production curve in plot
Implemented the stochastic reference points and equilibrium expression from Bordet and Rivest 2014 and now using these as defualts