Is the solution that SQP produces always the best or unique? How can I seek for other solutions?
By varying the start point, $x_0$?
Also, any idea why the start point may be "absurdly" large compared to the solution. I'm trying e.g. $x0=[50000.0,1000.0,3000.0,5000.0]$. When the sol is $[1, 0, 0, 1]$ and at step 2. the algorithm is already around $[1,1,1,1]$.