scipy.optimize.leastsq con restricciones vinculadas

Estoy buscando una rutina de optimización dentro de scipy / numpy que pueda resolver un problema de tipo de mínimos cuadrados no lineales (por ejemplo, ajustar una función paramétrica a un gran conjunto de datos) pero que incluya límites y ...