Changed ab and cd parameterizations to use updated guess values from AFTER the fit occurs that crosses the constraints. This seems to work better in terms of coming close to unbounded fit.
Removed some debugging prints.
Now testing for convergence in both AB and CD fits for Cobb-Douglas model. Also, re-running fits if first re-run forces other parameter beyond bounds.
Now including more tests of Cobb-Douglas fits.