library(lpSolve) f.obj <- c(6, 3, 4) f.con <- matrix(c(4,3,6, 1,0,0, 0,1,0, 0,0,1), ncol=3, byrow=TRUE) # 制約式の左辺の係数 2行目から4行目は変数が1以下の制約 f.dir <- c("<=", "<=", "<=", "<=") # 制約式の等号・不等号 f.rhs <- c(9, 1, 1, 1) # 制約式の右辺 res <- lp("max", f.obj, f.con, f.dir, f.rhs, int.vec=1:3) # 整数計画 #res <- lp("max", f.obj, f.con, f.dir, f.rhs) # 線形計画