Quote:
Originally Posted by patrickjtierney
I came across a clever trick in the previous classes discussions for improving the stability and results of qp() in Octave. Essentially you just add a tiny amount to the H matrix's diagonal (10^15) and then solve for alpha0 using qp. Then you rerun qp using alpha0 as the initial value and use your original H to obtain the value of alpha.

Your trick indeed looks very clever. I find it difficult to believe the huge difference between your solutions, though. So you may want to doublecheck on which one is correct. Hope this helps.