How to find the global optimization for a non-convex optimization problem in Pyhon?

Question

I have a portfolio optimization problem. Let\'s say R and P are expected return vector and correlation matrix of N stocks. I want to f