library(lpSolve)

f.obj <- c(4,1,5,0.6,3.5,0.7,4) 

#restriciones
R1 <- c(5,0,2,0,3,1,2)
R2 <- c(3,1,5,0,2,0,1)
R3 <- c(1,0,3,1,2,0,6)
R4 <- c(1,0,3,1,2,0,6)


f.con <- matrix (c(R1,R2,R3,R4), nrow=4, byrow=TRUE)
f.dir <- c("=", "=",">=", "<=")

f.rhs <- c(100,80,120,160)

respuesta <- lp ("min", f.obj, f.con, f.dir, f.rhs ) 

vvv<-respuesta$solution

print(vvv)
respuesta