library(lpSolve) f.obj <- c(10, 6, 16, 8, 8, 10) # 目的関数の係数 f.con <- matrix (c(1,1,1,0,0,0, 0,0,0,1,1,1, 1,0,0,1,0,0, 0,1,0,0,1,0, 0,0,1,0,0,1), ncol=6, byrow=TRUE) # 制約式の左辺の係数 f.dir <- c("=", "=", "=", "=", "=") # 制約式の等号・不等号 f.rhs <- c(80, 160, 120, 40, 80) # 制約式の右辺 res <- lp("min", f.obj, f.con, f.dir, f.rhs)