How to fit a polynomial with some of the coefficients constrained? How to fit a polynomial with some of the coefficients constrained? numpy numpy