Changed eig  to eigh for symmetric matrices, added function to solve the non-inequality constraints