I want to optimize the function with inequality constraints in Stata-mata, how should I do it.