I then create a matrix for Omega and a vector for Xi, give it the appropriate dimensions, and subsequently I introduce the constraint that the initial positions have to be world-size/2.0, with a strength value of 1.0, 