library(mco) binh2 <- function (x) { c1 <- 4 * x c2 <- (x - 5) c(crossprod(c1, c1), crossprod(c2, c2)) } binh2.constr <- function (x) { c1 <-