Adding a function to solve convex quadratic programs to scipy.optimize

By now, there exists no function tailored to solve convex quadratic programs in the module scipy.optimize. It is strongly suggested to add a function like “quadprog“ in Matlab. Although some other modules like cvxpy can be alternatives, it is still desirable that scipy can make it.