Bernstein Functions with p0 = 1e-6

This is done because the Bernstein is under-constrained.  I'm not sure if it changes the best fit much