added quadratic constraints, so far does not work....