library(lpSolve)
beneficio=c(15000,12500)
a=c(1.6,1,1,2,1,0,0,1,1,0,0,1)
matriz=matrix(a,nrow=6,ncol=2,byrow=TRUE)
print(matriz)
b=c(40000,33333,12000,8000,22500,15000)
restricciones=c("<=","<=",">=",">=","<=","<=")
resultado=lp("max",beneficio,matriz,restricciones,b)
resultado$objval
resultado$solution